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PREFACE 


This  is  Volume  2  of  a  two  volume  report  documenting  the  results  achieved  In  an 
experimental  and  theoretical  program  designed  to  determine  the  effect  of  propellant 
composition.  The  overall  program  objective  is  to  identify,  evaluate,  and  develop 
ammonium  nitrate  propellant  tailoring  techniques.  The  objective  of  this  subtask  of  the 
program  has  been  to  develop  a  model  to  describe  the  combustion  of  ammonium 
nitrate  based  propellants,  including  auxiliary  oxidizers  such  as  ammonium  perchlorate 
and  sodium  nitrate. 


This  volume  of  the  report  describes  the  separate  surface  temperature  model  for 
predicting  combustion  characteristics  of  ammonium  nitrate  solid  propellants.  This 
volume  contains  a  description  of  the  theoretical  development  of  the  model,  and  a 
description  of  the  results.  A  description  of  the  corresponding  computer  program  has 
been  published  and  is  contained  in  AFAL-TR-88-109,  Temperature  Sensitivity  Study, 
Volume  III,  Computer  Users  Manual.  The  test  cases  corresponding  to  ammonium 
nitrate  propellants  are  contained  in  Appendix  A  of  this  report. 


NOTICE 


When  U*S.  Government  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  a  definitely  related  Government  procurement  operation, 
the  fact  that  the  Government  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawi‘.ngs,  specifications,  or  other  data,  is  not  to  be 
regarded  by  implication  or  otherwise,  or  in  any  way  licensing  the  holder  or 
any  other  person  or  corporation,  or  conveying  any  rights  or  permission  to  ^ 
manufacture,  use,  or  sell  any  patented  invention  that  may  be  related  thereto. 
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INTRODUCTION 


When  the  space  shuttle  is  launched,  each  Solid  Rocket  Motor  (SRM)  booster 
generates  approximately  100  tons  of  HCI  in  the  exhaustT  Thus,  with  the  two  boosters, 
over  200  tons  of  HCI  is  generated  from  each  space  shuttle  launch.  The  HCI  presents 
several  potential  problems.  The  HCI  provides  nucleation  sites  for  moisture  in  the  air 
with  the  result  being  the  condensed,  hydrated  form  of  hydrochloric  acid  {HCI.H20). 
The  resultant  acid  cloud  results  in  a  very  visible  cloud  which  through  further 
condensation  can  result  in  acid  rain.  The  acid  rain  is  obviously  detrimental  to  both 
plant  and  animal  life  dependent  on  the  concentration  levels  it  achieves  as  it  reaches 
the  ground.  In  addition,  at  higher  elevations,  the  HCI  can  react  with  the  protective 
layer  of  ozone  causing  a  degradation  of  some  of  the  ozone.  The  HCI  comes  from  the 
AP  (ammonium  perchlorate)  which  is  the  oxidizer  in  the  solid  propellant  for  the  SRM 
boosters.  The  SRM  propellant  contains  approximately  70%  AP"'  which  generates 
slightly  more  than  21  weight  percent  HCI  in  the  exhaust  gases.  For  these  reasons  it  is 
desirable  to  find  a  replacement  for  the  AP  to  eliminate  the  HCL  in  the  rocket  exhaust. 

The  two  oxidizers  used  most  often  in  modern  solid  propellants  are  AP,  and 
cyclotetramethyline  tetranitramine,  HMX.  Although  HMX  does  not  contain  chlorine  and 
generates  "clean"  exhaust,  it  is  relatively  expensive  and  sensitive  which  makes  it 
undesirable  for  a  potential  replacement  for  AP.  To  develop  an  alternate  oxidizer  for 
AP  that  would  be  cheaper  and  would  be  "clean",  i.e.  would  not  produce  HCI  in  the 
exhaust,  ammonium  nitrate,  AN,  and  sodium  nitrate,  SN,  are  being  considered. 
Neither  contain  chlorine  and  both  are  potentially  less  expensive  than  AP.  However, 
neither  are  as  energetic  as  AP  and  the  burning  characteristics  are  not  well  known. 

The  objective  of  this  study  has  been  to  develop  a  model  to  describe  the 
combustion  of  AN  based  propellants,  including  auxiliary  oxidizers  such  as  AP  and  SN. 

No  burn  rate  modeling  research  for  AN  propellants  has  been  conducted  during 
the  last  20  years.  In  the  late  1950s,  AN  propellant  studies  were  performed  by 
researchers  at  Aerojet^-^  that  postulated  a  two-temperature  model  similar  in  concept  to 
the  Beckstead  SST  (Separate  Surface  Temperature)  mode|5.6.  Since  that  time,  the 
lack  of  interest  in  AN  propellants  has  precluded  any  serious  modeling  attempts 
specificaily  related  to  AN. 

The  status  of  combustion  modeling  of  composite  propellants  was  reviewed  by 
Cohen  in  1980^  which  resulted  from  a  1978  JANNAF  Workshop.  In  1982  Cohen^ 
revised  his  modei  to  incorporate  most  of  the  recommendations  that  were  made  during 
the  JANNAF  workshop.  More  recently,  1986,  Gusachenko  and  Zarko^  have  reviewed 
the  status  of  combustion  modeling  of  composite  propellants  including  a  perspective 
from  the  Russian  literature.  Appendix  A  consists  of  an  annotated  bibliography  of 
papers  relevant  to  modeling  monopropellants  and  composite  propeilants  (double 
base  propellant  modeling  papers  were  excluded). 

In  Cohen's  review  article,  he  compared  the  four  models  of  Cohenio-12, 
Beckstead®-''^,  King‘*4-16^  and  the  PEM  model  of  Glick,  Condon,  Renie,  etc."' ^'25  The 
review  concentrated  on  these  four  models  as  being  the  most  representative  of  the 
state  of  the  art  and  as  having  the  most  widely  accepted  physical  basis.  Each  of  these 
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models  is  based  on  the  BDP  multiple  flame  concept26.27  that  is  most  representative  of 
the  physical  mechanisms  that  control  propellant  combustion.  Table  1  is  a  brief 
summary  of  the  contents  and  physical  basis  of  the  various  models  that  have  been 
reviewed  and  considered  for  application  to  the  present  study. 

In  the  development  of  the  various  models  there  has  been  considerable 
discussion  relative  to  the  proper  method  of  averaging  the  burn  rate  of  the  various 
ingredients  in  that  make  up  the  propellant.  Cohen,  the  original  BDP  approach,  and  the 
PEM  models  all  use  an  area  or  surface  averaging  approach  that  nominally  appears  to 
work  well  for  AP  propellants.  Beckstead's  SST  model  and  King's  model  both  use  a 
time  averaging  technique  that  appears  to  work  better  with  HMX  propellants,  and  thus, 
might  be  better  with  mixed  oxidizers  having  widely  varying  burn  rates. 

All  of  the  models  are  based  on  the  BDP  multiple  flame  concept,  considering  a 
premixed  monopropellant  flame,  an  initial  diffusion  flame  between  the  oxidizer  and 
binder  decomposition  products,  and  a  final  diffusion  flame  between  the  products  of  the 
oxidizer  monopropeliant  flame  and  the  binder  decomposition  products.  In  each  the 
analysis  of  the  diffusion  flame  is  based  on  the  Burke-Schumann^s  approach  with 
minor  modifications.  All  are  based  on  global  kinetics.  None  of  the  models  consider 
an  explicit  interaction  of  flames  between  adjacent  particles  of  differing  sizes.  Most  of 
the  models  have  been  applied  to  aluminized  propellants,  but  only  in  a  superficial 
manner. 

The  Beckstead  SST  model  was  selected  as  the  basis  for  development  because 
of  its  flexibility  and  its  apparent  potential  for  handling  mixed  oxidizers  with  significantly 
varying  surface  temperatures.  Therefore,  the  current  work  is  an  extension  of  the 
Beckstead  SST  model  with  some  modifications  taken  from  Cohen's  model^. 
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THERMOCHEMICAL  CALCULATIONS 
General  Performance  Calculations 

A  general  investigation  was  conducted  evaluating  the  performance  of  severai 
alkali  metal  nitrates  as  potential  combustion/scavenger  additives  for  low  cost,  clean 
propellants.  Lithium  nitrate  (LN),  sodium  nitrate  (SN),  and  potassium  nitrate  (KN)  al! 
have  the  potential  to  scavenge  the  undesirable  chloride  from  AP,  yielding  the 
corresponding  alkali  metal  chloride.  In  addition,  KNO3  acts  as  a  phase  stabilizer  for 
AN,  and  NaNOs  has  been  studied  extensiveiy  by  the  Soviets29-3i  in  conjunction  with 
aluminum  agglomeration.  Both  of  these  ingredients  seem  to  show  potential  within  the 
context  of  a  clean  propellant.  Some  of  the  thermochemical  and  physical  properties  of 
these  oxidizers  have  been  tabulated  in  Table  2.  Although  AP,  AN  and  HMX  all  burn 
as  monopropellants,  none  of  the  alkali  metai  nitrates  burn  as  monopropellants  and 
thermochemical  calculations  do  not  give  flame  temperatures,  Isp  values,  etc.  Both  AP 
and  AN  burn  as  monopropellants  providing  excess  oxygen  in  their  products,  29  mole 
%  for  AP  and  14%  for  AN.  HMX  on  the  other  hand  is  actually  fuel  rich  with  CO  being 
the  dominant  carbon  product  with  virtually  no  CO2.  All  of  the  alkali  metai  nitrates  have 
relatively  high  densities  which  is  desirable,  and  AN  has  the  lowest  density  which  is 
undesirabie.  AN  also  has  a  relatively  low  melting  point  and  AP  has  the  highest. 

Table  2  Oxidi2er  Thermochemical  and  Physical  Properties 


AP 

AN 

HMX 

LM 

RN 

m 

Density  (gm/cc) 

1.95 

1.73 

1.91 

2.38 

2.26 

2.11 

Melting  Point  (K) 

865 

443 

553 

534 

580 

607 

"  "  (C) 

592 

170 

280 

261 

307 

334 

Enthalpy  (Kcai/mole) 

-70.7 

-87.4 

18.06 

-115.3 

-111.8 

-118.2 

MW  (gm/mole) 

117.5 

80 

256 

69 

85 

101 

Flame  Temperature  (K) 

1,405 

1,247 

3,277 

DNR 

DNR 

DNR 

Isp  (sec) 

157 

152 

265 

Products  (mole  fraction): 

H20 

0.376 

0.571 

0.227 

N2 

0.119 

0.286 

0.326 

02 

0.287 

0.143 

HCI 

0.197 

CI2 

0.019 

CO 

0.246 

C02 

0.082 

H2 

0.088 

OH 

0.013 

H 

0.012 

4 


Figure  1  shows  the  calculated  performance,  adiabatic  flame  temperature  and 
molecular  weight  of  each  of  the  three  alkali  metal  nitrates  (compared  to  an  AN 
propellant  as  a  reference  condition).  The  calculations  were  performed  assuming  a 
10%  HTPES  binder  and  the  percentage  of  aluminum  was  varied  with  the  percentage  of 
the  nitrate. 

Calculated  values  of  Isp  for  the  corresponding  AP  propellant  were  included  in 
the  Isp  plot.  The  Isp  values  are  very  low  (as  expected)  with  the  Isp  values  in  the  order 
LN>SN>KN.  The  Isp  of  the  AN  propellant  is  significantly  higher  than  that  of  the  three 
alkali  metal  nitrates  due  to  the  much  lower  molecular  weight  of  its  products.  The 
ammonium  ion  provides  a  source  of  hydrogen  atoms  which  results  in  20  to  30% 
hydrogen  in  the  exhaust  products,  and  the  very  low  molecular  weight.  The  large 
fraction  of  hydrogen  in  the  AN  propellant  exhaust  is  part  of  the  attractiveness  of  AN 
propellants.  By  contrast  the  alkali  metals  form  their  corresponding  carbonates  as 
exhaust  products,  which  have  very  high  molecular  weights  and  contribute  to  their  low 
Isp  values. 

The  comparison  between  AN  and  AP  propellants  is  interesting.  The  AP 
propellant  has  a  maximum  Isp  at  22%  aluminum  but  the  curve  is  a  very  broad  curve. 
On  the  other  hand,  the  AN  propellant  has  a  maximum  Isp  at  28%  aluminum  and  the 
curve  has  a  much  sharper  peak.  The  two  have  equivalent  Isp  values  at  26% 
aluminum.  Unfortunately,  it  is  difficult  to  get  good  combustion  efficiency  above  about 
20%  aluminum,  which  reduces  the  potential  practical  usage  of  AN  propellants. 

It  is  significant  that  the  LN  flame  temperature  is  significantly  higher  than  the 
other  flame  temperatures.  It  could  provide  a  source  of  high  temperature  near  the 
surface  to  help  as  an  ignition  aid.  However,  it  is  very  hygroscopic  as  are  most  of  these 
compounds,  which  would  probably  preclude  its  practical  usage. 

It  is  anticipated  that  AP  will  be  required  for  a  practical  AN  containing  propellant 
in  order  to  achieve  the  desired  combustion  efficiency  and  burn  rate  control  (in  the 
same  way  that  AP  is  used  in  double  base  propellants).  With  this  in  mind,  a  series  of 
calculations  were  made  examining  the  performance  of  propellants  containing  a 
mixture  of  AP,  AN  and  an  alkali  metal  nitrate.  The  percent  alkali  metal  nitrate  was 
added  in  a  stoichiometric  ratio  of  AP  to  nitrate,  calculated  to  scavenge  all  of  the 
chloride  ion.  For  example. 


NH4CIO4  +  XNO3  =  xCL  +  . 

Thus,  the  mole  ratios  of  0.59,  0.72,  and  0.86  were  used  for  AP  to  LiNOa,  NaNOs, 
and  KNO3,  respectively.  The  results  for  a  10%  HTPB  propellant  having  10  or  20%  AP 
with  varying  amounts  of  aluminum  are  shown  in  Figure  2.  As  expected  the  AP/LN 
propellant  family  shows  the  best  performance,  with  the  AP/KN  propellant  family  having 
the  worst  performance.  Many  of  the  propellants  have  calculated  Isp  values  in  the  250 
to  260  range  which  are  not  unreasonable.  The  calculated  percentage  of  HCI  in  the 
exhaust  for  all  of  the  propellants  was  below  2%  and  is  shown  in  Figure  3.  The 
dramatic  decrease  in  HCI  concentration  at  27%  aluminum  is  due  to  the  preferential 
formation  of  AlCl  at  the  very  high  aluminum  concentrations. 
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%  Aluminum 

Figure  2.  Comparison  of  Isp  values  for  mixed  oxidizer  systems  containing  ammonium 
nitrate  plus  stoichiometric  amounts  of  alkali  metal  nitrates  needed  to  scavenge 
chloride  ion  from  different  levels  of  AP. 


%  Aluminum 

Figure  3.  Comparison  of  HCI  exhaust  concentrations  for  mixed  oxidizer  systems 
containing  ammonium  nitrate  plus  stoichiometric  amounts  of  alkali  metal  nitrates 
needed  to  scavenge  chloride  ion  from  AP. 
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Although  SN  did  not  appear  to  be  as  effective  as  lithium  nitrate,  it  is  less 
hygroscopic  and  more  compatible  with  the  other  composite  propellant  ingredients  than 
the  LN.  Therefore,  it  appeared  to  be  the  most  desirable  of  the  scavenging  nitrates,  and 
ultimately  became  the  additive  of  choice  for  the  study.  In  addition  to  these  properties 
the  SN  also  made  a  significant,  positive  impact  on  the  agglomeration  tendencies  of  the 
propellants,  • 

Additional  calculations  were  made  to  examine  the  effects  of  SN  on  propellants 
of  interest.  The  calculations  were  made  for  a  12%  HTPB  binder  with  20%  aluminum 
and  10  or  20%  AP.  The  amount  of  SN  was  varied  from  0  to  10%  with  the  remainder  of 
the  formulation  being  AN.  The  calculations  show  that  both  the  isp  and  the  percentage 
HCI  decrease  with  increasing  SN  content.  The  results  are  summarized  in  Figure  4 
where  the  Isp  has  been  plotted  versus  the  percent  HCI  in  the  exhaust.  Although  the 
SN  is  a  very  effective  scavenger,  it  also  reduces  the  Isp  significantly. 


%  HCI  in  Exhaust 

Figure  4.  A  nnmnari.son  of  calculated  Iso  values  versus  the  concentration  of  HCI  in  the 
exhaust  of  orooellants  containing  AN.  AP  and  SN. 

The  calculated  Isp  values  are  on  a  per  unit  mass  basis  and  the  figure  of  merit  for 
a  large  space  booster  is  m'ore  logically  based  on  the  volume  of  the  booster,  which 
corresponds  to  the  product  of  the  Isp  and  the  density.  Therefore,  the  Isp-density  has 
been  examined  to  evaluate  the  performance  of  SN  in  an  AN  propellant  on  this  basis. 
Figure  5  shows  the  calculations  from  Figure  4  replotted  as  Isp-density  versus  percent 
HCI  in  the  exhaust.  Sodium  nitrate  has  a  much  higher  density  (see  Table  2)  than 
either  AN  or  AP  (2.26  versus  1.73  and  1.95,  respectively),  and  the  increase  in 
propellant  density  due  to  the  trade-off  between  SN  and  the  lower  density  AN  makes 
the  SN  look  much  more  attractive  on  the  Isp-density  basis.  However,  it  is  also 
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apparent  that  the  increased  amount  of  AP  also  improves  the  performance  due  to  its 
higher  density  and  inherently  higher  performance. 


%  HC!  in  Exhaust 


Figure  5.  A  comparison  of  calculated  Iso-densitv  values  versus  the  concentration  of 
HCI  in  the  exhaust  of  propellants  containing  AN.  AP  and  SN. 


During  the  course  of  the  contract^^  curing  problems  were  encountered  with  the 
nickel  oxide  phase  stabilized  AN.  Ultimately  a  potassium  nitrate  phase  stabilized  AN 
(ANPN)  was  selected  as  the  standard  AN  to  be  used  throughout  the  contract.  The 
ANPN  contains  approximately  11%  KNO3  (KN),  added  for  phase  stabilization,  but 
which  obviously  will  provide  some  scavenging  capability  for  hydrogen  chloride. 

A  series  of  calculations  were  made  to  evaluate  the  performance  of  what  might 
be  candidate  type  propellants  using  the  ANPN  material.  The  baseline  propellant 
formulation  was  assumed  to  be  10%  HTPB  with  20%  aluminum.  AP  content  was 
varied  with  SN  and  with  the  ANPN.  Figure  6  is  a  plot  similar  to  Figure  5  where  Isp- 
density  has  been  plotted  versus  percent  HCI  in  the  exhaust.  Different  levels  of  AP 
were  assumed  and  the  amount  of  SN  was  traded  off  versus  the  amount  of  ANPN.  For 
all  levels  of  AP  up  to  40%,  it  was  possible  to  add  sufficient  SN  to  scavenge  the  HCI  to 
less  than  2%  of  the  exhaust  gases.  The  40%  AP  calculation  required  30%  SN  as  a 
scavenger,  thus  precluding  any  ANPN  at  all.  It  would  not  be  possible  to  keep  the  HCI 
content  below  2%  for  AP  levels  above  40%. 
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Figure  6.  A  comparison  of  calculated  Isp-densitv  values  versus  the  concentration  of 
HOI  in  the  exhaust  of  Drooellants  containing  ANPN.  AP  and  SN. 


The  same  calculations  were  used  to  identify  the  amount  of  SN  necessary  to 
keep  the  HCI  at  or  below  2%  and  the  results  are  plotted  in  Figure  7  as  Isp-density 
versus  percent  AP.  Increasing  AP  requires  increasing  SN  up  to  the  40%  AP  level 
where  all  of  the  ANPN  is  eliminated.  The  performance  is  impressive  in  that  the 
elimination  of  the  ANPN  dramatically  improves  the  performance.  The  low  inherent 
performance  of  AN  plus  its  low  density  result  in  very  unimpressive  performance 
numbers  for  the  propellant.  The  efficiency  of  the  SN  as  a  scavenger  is  shown  in 
Figure  8,  where  the  data  from  Figure  6  have  been  replotted  as  percent  HCI  versus  SN 
content  for  increasing  levels  of  AP.  The  ratio  of  SN  to  AP  needed  to  scavenge  the  HCI 
increases  as  the  AP  content  increases  due  to  the  decreasing  amount  of  KN  present  as 
the  ANPN  content  is  reduced.  From  a  practical  point  of  view,  AN  does  not  appear  to 
be- a  very  good  candidate  for  space  booster  applications.  The  combination  of  AP  and 
SN  appears  to  be  much  better. 
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Figure  7.  A  comparison  of  calculated  Isp-densitv  values  versus  the  concentration  of 
AP  for  propellants  containing  ANPN.  AP  and  SN. 


0  10  20  30  40 

%  SODIUM  NITRATE 


Figure  8.  A  comparison  of  the  HCI  exhaust  concentration  versus  the  concentration  of 
SN  for  propellants  containing  ANPN.  AP  and  SN. 
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Due  to  the  nitrate  in  AN,  it  had  been  speculated  that  the  AN  flame  structure' 
might  be  similar  to  double  base  propellants  or  HMX,  where  a  dual  stage  flame  exists. 
In  an  attempt  to  determine  what  the  flame  structure  might  be,  a  series  of 
thermochemical  calculations  were  made  that  correspond  to  different  potential  products 
in  the  flame.  The  results  are  summarized  in  Table  3  where  the  heat  of  reaction  and 
adiabatic  reaction  temperature  are  recorded  for  five  possible  sets  of  products.  The  key 
molecule  that  characterizes  the  energy  release  is  nitrogen,  as  it  is  reduced  from  the 
nitrate  ion  through  the  various  nitrogen  oxides  to  N2. 


Tahip  3  AN  Monopropellant  Reactions 


AH  (Kean 

cal/am 

im 

(1) 

NH4NO3  =  HNO3  +  NH3 

44.5 

556 

- 

(2) 

NH4NO3  =  NO2  +  H2O  +  H2  +I/2N2 

37.6 

470 

- 

(3) 

NH4NO3  =  NO  +2H2O  +  1/2  N2 

-6.6 

-83 

550 

(4) 

NH4NO3  =  N2O  +  2H2O 

-8.7 

-109 

625 

(5) 

NH4NO3  =  N2  +  2H2O  +I/2O2 

-28.2 

-352 

1247 

The  initial  decomposition  of  AN  to  HNO3  and  NH3  and  the  reaction  leading  to 
NO2  are  both  very  endothermic  and  unlikely  to  occur  in  a  flame  but  are  more  likely  to 
occur  in  the  condensed  phase.  The  other  three  reactions  were  considered  as 
potential  flame  reactions.  The  most  obvious  would  be  to  assume  that  the  flame  is  in 
complete  equilibrium  with  final  products  of  N2,  H2O,  and  02-  Initially,  this  was  not 
considered  to  be  very  likely  because  of  the  dark  zone  flame  structure  in  double  base 
flames,  where  NO  forms  in  a  pseudo  equilibrium  condition.  The  third  reaction  would 
be  similar  to  double  base  propellants  with  a  dark  zone  flame  where  the  principal 
nitrogen  product  would  be  NO.  This  leads  to  final  products  of  NO  and  H2O.  However, 
this  possibility  has  been  rejected  for  two  reasons.  First,  experimental  decomposition 
data32  Indicate  that  N2O  appears  to  be  the  principal  nitrogen  product  coming  off  the 
propellant  surface.  Apparently  the  melting  and  subsequent  condensed  phase 
reactions  reduce  the  nitrogen  past  the  NO  step  to  N2O.  Secondly,  the  resultant  flame 
temperature  calculated  for  NO  as  a  product  is  less  than  600  K,  which  is  the  norninally 
measured  surface  temperature  of  burning  AN.  The  fourth  reaction  has  been  rejected 
for  the  same  two  reasons.  Therefore,  it  would  appear  that  the  nitrate  is  reduced  at 
least  to  NO  in  the  condensed  phase  and  probable  some  of  the  reaction  proceeds  to 
N2O,  which  is  the  predominant  nitrogen  product  leaving  the  surface.  Thus  it  seerns 
unlikely  that  there  is  an  inner  flame  for  AN  combustion  where  either  NO  or  N2O  is  in 
equilibrium.  Based  on  these  various  reasons  it  has  been  assumed  that  the  flame  is  in 
complete  equilibrium  with  final  products  of  N2.  H2O,  and  O2  and  at  a  monopropellant 
flame  temperature  of  1247  K.  The  need  for  flame  temperature  measurements  is 
readily  apparent.  It  would  be  of  great  benefit  to  have  at  least  one  thermocouple 
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measurement,  or  a  concentration  measurement  to  determine  the  flame  temperature  or 
the  species  present  in  the  AN  monopropellant  flame. 


A  critically  important  part  of  the  model  is  the  primary  diffusion  flame  that  occurs 
between  the  decomposition  products  of  the  binder  and  the  oxidizer.  For  AP  composite 
propeilants  the  primary  diffusion  flame  is  the  dominant  combustion  mechanism  in 
determining  the  burning  rate  of  the  propellant33.  Understanding  the  characteristics  of 
the  diffusion  flame  in  AN  propellants  is  prerequisite  to  the  development  of  a  realistic 
propeilant  combustion  modei..  Flame  temperature  calculations  have  been  made  for 
HTPB/AN  propellants  to  aid  in  identifying  the  flame  structure  and  for  incorporation  into 
the  computer,  program.  As  it  has  been  discussed  in  the  previous  section,  it  had  been 
speculated  that  the  AN  flame  structure  may  be  similar  to  double  base  propellants  due 
to  the  nitrate  in  AN.  in  an  attempt  to  determine  what  the  flame  structure  might  be,  a 
series  of  thermochemical  calculations  were  made  corresponding  to  different  potential 
flames.  The  calculations  were  made  varying  the  percent  AN  in  an  HTPB  binder  over  a 
range  of  pressures  from  100  to  10,000  psi.  Typical  results  are  shown  in  Figure  9. 


Weight  %  AN 


Initially  four  separate  flame  structures  were  postulated.  The  most  obvious  is  to 
assume  that  the  flame  is  in  complete  equilibrium  (the  solid  curve  in  the  figure).  This 
was  not  considered  to  be  very  likely.  The  other  extreme  was  assumed  to  be  similar  to 
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double  base  propellants  with  a  dark  zone  flame  where  the  principal  nitrogen  product 
would  be  NO.  However,  this  possibility  has  been  rejected  for  the  reasons  discussed 
above  pertinent  to  the  AN  monopropellant  flame.  Therefore,  it  was  anticipated  that  if 
there  were  an  inner  flame  for  AN  combustion  it  would  correspond  to  the  case  where 
the  N2O  is  in  equilibrium  (this  is  represented  by  the  curve  marked  'suppress  CO2  and 
N2'  in  the  figure),  and  corresponds  to  a  calculated  monopropellant  flame  temperature 
of  650  K.  As  discussed  above  this  possibility  of  an  inner  monopropellant  flame 
reacting  to  N2O  seems  highly  unlikely.  The  corresponding  flame  between  AN  and 
HTPB  seems  equally  unlikely  due  to  the  very  low  temperatures  that  are  calculated  (i.e. 
the  temperatures  are  consistently  below  1000  K).  The  final  alternative  to  be 
considered  would  be  the  case  where  it  is  assumed  that  the  nitrogen  reacts  completely 
to  N2  and  there  is  no  inner  flame  except  for  the  delayed  reaction  of  CO  to  CO2.  This  is 
shown  in  the  figure  as  the  curve  marked  'suppress  CO2'.  It  was  finally  concluded  that 
the  primary  flame  between  AN  and  the  binder  occurs  with  N2,  CO  and  H2O  as  the  final 
products.  Previous  work  on  AP6.33  and  other  systems^^  leads  to  the  conclusion  that 
the  formation  of  CO2  is  relatively  slow  and  that  CO  is  the  product  of  importance  near 
the  burning  surface. 

The  flame  temperatures  were  calculated  as  a  function  of  pressure  and 
concentration  with  the  CO2  suppressed.  The  results  are  presented  in  Figure  10.  The 
peak  temperature  occurs  at  91  percent  AN  or  at  an  oxidizer  to  fuel  ratio  of  10.1 .  The 
effect  of  pressure  is  to  cause  the  temperature  to  increase  with  pressure  at  low  oxidizer 
concentrations,  but  virtually  no  effect  of  pressure  is  observed  above  82  percent  AN. 
These  are  the  data  that  have  been  programmed  into  the  computer  program  to  provide 
the  primary  flame  temperature  as  a  function  of  pressure  and  oxidizer  to  fuel  ratio. 


Figure  10.  Calculated  primary  flame  temperatures  for  the  CO  equilibrium  flame  for 
AN/HTPB  propellants  for  use  in  the  model. 
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It  is  also  informative  to  look  at  the  composition  of  the  gas  for  the  same 
conditions.  Figure  11  shows  the  variation  of  gas  composition  at  1000  psi.  Below 
about  70  percent  AN  a  significant  amount  of  solid  carbon  is  predicted  to  form,  and  the 
numbers  in  the  figure  do  not  reflect  the  free  carbon  in  the  products.  At  50  percent  AN 
29  percent  of  the  exhaust  products  are  predicted  to  be  solid  carbon.  Above  70  percent 
no  free  carbon  is  predicted.  The  very  high  percentage  of  H2  in  the  products 
contributes  to  a  very  low  molecular  weight  of  the  gases  which  makes  the  Isp 
reasonable  for  AN  propellants. 


Weight®/*  AN 

Figure  1 1 .  Exhaust  gas  composition  for  AN/HTPB  propellants  varvina  percent  AN. 


Other  Primary  Flame  Temperatures 

Because  the  propellants  to  be  considered  were  to  have  multiple  oxidizers, 
flame  temperature  calculations  were  also  made  for  AP  and  $N  propellants.  The  same 
approach  was  taken  as  described  above,  i.e.  the  formation  of  CO2  was  suppressed 
and  the  calculations  were  made  as  a  function  of  oxidizer  fraction  and  pressure. 
Typical  results  are  shown  graphically  in  Figure  12.  The  peak  of  the  curves  for  AP  and 
SN  are  nearly  the  same  weight  fraction  of  oxidizer  while  the  AN  peak  occurs  at  a 
higher  oxidizer  fraction.  The  AP  system  has  a  peak  temperature  of  2750  K  and  SN 
has  a  peak  temperature  slightly  greater  than  2300  K  while  the  peak  temperature  for 
AN  is  barely  2000  K.  It  is  apparent  from  this  very  low  temperature  that  an  AN 
propellant  would  have  a  much  lower  burn  rate  and  its  propensity  to  ignite  and  burn 
aluminum  would  be  very  low. 
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Figure  12.  Cnmnarison  of  calculated  flame  temperatures  for  primary  flames  in  AP. 

AN.  and  SN/HTPB  prooellants. 

Performance  Characteristics  of  AN/Silicone  Rubber  Propellants 

In  a  parallel  contract  at  Hercules  (the  Aluminum-Water  Rocket  Program)35, 
propellant  evaluation  work  was  performed  developing  a  highly  metallized  AN  based 
propellant.  Under  that  contract  it  was  decided  to  use  a  silicone  rubber  binder  (Sylgard 
182)  due  to  the  improved  combustion  characteristics  that  were  obtained  using  that 
binder  within  a  highly  metallized  propellant. 

For  comparative  purposes  calculations  were  performed  for  a  series  of 
propellants  identical  to  those  in  Figure  12,  but  with  the  silicone  rubber  binder  in  place 
of  HTPB.  The  silicone  rubber  has  a  nominal  formulation  of  SiCHaO  with  a  heat  of 
formation*  of  -24.9  Kcal/mole.  The  results  are  displayed  in  Figure  13.  The  flame 
temperatures  of  the  silicone  based  propellants  are  dramatically  different  for  each  of  the 
three  oxidizers,  giving  higher  flame  temperatures  and  much  broader  peaks  in  the 
curves  than  the  corresponding  HTPB  propellants.  The  extra  oxygen  in  the  silicone 
rubber  broadens  the  curve  and  causes  the  peak  to  shift  to  a  much  lower  oxidizer 
fraction.  The  increase  in  temperature  is  apparently  due  to  the  differeht  exothermicity 
of  the  products  from  the  two  binders.  The  Sylgard  binder  forms  Si02  in  addition  to  the 
CO  and  CO2  formed  by  the  HTPB  binder.  The  higher  heat  of  formation  of  Si02  relative 
to  CO  and  CO2  apparently  accounts  for  the  increased  temperature.  The  net  result 
being  the  observed  shift  in  both  the  flame  temperature  and  the  stoichiometric  ratio  of 
oxidizer  to  binder. 


16 


The  propellants  used  in  the  Aluminum-Water  Rocket  Program's  were  ' 
formulated  with  approximately  18  percent  binder.  From  the  figure  it  is  apparent  that  at 
18  percent  binder  loading  the  silicone  propellant  flame  temperature  is  significantly 
higher  than  a  corresponding  HTPB  binder  propellant  and  one  would  expect  the 
silicone  propellant  to  burn  better,  both  faster  and  more  efficiently,  which  was  observed 
experimentally. 


Figure  13.  Comparison  of  calculated  flame  temperatures  for  primary  flames  in  AP. 

AN,  and  SN  propellants  with  HTPB  and  Sviaard  binders. 

Performance  calculations  were  also  made  comparing  the  different  binder 
systems,  and  are  presented  in  Figures  14  and  15.  Figure  14  displays  the  calculated 
Isp  values  for  AP  and  AN  systems  with  both  the  Sylgard  binder  as  well  as  an  HTPB 
binder  for  reference.  Although  the  Sylgard  binder  produces  Si02  which  gives  the 
Sylgard  propellants  an  increased  molecular  weight,  the  flame  temperature  is 
sufficiently  higher  than  the  HTPB  propellants  that  the  AN/Sylgard  propellant  has  a 
much  superior  Isp  than  the  AN/HTPB  propellant.  The  AP  based  calculations  reach 
comparable  Isp  values  but  at  very  different  binder  loadings.  Figure  15  compares  the 
four  different  cases  plotted  versus  the  aluminum  content.  The  additional  oxygen  in  the 
silicone  binder  apparently  reduces  the  overall  performance  slightly,  so  that  the 
calculated  Isp  values  for  the  AP  propellant  are  about  10  sec  below  the  HTPB  values. 
However,  the  AN  silicone  system  is  very  comparable  to  the  HTPB  system  (they  are 
equal  at  ~17%  binder),  and  considering  combustion  efficiency,  the  silicone  binder 
would  appear  to  be  a  better  choice  than  the  HTPB  binder.  Above  20%  aluminum  the 
HTPB  binder  gives  significantly  improved  performance,  but  obtaining  good 
combustion  efficiency  in  that  range  of  propellant  formulations  could  be  very  difficult. 
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Figure  14.  Comparison  of  calculated  Isp  for  AP  and  AN  Dropellants  with  HTPB  ancj 
Svloard  binders  varying  binder  level. 


Figure  15.  Comnari.son  of  calculated  Iso  for  AP  and  AN  oropellantswith  MTPB  aad 
Svigard  binders  varying  aluminum  level.. 
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COMBUSTION  MECHANISMS  OF  SOLID  PROPELLANTS 


The  Condensed  Phase 


The  two  oxidizers  used  most  often  in  modern  solid  propellants  are  ammonium 
perchlorate,  AP,  and  cyclotetramethyline  tetranitramine,  HMX.  In  this  study  two 
additional  oxidizers  have  been  considered.  These  are  ammonium  nitrate,  AN.  and 
sodium  nitrate,  SN.  The  motivation  for  studying  them  is  to  develop  an  alternate 
oxidizer  for  AP  that  would  be  cheaper  and  would  be  "clean",  i.e.  would  not  produce 
HCI  in  the  exhaust.  The  two  most  common  binder  systems  consist  of  either  an  inert 
binder  system  based  on  a  polybutadiene  polymer  (almost  always  HTPB),  with  various 
plasticizers,  additives  and  cross  linking  agents  or  an  active  binder  system  based  on 
energetic  components  such  as  nitroglycerin,  NG,  and  nitrocellulose,  NC.  As 
referenced  in  the  previous  chapter,  a  silicone  based  binder  has  also  been  used  rather 
successfully  in  a  program  using  AN  as  the  propellant  oxidizer.  However,  the  current 
study  is  limited  to  the  more  common  HTPB  binder  propellants. 

As  a  propellant  burns,  both  oxidizer  and  binder  are  preheated  by  conduction  as 
they  approach  the  surface.  The  initial  step  in  the  overall  process  is  probably  melting 
which  results  in  a  thin  liquid  layer  on  the  oxidizer  surface  (probably  on  the  binder 
surface  also).  AP  and  HMX  melting  points  are  approximately  865  K  and  553  K, 
respectively  (see  Table  1).  AN  melts  at  443  K  and  SN  melts  at  580  K,  both  of  which 
are  much  lower  than  the  melting  point  of  AP.  The  very  low  melting  point  of  AN  will 
result  in  a  relatively  long  residence  time  for  the  molten  material,  a  thicker  melt  layer, 
and  could  allow  for  some  lateral  mixing  at  the  propellant  surface.  In  the  liquid  layer 
and  immediately  below  the  propellant  surface  the  initial  decomposition  step  is 
apparently  the  endothermic  dissociation  of  the  oxidizer.  Ammonia  and  perchloric  acid 
are  the  initial  products  for  AP,  while  AN  yields  ammonia  and  nitric  acid.  While  still  near 
the  surface  in  the  liquid  layer,  these  decomposition  products  apparently  undergo 
further  exothermic  condensed  phase  reactions  before  passing  through  the  surface  into 
the  gas  phase  (see  Table  3  in  the  previous  chapter).  These  additional  reactions  are 
exothermic,  to  the  extent  that  the  overall  condensed  phase  decomposition  is 
exothermic.  The  possibility  of  similar  exothermic  reactions  occurring  at  the  binder 
surface  must  also  be  considered.  These  could  occur  by  reaction  with  oxygenated 
species  coming  from  within  the  binder  itself  or  by  diffusion  of  reactive  species  to  the 
binder  from  the  adjacent  oxidizer. 

The  products  observed  from  the  low  pressure,  high  heating  rate  decomposition 
of  ammonium  perchlorate  include  chlorine,  nitrogen  oxides,  water,  hydrogen  chloride 
as  well  as  ammonia  and  perchloric  acid.  The  corresponding  tests  for  AN  yield  a 
mixture  of  NO,  N2O,  N2  and  water.  As  was  discussed  in  the  previous  section,  the 
gases  coming  from  the  oxidizer  surface  are  a  mixture  of  partially  decomposed 
products.  A  significant  amount  of  the  available  energy  is  apparently  released  in  the 
condensed  phase.  During  combustion  it  is  anticipated  that  the  gases  coming  from  the 
oxidizer  surface  contain  little  or  none  of  the  initial  decomposition  products  but  are 
almost  exclusively,  the  partially  decomposed  products. 

A  question  to  be  addressed  is,  how  important  or  relevant  are  these  initial 
reactions  to  the  ultimate  burning  rate  of  the  propellant.  Boggs,  et.al.36  jn  a  recent 
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paper  on  AP(ammonium  perchiorate),  concludes  that  low  temperature  decomposition 
studies  can  and  do  lead  to  conflicting  predictions  relative  to  actual  deflagration.  The 
authors  quote  many  instances  showing  that  low  temperature  decomposition  tests  do 
not  correlate  with  actual  combustion  conditions.  The  initial  decomposition  reactions 
apparently  only  provide  a  boundary  condition  of  gaseous  products  and  temperature 
for  the  energy  feedback  from  the  flame.  As  the  energy  feedback  varies,  the  reactions 
at  the  surface  apparently  adjust  to  provide  a  different,  compatible  composition  and 
temperature  for  the  given  energy  feedback. 

Surface  Temperature 

The  surface  temperature  is  related  to  the  burning  rate  in  most  models  as  a  one- 
step,  zero-order  Arrhenius  pyrolysis  law 

r  =  (1 ) 

Based  on  Equation  (1)  the  surface  temperature  is  directly  dependent  on  the  burning 
rate  but  only  indirectly  dependent  on  the  pressure  (due  to  the  fact  that  the  burning  rate 
is  dependent  on  pressure). 

Only  two  sources  of  surface  temperature  data  for  AN  were  found  in  the 
literature.  Anderson,  et  al^^  reported  hot  plate  data  for  pellets  of  AN  in  1958.  in  1964 
Whittaker  and  Barham38  measured  temperature  profiles  using  thermocouples  in 
burning  pressed  pellets  of  AN  containing  up  to  2.5%  CraOa.  Their  data  were  taken 
over^a  range  of  pressures  from  1000  to  4800  psi.  It  is  curious  that  they  concluded  that 
the  Anderson  data  were  incorrect.  Figure  16  shows  the  combined  data  plotted  as  rate 
versus  1/K.  The  Whittaker/Barham  combustion  data  are  much  more  scattered  that  the 
hot  plate  data.  Within  the  data  scatter,  the  data  are  actually  in  very  reasonable 
agreement,  with  an  apparent  activation  energy  of  10.4  Kcal/mole.  It  is  significant  that 
the  burn  rate  catalyst  did  not  appear  to  have  a  significant  influence  on  the  surface 
decomposition  characteristics. 

Two  sources  of  data  were  found  for  the  surface  temperature  of  SN.  Frolov,  et 
a|30  report  three  datum  points  for  composite  mixtures  of  SN,  aluminum  and  an 
additive,  but  do  not  say  what  method  was  used  to  measure  the  surface  temperature. 
Girdhar  and  Arora39  report  surface  temperatures  for  SN  composite  propellants,  but 
they  infer  the  surface  temperatures  from  temperature  sensitivity  data  and  an  empirical 
relationship  between  the  rate,  the  surface  temperature  and  the  temperature  sensitivity. 
Their  data  indicate  an  activation  energy  less  than  zero  which  is  considered  to  be 
inconsistent  and  unreliable.  Therefore,  their  data  were  not  used.  Frolovs  data  are 
plotted  in  Figure  17.  The  data  indicate  a  very  high  surface  temperature  for  SN,  in  the 
neighborhood  of  1000  K,  and  a  relatively  low  activation  energy,  4.5  Kcal/mole. 
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Figure  16.  Surface  temperature  data  for  pressed  pellets  of  AN.  The  hot  plate  data  are 
from  Anderson  and  the  combustion  data  from  Whittaker. 
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Figure  17.  Surface  temperature  data  bv  Frolov,  et  al  for  a  composite  mixture 
containing  SN  and  aluminum. 
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AN  and  SN  are  to  be  used  in  composite  propellants  with  AP  as  a  co-oxidizer. 
For  comparative  purposes  the  data  from  Figures  16  and  17  have  been  combined  with 
data  for  AP40.41.42  and  are  plotted  in  Figure  18.  In  addition  the  calculated  values  that 
are  used  for  HTPB  binder  in  the  model  are  also  included  in  the  figure.  The  wide 
diversity  in  surface  temperatures  poses  some  interesting  possibilities.  For  example,  if 
a  propellant  were  burning  at  a  rate  of  one  cm/sec,  for  all  of  the  components  to  be  in 
equilibrium,  the  surface  temperatures  would  have  to  be  700  K  for  AN,  1000  K  for  AP , 
1400  K  for  HTPB,  and  much  higher  for  SN.  Obviously  this  is  not  realistic.  However, 
there  must  be  a  wide  diversity  of  surface  temperatures  and  individual  rates  occurring 
simultaneously  on  a  propellant  surface.  To  model  such  a  situation,  one  must  allow  for 
separate  surface  temperatures  of  each  ingredient  and  individual  rates  for  each 
ingredient.  An  ingredient  such  as  SN  should  burn  very  rapidly  once  it  ignites,  possibly 
leaving  an  empty  pocket  of  binder  or  possibly  disrupting  the  adjacent  surface  regions 
with  an  explosive  type  of  flow  due  to  the  very  high  relative  rate.  AN  on  the  other 
extreme  could  possible  puddle  and  flow  along  a  surface  due  to  its  high  residence  time 
and  low  combustion  rate.  The  physical  picture  for  such  an  environment  is  obviously 
very  complicated. 
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Figure  18.  A  r.nmr^arison  of  surface  temperature  data  for  AP.  AN.  SN  and  HTPB.  the 
principal  ingredients  of  interest  in  this  study. 
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Condensed  Phase  Heat  RelPasP 


A  critical  step  in  the  combustion  process  is  the  energy  release  in  the  condensed 
phase.  Attempts  have  been  made  during  this  contract  to  calculate  this  condensed 
phase  heat  release,  but  there  isn’t  sufficient  information  about  the  condensed  phase 
reactions  to  allow  a  general  calculation.  To  calculate  a  heat  of  reaction,  the  enthalpies 
of  the  reactants,  and  products  must  be  known  at  the  initial  temperature  and  the 
reaction  temperature.  Based  on  the  calculations  reported  in  Table  3  of  the  previous 
chapter,  the  heat  release  in  the  condensed  phase  should  be  on  the  order  of  80  to  1 10 
cal/gm  exothermic.  In  a  parallel  contract  studying  double  base  propellants^s,  it  was 
determined  from  an  examination  of  experimental  data  that  the  heat  of  reaction’will  be 
proportional  to  the  reaction  temperature  and  decreases  with  increasing  initial 
temperature.  The  data  were  correlated  by  plotting  the  surface  heat  release  versus  the 
difference  of  Tg-To,  which  provided  an  excellent  correlation  of  the  data.  In  the  model 
the  condensed  phase  heat  release  was  assumed  to  be 


Ql  =  Qc(Ts-To)  (2) 

where  Qq  is  a  proportionality  constant.  The  reader  is  referred  to  Reference  43  for  a 
complete  presentation  of  the  data  and  logic  supporting  this  equation. 

Oxidizer  Burning  Characteri.qtir:.g; 

The  relative  burning  rates  of  HMX,  AP.  and  AN  as  monopropellants  and  in 
composite  propellant  mixtures  are  shown  in  Figure  19.  As  monopropellants  the  rate  of 
HMX  is  slightly  higher  than  that  of  AP  and  they  both  have  similar  burn  rate  exponents 
for  typical  solid  propellant  rocket  pressures.  The  burning  rate  of  AN  is  approximately  a 
factor  of  five  lower  than  the  other  two  monopropellants.  Although  the  difference 
between  the  linear  burning  rates  of  AP  and  HMX  does  not  appear  to  be  significant,  but 
when  the  two  oxidizers  are  mixed  into  a  composite  propellant,  the  resultant  burn  rates 
vary  by  as  much  as  an  order  of  magnitude.  AP  composite  propellants  typically  have  a 
rate  higher  than  the  monopropellant,  while  HMX  propellants  typically  have  a  rate  lower 
that  the  monopropellant.  When  AN  is  added  to  a  binder,  -the  propellant  rate  is 
comparable  to  that  of  the  monopropellant.  Thus,'  it  would  appear  that  AN  composite 
propellants  burn  by  a  mechanism  similar  to  AP  propellants.  Also,  because  of  the  very 
dramatic  effects  observed  with  each  different  oxidizer,  it  must  be  assumed  that  the 
primary  diffusion  flame  between  the  oxidizer  and  binder  decomposition  products  is  the 
dominant  mechanism  controlling  the  burning  rate. 
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Composite  Propellant  Flame  Structure 


As  has  been  discussed  above,  when  an  oxidizer  decomposes  during 
combustion,  a  mixture  of  partially  oxidized  products  are  formed  in  the  thin  liquid  layer 
on  the  propellant  surface.  In  a  propellant  this  mixture  of  partially  decomposed 
products  coming  from  the  oxidizer  surface  is  available  to  either  react  and  form  the 
monopropellant  flame,  or  react  with  the  fuel  decomposition  products  and  form  the 
primary  diffusion  flame.  In  an  AP  propellant  the  reactive  gases  leaving  the  surface  are 
apparently  a  mixture  of  chlorine  oxides  probably  CIO  and  CIO2.  These  species  are 
very  reactive  and  react  readily  with  fuel  species  from  the  binder  to  form  a  very  powerful 
prima|7  diffusion  flame.  When  HMX  decomposes  formaldehyde,  hydrogen  cyanide, 
and  nitrogen  oxides  are  formed  which  are  actually  fuel  rich  (relative  to  CO2).  None  of 
these  products  is  a  very  strong  oxidizer.  They  can  react  to  form  a  monopropeliant 
flame,  or  when  combined  with  fuel  decomposition  products  they  form  a  low 
temperature  diffusion  flame  (due  to  their  fuel  rich  nature).  When  AN  decomposes 
nitrogen  oxides  (somewhat  similar  to  HMX  or  double  base  propellants)  are  formed 
which  are  not  very  strong  oxidizers  compared  to  the  chlorine  oxides.  They  can  react.to 
form  a  monopropellant  flame,  or  when  combined  with  fuel  decomposition  products 
they  also  form  a  low  temperature  diffusion  flame  due  to  their  low  inherent  energy.  SN 
is  somewhat  unique  in  that  it  does  not  form  a  monopropellant  flame.  Only  the  diffusion 
flame  is  formed. 

Figure  20  shows  the  postulated  flame  structure  for  AP,  AN,  and  SN  within  a 
composite  propellant  with  the  various  flame  and  surface  temperatures  indicated. 

The  relative  surface  temperatures  for  each  oxidizer  are  significantly  different 
and  have  been  discussed  previously.  SN  has  the  highest  apparent  surface 
temperature,  with  AP  in  between  and  AN  has  a  very  low  value.  Each  apparently  burns 
with  a  thin  liquid  layer  on  the  surface.  The  characteristics  of  the  liquid  layer  have  not 
been  determined  experimentally,  but  have  been  inferred  from  calculations,  quench 
samples  and  microcinematography.  The  liquid  layer  is  inferred  to  be  on  the  order  of  a 
few  microns  thick  for  AP,  but  should  be  much  thicker  for  AN  perhaps  as  thick  as  100 
microns  due  to  the  lower  melting  point  and  lower  rate.  Because  the-  rate  increases 
with  increasing  pressure,  the  liquid  layer  thickness  will  decrease  with  increasing 
pressure.  Reactions  in  the  liquid  layer  are  probably  relatively  unimportant  at 
pressures  of  interest  and  only  become  significant  at  low  pressures  where  the  flame 
heights  become  large. 
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Figure  20.  Flame  structure  for  AP.  AN.  and  SN  composite  propellants 
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At  low  pressures  (ambient  to  -100  psi)  composite  propellants  probably  burn  as 
a  premixed  flame  with  the  oxidizer  and  binder  decomposition  products  mixing 
completely  before  reaction  occurs.  This  would  especially  be  true  for  small  oxidizer 
particles.  With  increasing  pressure,  it  becomes  more  difficult  for  the  fuel  products  to 
diffuse  into  the  stream  above  the  oxidizer  (and  vice  versa)  aad  the  two  reaction  paths 
(i.e.  the  monopropellant  flame  versus  the  primary  diffusion  flame)  become  competitive. 
At  higher  pressures  where  reactions  occur  more  rapidly,  the  oxidizer  products  react 
preferentially  to  yield  the  normal  monopropellant  flame  before  the  binder  products  can 
diffuse  into  the  oxidizer  stream  and  react.  However,  near  the  binder-oxidizer  interface, 
a  fraction  of  the  oxidizing  products  will  always  react  with  the  binder  products.  To 
properly  describe  the  trade-off  between  these  two  flames  and  the  relative  heat  transfer 
from  them,  is  a  major  challenge  in  modeling  solid  propellants. 

In  addition  to  the  monopropellant  flame  and  the  primary  diffusion  flame  there  is 
also  a  final  diffusion  flame  for  AP  and  AN  propellants.  As  indicated  in  Table  2  in  the 
preceding  chapter,  the  AP  and  AN  monopropellant  products  are  approximately  30% 
and  14%  O2,  respectively.  When  either  the  AP  or  AN  monopropellant  flame  occurs,  it 
is  followed  by  a  diffusion  flame  (a  2-stage  flame)  where  the  oxygen-rich  products  of  the 
flame  can  react  with  the  binder  pyrolysis  products.  This  reaction  should  depend  on 
mixing  only,  as  the  oxygen  is  heated  to  the  temperature  of  the  monopropellant  flame, 
1400  K  or  1250  K,  and  at  this  temperature  oxygen  would  be  expected  to  be  very 
reactive.  All  the  characteristics  of  a  typical  diffusion  flame  can  be  expected  for  this  final 
flame.  In  Figure  20  a  higher  temperature  is  indicated  for  the  final  diffusion  flame  than 
for  the  primary  diffusion  flame.  This  is  based  on  the  assumption  that  the  preheated 
oxygen  will  react  completely  to  CO2,  while  the  primary  diffusion  flame  is  only  expected 
to  react  to  CO. 

Allocation  of  Binder  to  Oxidizer 

A  critical  assumption  in  developing  the  SST  model  relates  to  how  the  binder  is 
proportioned  to  the  various  oxidizer  fractions.  Two  approaches  have  been  considered 
in  the  past.  In  the  SST  approach,  it  is  assumed  that  that  the  primary  flame  temperature 
is  the  adiabatic  flame  temperature  for  the  stoichiometric  mixture  ratio  of  binder  and 
oxidizer.  However,  some  relationship  between  the  decomposing  binder  and  the 
amount  of  these  decomposition  products  that  react  in  the  primary  flame  must  be 
obtained,  which  amounts  to  estimating  how  much  binder  reacts  in  the  primary  flame 
and  how  much  does  not  react.  To  accomplish  this,  a  term  accounting  for  the  fraction  of 
binder  that  reacts  must  be  introduced  into  the  equations.  In  the  approach  used  in  both 
the  original  BDP  model  and  the  PEM  type  models,  all  the  binder  associated  with  a 
given  crystal  is  considered  to  react  with  the  available  oxidizer  in  the  diffusion  flames. 
This  assumption  leads  to  a  variable  flame  temperature  per  oxidizer  size  fraction. 

The  theoretical  concept  of  a  diffusion  flame  as  well  as  experimental  data  dictate 
that  the  flame  burns  at  "stoichiometric”  conditions.  In  other  words,  when  sufficient  fuel 
and  oxidizer  have  diffused  together,  they  will  react,  independent  of  the  overall  0/F 
ratio.  However,  if  there  is  excess  fuel  or  oxidizer  present,  the  excess  material  will 
ultimately  reduce  the  temperature  in  the  rocket.  This  final  mixing  and  equilibration 
probably  occurs  far  enough  from  the  surface  that  it  does  not  effect  the  burn  rate 
significantly. 
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It  must  also  be  noted  that  the  stoichiometric  condition  within  the  SST  model  was 
not  assumed  to  be  that  corresponding  to  final  production  of  CO2  Because  the  CO  to 
CO2  reaction  is  a  relatively  slow  reaction,  an  intermediate  "stoichiometric"  condition 
was  assumed  for  CO  products.  Further  reaction  .and  formation  of  CO2  undoubtedly 
occurs,  and  the  final  temperature  of  the  gases  far  from  the  surface  will  be  the  flame  -  - 
temperature  corresponding  to  the  overall  propellant  formulation.'  However,  it  is  most 
likely  that  this  equilibration  and  'after  burning’  also  occurs  too  far  from  the  surface  to 
have  an  effect  on  the  burning  rate. 

The  assumption  of  a  variable  flame  temperature  is  the  common  assumption 
used  in  most  models,  but  it  is  not  necessarily  on  a  sound  mechanistic  basis. 
Furthermore,  it  leads  to  peculiar  results  due  to  the  fact  that  the  AP/HTPB  flame 
temperature  has  a  maximum  with  respect  to  AP  concentration.  The  PEM  models 
assume  an  unrealistic  distribution  of  binder  per  oxidizer  distribution  (and 
corresponding  flame  temperature)  in  order  to  match  experimental  data.  The  crux  of 
the  problem  is  really  in  how  the  binder  is  assumed  to  be  distributed  in  the  propellant. 

CohenS  has  proposed  and  used  a  compromise  position  which  appears  to  be 
quite  realistic  on  a  mechanistic  basis.  He  assumes  that  the  fuel  left  over  from  the  fine 
particle  size  is  available  to  be  burned  in  the  diffusion  flame  of  the  next  larger  size 
fraction.  This  makes  sense.  Typical  binder  thicknesses  associated  with  the 
intermediate  and  coarse  oxidizer  fractions  are  on  the  order  of  0.2|i,  for  20|i  particles 
and  2p.  for  200  |i  particles  (i.e.  1%  of  the  particle  diameter)44.  This  is  smaller  than  the 
roughness  and  non-sphericity  of  typical  crystals.  Thus,  it  can  be  anticipated  that  there 
will  be  a  significant  amount  of  binder  buried  within  the  surface  irregularities  of  a  crystal 
in  excess  of  what  would  be  calculated  for  a  purely  spherical  particle.  Cohen's 
approach  reflects  this  reality.  In  his  most  recent  model®  he  assumes  that  the  finest 
oxidizer  fraction  burns  at  the  stoichiometric  0/F- condition.  The  amount  of  binder 
available  to  the  fine  fraction  will  be  calculated  to  be  greater  than  the  stoichiometric  0/F 
ratio  (i.e.  the  binder  is  assumed  to  be  distributed  according  to  the  surface  area  of  the 
particles).  Therefore,  there  is  binder  left  over  from  the  fine  fraction  which  Cohen 
assumes  is  available  for  combustion  in  the  next  larger  size  fraction.  That  size  fraction 
is  also  assumed  to  burn  at  stoichiometric  conditions  if  there  is  sufficient  fuel,  and  any 
excess  fuel  is  then  assumed  to  be  available  to  burn  with  the  coarsest  fraction.  This 
approach  appears  to  be  very  realistic  and  has  been  incorporated  into  the  SST  model 
for  application  with  both  AP  and  AN  propellants.  Furthermore,  a  preference  order  of 
oxidizer  must  be  established.  Because  AP  appears  to  be  the  most  reactive  of  the 
oxidizers  considered,  it  has  been  assumed  that  the  AP  will  react  whatever  binder  is 
available  to  achieve  stoichiometric  conditions  and  that  the  AN  is  left  to  use  whatever 
binder  is  left  over. 
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AN  MONOPROPELLANT  MODELING 


The  model  discussed  in  this  section  is  essentially  the  BDP  monopropellant 
mode|27.45,46  that  has  previously  been  applied  to  AP,  HMX  and  double  base  binders 
burning  as  monopropellants.  The  overall  model  is  being  developed  with  the  intent  of 
applying  it  to  composite  propellant  systems.  The  monopropellant  model  is  a  limiting 
condition  of  the  model,  with  no  binder. 

Basic  Model  Equations 

The  BDP  monopropellant  mode|27.45,46  reduces  to  essentially  three  equations: 
one  describing  the  surface  kinetics;  one  resulting  from  an  energy  balance  at  the 
burning  surface;  and  one  related  to  the  gas  phase  kinetics.  The  surface  regression  is 
described  as  a  one-step,  zero-order  Arrhenius  pyrolysis  law  (see  Equation  (1)). 

The  surface  temperature  is  assumed  to  be  uniform  over  the  burning  surface  and 
is  determined  from  an  energy  balance.  The  energy  balance  can  be  expressed  as 

Cp  (Ts  -  To)  -  Ql  =  Qf  6'^ 

where  Tq  is  the  initial  temperature,  Ql  is  the  condensed  phase  heat  release  and  the 
final  term  represents  the  energy  transported  from  the  flame  back  to  the  surface.  For  a 
premixed  laminar  flame,  the  nondimensional  flame-standoff  distance  is  represented 
by  the  expression 

^  =  Cpm2/^^  1^  p5j  (4) 


where  the  5  is  the  reaction  order  and  the  gas-phase  rate  constant  k  is  expressed  by 
an  Arrhenius  function  ofthe  monopropellant  flame. 

k=  (5) 


The  actual  flame-standoff  distance  is  then 


X-  =  m/(kp8)  ,  (6) 

In  the  previous  applications  of  the  Beckstead  mode|2  for  monopropellants,  the 
value  of  Ql  was  taken  to  be  constant  with  varying  pressure,  which  was  considered  to 
be  a  deficiency  in  the  model.  As  a  result  of  the  concurrent  work  being  done  on  the 
Temperature  Sensitivity  Program's,  available  data  for  double  base  propellants  and  for 
HMX  monopropellant  indicate  that  the  condensed  phase  heat  release  varied  with 
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surfacG  tempGratuTG  and  invGrsGly  with  initial  tamparaturG.  Tha  Gxprassion  that  has 
baan  programmad  into  tha  modal  for  Ql  is  shown  in  Equation  (2). 

This  introducGS  the  condensad  phasa  haat  ralaasa  into  the  calculation  as  a 
variable  which  increases  with  increasing  pressure  which  is  in  general  agreement  with 
what  limited  data  are  available.  Pressure  does  not  appear  explicitly  in  Equation  (2), 
but  the  surface  temperature  varies  with  burn  rate  which  is  a  function  of  pressure. 
Although  there  are  no  specific  data  for  AN,  Ql  will  be  assumed  to  follow  the  form  of 
Equation  (2). 

The  predictions  of  the  most  recent  application  of  the  Beckstead  monopropellant 
mode|45  were  very  reasonable.  A  comparison  was  made  for  both  AP  and  HMX,  with 
data  at  several  pressures  and  three  initial  temperatures.  There  was  some  deviation  of 
the  model  from  experimental  data  at  high  pressure,  and  there  was  criticism  of  the 
activation  energy  values  that  were  used.  However,  all  of  the  calculated  trends  were 
consistent  with  observed  experimental  data. 


Pressure  Exponent  and  Temperature  Sensitivity 

One  of  the  advantages  of  the  BDP  approach  is  its  simplicity.  The  equations  are 
sufficiently  tractable  that  they  can  be  differentiated  to  obtain  a  closed-form  solution  for 
both  pressure  and  temperature  dependence.  The  pressure  exponent  can  be 
calculated  to  be 


n  =  [5-»-( 


CpTf  Ef  \dlnTf-i  /r2 


(7) 


For  most  situations  the  pressure  exponent  will  be  ~  half  of  the  reaction  order,  6. 
However,  that  assumes  an  equilibrium  flame  temperature,  which  is  essentially 
independent  of  pressure.  The  monopropellant  flame  is  most  likely  in  equilibrium  and 
therefore,  independent  of  pressure.  For  double  base  propellants  this  is  not  true. 

The  temperature  sensitivity  can  be  calculated  as 


ap=[— + 

Q,J  e-5  RT,  dTo  P-  - 


(8) 


Qf^  e-< 

The  denominator  of  the  expressions  for  n  and  Gp  is  identical  and  therefore,  predicted 
dependencies  will  be  similar.  From  thermochemical  calculations,  the  rate  of  change 
of  flame  temperature  with  initial  temperature  (dTf/dTo)  is  approximately  0.6  for  both  AP 
and  AN.  For  a  nominal  flame  temperature  of  1500  K  and  an  activation  energy  of  15 
Kcal/mole,  the  second  term  in  the  numerator  gives  a  value  of  0.002  1/K  which,  divided 
by  the  value  of  2  in  the  denominator  gives  a  contribution  to  the  Cp  which  is  the  correct 
order  of  magnitude  for  AP  monopropellant  (op  values  are  typically  0.001  to  0.002  1/K). 
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Comparison  of  Monopropellant  Model  with  AP  Data 

Because  of  the  apparent  similarities  between  AP  and  AN,  and  because  of  the 
extensive  previous  work  done  with  AP,  it  was  decided  to  use  AP  as  a  basis  of 
comparison  for  the  monopropeliant  model.  The  key  parameters  that  are  required  to 
make  the  calculation  are  tabulated  in  Table  4  with  a  comparison  of  the  values  used  for 
both  AP  and  AN.  Some  of  the  parameter  values  are  not  available  from  experimental 
data,  and  the  rationale  for  selecting  the  specified  value  will  be  discussed  in  the 
discussion  of  the  figures  that  follow. 


Table  4  Model  Thermochemical  and  Physical  Parameters 


AN 

Density  (gm/cc) 

1.95 

1.73 

Melting  Point  (K) 

865 

443 

Heat  Capacity,  (cal/gm) 

0.30 

0.30 

Gas  Thermal  Conductivity  (cal/cm  sec  K) 

0.0003 

0.0003 

Condensed  Phase  Heat  Release  Coef. 

0.20 

0.26 

Reference  Flame  Temperature  (K) 

1405 

1247 

Reference  Surface  Temperature  (K) 

950 

530 

Reference  Burn  Rate  @  68  atm  (cm/sec) 

0.81 

0.19 

Reaction  Order 

1.7 

2.0 

Flame  Activation  Energy  (Kcal/mole) 

15 

25 

Surface  Activation  Energy  (Kcal/mole) 

30 

14 

The  burning  rate  and  the  burning  rate  temperature  sensitivity,  Op,  have  been 
calculated  as  a  function  of  pressure  and  compared  to  Boggs'  data  for  pressed  pellets 
of  AP47,48,  7he  calculated  results  are  compared  with  the  ambient  temperature  data  in 
Figure  21 ,  and  the  calculated  values  of  ap  are  shown  in  Figure  22.  Equation  (8)  of  the 

model  indicates  that  ap  is  very  dependent  on  the  gas  phase  activation  energy. 
Activation  energy  values  of  10, 15  and  20  Kcal/mole  were  used  to  calculate  the  results 
presented  in  Figure  22.  The  precise  value  of  the  activation  energy  does  not  have  a 
significant  effect  of  the  calculated  ambient  burning  rate  due  to  the  way  that  the 
reference  burning  rate  is  incorporated  into  the  model.  Within  the  scatter  of  the  data, 
any  of  the  three  values  would  be  reasonable.  The  value  of  15  was  used  previously  in 
the  mode|45  as  this  is  the  value  determined  experimentally  by  Pearson,  et  ai^s,  and  it 
is  also  the  value  calculated  by  Guirao  and  WilliamsSO. 
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Three  sources  of  surface  temperature  data  for  AP  were  found  in  the  literature 
and  used  for  a  basis  of  comparison.  Three  sources  of  linear  pyrolysis  hot  plate 
data^t  .42,51  using  burning  pressed  pellets  of  AP  are  compared  with  the  extensive  data 
of  Powling40  using  an  IR  detector  with  AP  composite  propellants.  These  data  are 
compared  with  the  calculated  surface  temperatures  of  Beckstead  and  Hightower^s 
based  on  the  cubic  phase  thickness  of  quenched  single  crystals  of  AP.  Figure  23 
shows  the  combined  data  plotted  as  log  rate  versus  the  reciprocal  of  the  surface 
temperature,  compared  with  the  modeling  calculations.  There  appears  to  be  a 
consistent  bias  in  the  linear  pyrolysis  data  relative  to  the  IR  detector  data.  The  hot 
plate  data  consistently  show  higher  surface  temperatures.  Variations  in  the  surface 
activation  energy  indicate  that  a  value  between  ~22  and  30  Kcal/mole  would  fit  the 
data.  The  calculated  surface  temperature  values  from  the  Price,  Boggs,  Derr 
mode|47,48  have  also  been  plotted  for  reference.  Their  calculations  agree  well  with 
the  surface  temperatures  of  Beckstead  and  Hightower.  It  is  also  very  apparent  that 
most  of  the  data  were  obtained  at  burning  rates  and  surface  temperatures  well  below 
those  of  normal  interest.  Again,  within  the  scatter  of  the  data,  any  of  the  calculated 
surface  temperature  curves  would  appear  to  be  reasonable.  The  need  for 
experimental  surface  temperature  for  rates  above  0.1  cm/sec  is  also  very  apparent. 


1/Ts  (1/K) 


Figure  23.  Comparison  of  AP  surface  temperature  data  and  model  calculations. 


For  both  AP  and  HMX  a  significant  amount  of  reaction  apparently  occurs  in  the 
melt  layer  (up  to  75%)27.45,50.  jhe  calculated  heat  release  is  shown  in  Figure  24. 
Although  there  are  no  experimental  data  for  the  heat  release  in  AP,  the  calculated 
values  are  very  comparable  to  measured  values  for  double  base  propellants43  and 
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HMX  monopropGlIant^^’^'^i  In  the  original  BDP  monopropellant  calculations,  values  of 
-120  and  -150  were  used  for  consistency.  The  values  calculated  with  the  value  of 
of  0.20  are  in  this  same  range  and  therefore,  a  value  of  Qq  equal  to  0.20  was  used  in 
the  model. 


20  50  100  atm 


Figure  24.  Calculated  condensed  phase  heat  release  (values  are  exothermicl- 


Another  variable  of  interest  is  the  flame  standoff  distance.  Figure  25  is  a  plot  of 
the  flame  standoff  distance  calculated  by  the  model.  The  flame  distances  are 
comparable  to  those  calculated  by  other  models.  The  calculated  values  for  the 
original  BDP  mode|27.45^  from  the  Price,  Boggs,  Derr4T48  model  and  the  Guirao  and 
Williams  model^o  have  all  been  included  for' comparison.  At  low  pressure  there  is  a 
considerable  discrepancy  between  the  different  calculations  but  at  increasing 
pressures  the  BDP,  SST  and  PBD  models  all  come  together.  The  flame  standoff 
distance.  Equation  (6),  is  very  dependent  of  the  value  of  the  thermal  conductivity 
assumed.  Guirao  and  Williams  used  a  thermal  conductivity  of  0.0001  while  a  value  of 
0.0003  has  been  used  in  the  BDP  modeling.  This  factor  of  three  accounts  for  most  of 
the  difference.  In  the  Price,  Boggs,  Derr  model  a  variable  thermal  conductivity  value  is 
used  as  they  integrate  through  the  flame.  Their  calculations  would  appear  to  be  the 
most  accurate.  For  the  simplified  model  used  here  with  a  constant  thermal 
conductivity,  the  standoff  distances  calculated  appear  to  be  reasonable.  Again,  it  is 
unfortunate  that  there  are  no  experimental  data  for  comparison  purposes. 
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Figure  25.  Calculated  flame  stand-off  distance  for  AP  monoDropellant. 


AN  Monopropellant  Calculations 

Monopropellant  parameters  have  been  incorporated  into  the  model  to  calculate 
the  AN  monopropellant  burning  (see  Table  4).  The  burning  rate  has  been  calculated 
as  a  function  of  pressure  and  compared  to  Whittaker’s  data^s  for  pressed  pellets  of  AN 
containing  up  to  2.5%  CraOa.  The  calculated  burn  rates  are  compared  in  Figure  26, 
and  the  calculated  values  of  ap  are  shown  in  Figure  27.  Using  the  available  physical 
properties  for  AN,  the  following  combustion  parameters  (see  Table  4)  were  selected:  a 
reference  burning  rate  of  0.19  cm/sec  at  68  atm,  a  reference  surface  temperature  of 
530  K,  reaction  order  of  2,  Es  of  14  Kcal/mole,  Ef  of  25  Kcal/mole  and  a  Qc  value  of 
0.26.  The  reference  burning  rate  is  used  to  calculate  the  flame  kinetic  prefactor  at  a 
given  condition,  thus  ensuring  that  the  model  will  calculate  the  correct  burning  rate  at 
the  reference  condition.  The  burning  rate  exponent  is  influenced  most  by  the  selected 
reaction  order.  An  order  of  2  gives  the  observed  pressure  exponent  of  slightly  less 
than  one. 
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Unfortunately  there  are  no  experimental  data  available  on  the  burning  rate 
temperature  sensitivity,  ap.  The  model  indicates  that  the  ap  is  very  dependent  on  the 
gas  phase  activation  energy.  Four  different  values  of  activation  energy  were  used  to 
calculate  the  observed  results  in  Figure  27.  With  the  selected  flame  activation  energy 

of  25  Kcal/mole,  a  crp  of  ~0.3  was  calculated.  The  precise  value  of  the  activation 
energy  does  not  have  a  significant  effect  of  the  calculated  ambient  burning  rate  due  to- 
the  way  that  the  reference  burning  rate  is  incorporated  into  the  model.  Although  AP 

and  HMX  have  ap  values  on  the  order  of  0.16  and  0.2%/K  respectively,  it  is  anticipated 
that  AN  would  have  a  relatively  high  ap  due  to  its  lower  surface  temperature  and  low 

burning  rate  (similar  to  double  base  propellants  which  have  ap  values  of  0.2  to  0.5 
%/K). 


Figure  28  shows  the  combined  surface  temperature  data^^-^s  for  AN  plotted  as 
log  rate  versus  the  reciprocal  of  the  surface  temperature,  compared  with  the  modeling 
calculations.  It  is  significant  that  the  burn  rate  catalyst  did  not  appear  to  have  a 
significant  influence  on  the  surface  decomposition  characteristics.  This  seems  to  verify 
that  the  gas  phase  reaction  is  the  controlling  mechanism,  and  the  surface 
decomposition  simply  accommodates  the  heat  flux  from  the  gas.  Variations  in  the 
surface  activation  energy  indicate  that  a  value  between  ~10  and  15  Kcal/mole  would  fit 
the  data  approximately  the  same  (the  four  curves  correspond  to  activations  energies  of 
11,  13, 14,  and  15).  The  measured  and  calculated  surface  temperature  varies  between 
500  and  600  K  for  the  AN  monopropellant.  This  puts  AN  in  the  same  range  of  surface 
temperatures  as  double  base  propellants,  but  significantly  below  AP. 
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Figure  28.  Comparison  of  AN  burning  rate  data  and  model  calculations  for  varying 
s.urface  temperature. 
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AN  is  known  to  melt  readily  at  443  K.  Due  to  such  a  low  melting  point,  it  is 
anticipated  that  there  could  be  a  significant  amount  of  condensed  phase  reaction 
occur  in  the  melt  phase,  similar  to  both  AP  and  HMX.  The  calculated  heat  release  is 
shown  in  Figure  29.  Although  there  are  no  known  experimental  data  for  the  heat 
release  in  AN,  the  calculated  values  are  the  same  order  of  magnitude  as  measured 
values  for  double  base  propellants  and  HMX  monopropellant.  Also,  the  calculated 
values  are  very  close  to  the  heat  of  reaction  calculated  for  reaction  (3)  in  Table  3  of  83 
cal/gm  which  corresponds  to  the  reduction  of  nitrogen  to  NO  in  the  condensed  phase. 
At  higher  pressure  and  burning  rate  the  calculated  heat  release  approaches  that  of 
reaction  (4)  in  Table  (3),  109  cal/gm,  corresponding  to  the  reduction  of  nitrogen  to 
N2O.  Therefore,  the  coefficient,  Qc,  value  of  0.29  was  used  in  the  modeling  work. 


Figure  29.  Calculated  condensed  phase  heat  release  (values  are  exothermic). 


Figure  30  is  a  plot  of  distance  parameters  calculated  by  the  model.  The  melt 
layer  thickness  is  calculated  assuming  an  exponential  temperature  distribution  in  the 
condensed  phase  and  an  average  thermal  diffusivity  of  1.6X10-3  for  molten  AN.  The 
calculated  thicknesses  of  15  to  50  microns  over  the  pressure  range  of  interest  are  in 
general  agreement  with  the  observed  window  bomb  movies  of  burning  propellant^^. 
The  calculated  flame  stand-off  distances  of  25  to  250  microns  are  much  larger  than 
those  values  previously  calculated  for  AP  and  HMX.  This  is  due  to  the  lower  burning 
rate  of  AN  (see  Equation  8).  The  flame  distances  are  comparable  to  those  observed 
for  double  base  propellants.  Again,  it  is  unfortunate  that  there  are  no  experimental 
data  for  comparison  purposes. 
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MODEL  EQUATIONS  FOR  AP  AND  AN  COMPOSITE  PROPELLANTS 


Surface  Temperature  Calculation 

In  the  original  BDP  model,  the  surface  temperature  was  averaged  over  the 
oxidizer  as  well  as  the  binder.  To  generalize  and  improve  the  model  the  Beckstead 
SST  mode!  was  developed  with  each  separate  oxidizer  fraction  having  its  own  surface 
temperature  and  allowing  the  binder  to  have  its  own  surface  temperature.  This 
approach  generalizes  the  model  so  as  to  apply  to  propellants  with  very  wide 
distributions  of  AP,  and  allows  application  to  double-base  propellants,  aluminized 
propellants,  and  HMX-oxidized  propellants. 

The  basis  of  the  current  model  is  the  Beckstead  SST  (Separate  Surface 
Temperature)  model^'®’'^^’^^.  The  SST  model  is  an  adaptation  of  the  original  BDP 
model  for  composite  propellants26,27_  but  with  two  major  modifications; 

1)  a  separate  energy  balance  is  written  for  each  ingredient  category, 
resulting  in  a  separate  surface  temperature  calculation  for  each  ingredient,  and 

2)  the  overall  burning  rate  is  calculated  on  a  time  averaged  basis  rather 
than  the  original  BDP  approach  of  using  a  space  averaged  approach. 

These  concepts  have  been  discussed  in  detail  in  previous  publications5.6.44,55 
and  only  the  basic  equations  will  be  repeated  here. 

An  energy  balance  is  written  for  each  type  and  size  of  oxidizer  ingredient.  The 
resultant  equation  for  the  oxidizer  surface  temperature  is 

Tsoxii  =  To  -  +  PpijP  Fij  (1  +■<  /^st)^  e-^‘pFij  +  pAiij ^  e-^  Aiij 

Cp  v-p  '■'P  (9^ 

+  (1-pFii^^2twie-5-MONOi|  -ATloss 

Cp 

where  the  terms  denoted  PF  refer  to  the  primary  flame  and  the  terms  denoted  MONO 
refer  to  monopropellant  flame.  The  various  terms  are  all  defined  in  the  nomenclature. 
The  j  subscript  represents  different  oxidizer  fractions  sizes,  and  i  represents  different 
oxidizer  types  (i.e.  AP,  HMX,  AN  or  SN).  It  should  be  noted  that  the  effect  of  a  final 
flame  does  not  appear  explicitly  in  the  equation,  but  appears  indirectly  through  the 

S  PF  and  S  MONO  terms  which  are  discussed  below. 

Because  of  its  relatively  large  distance  from  the  surface,  the  final  diffusion  flame 
has  a  small  effect  on  the  surface  temperature  (at  low  pressures).  Its  principle  influence 
is  likely  felt  as  an  effective  increase  in  the  monopropellant  flame  temperature. 
Therefore,  the  energy  feedback  from  the  final  flame  was  modified  in  the  SST  model  to 
represent  feedback  from  the  final  flame  to  the  oxidizer  flame  (i.e.,  its  influence  is  to 
raise  the  oxidizer  adiabatic  flame  temperature).  The  mass  of  reactants  in  the  final 
flame  is  determined  from  the  residual  oxidizer  species  in  the  monopropellant  flame 
and  from  the  residual  fuel  species  coming  from  the  primary  flame.  When  AP  and  AN 
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burn,  the  combustion  products  are  not  fully  oxidized  and  -34%  and  14%  oxygen  is 
produced  from  the  monopropeilant  flame.  This  is  available  for  further  reaction  with  the 
fuei  in  the  final  diffusion  flame.  The  fraction  of  the  oxidizer  flow  available  for  final 

reaction  is  denoted  as  Pap  and  has  the  values  of  0.34  for  AP,  0.14  for  AN,  and  0.0  for 
HMX.  The  oxidizer  will  react  with  the  stoichiometric  amount  of  fuel  from  the  binder 
decomposition  products.  Thus,  the  expression  for  the  effective  oxidizer  flame  energy 
becomes: 


OaPIj  =  (QAP)adij  +  pAPij  (1  +  1/<E>St)  QpFije'^''''*' (10) 


where  SFHjSAP  is  the  nondimensional  distance  between  the  two  flames.  This  term 
has  a  minor  effect  on  the  overall  burn  rate  of  propellants  at  low  pressures,  but  has  an 
increasing  effect  above  1000  psi. 


The  equation  for  the  surface  temperature  of  the  binder  is  developed  in  a  manner 
analogous  to  the  oxidizer  surface  temperature5.6.55 


Tsf  —  To  - 


X  {"*  'PP'j)  P  F'i  ^  4"’  +^St)  e'^’pFij 


_  rriAi  Qlm 
rrifSf  Cp 


+ 1 


(11) 


ij 


rrifSf  Cp 


For  a  composite  propellant  the  energy  feedback  from  the  double  base  binder  flame  is 
zero  and  the  equation  reduces  to: 


Tsf  =  To  -  +  X  ("*  'PP'i)  P  FijS  fij{l  +^st)  ^7^  e'^VFij 

Cp  ij  Cp 

■  +  +  XAT.ss:i 


QpR 


(12) 


mfSf  ^p 


rrifSf  Cp 


Equations  (9)  and  (12)  are  the  backbone  of  the  model  and  determine  the  surface 
temperature  of  each  ingredient.  Once  the  surface  temperatures  of  the  various 
ingredients  are  determined,  the  mass  burning  rate  of  each  ingredient,  oxidizers  and 
binder,  is  calculated  from  a  zero-order  Arrhenius  expression  (see  Equation  (1)). 


mj  =  Aje(’E/RT5j) 


The  overall  propellant  burn  rate  is  calculated  from  the  individual  mass  burning  rates, 
appropriately  weighted.  Other  aspects  of  the  model  and  approaches  for  weighting  the 
individual  rates  are  summarized  in  the  following  sections. 
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The  surface  temperature  equations  are  based  on  the  assumption  of  steady  state 
conditions.  However,  for  very  small  particles,  it  is  very  unlikely  that  a  true  steady  state 
condition  is  achieved.  Fiqure  7  of  Volume  III  of  reference  44  compares  the  rates  of 
monomodal  propellants  to  the  characteristic  thermal  wave  thickness  relative  to  the 

binder  thickness  between  particles.  In  general,  for  particles  greater  than  20  to  SOp,  the 
binder  will  achieve  a  steady  thermal  profile.  For  smaller  particles  (depending  on  the 
overall  rate)  the  binder  will  be  preheated  to  75%  of  its  surface  temperature.  This 
represents  a  significant  heat  loss  from  the  small  particles  to  the  binder.  Based  on  a 
steady  thermal  profile  in  the  oxidizer  as  it  burns  into  binder,  the  loss  is  (in  degrees  of 
temperature). 


AT.ss  =  ^(T.x-To)j2^e:^  (14, 

For  particles  above  50p  this  heat  loss  will  be  very  small,  but  for  fine  particles,  the  loss 
can  reach  several  hundred  degrees  resulting  in  a  very  reduced  rate.  It  is  felt  that  this 
mechanism  is  very  significant  in  interpreting  self  extinguishing  phenomena's  and  the 
so  called  "ledge"  mechanism  discussed  in  Reference  44  where  burning  rate  can 
actually  decrease  with  decreasing  particle  size  rather  than  increase. 

The  flame  standoff  distances  that  are  needed  in  Equations  (9)  and  (11)  are 
calculated  as  in  the  BDP  and  SST  models,  i.e. 

^*PF  =  CpHTr  (xVd  +  x*pf)  (‘•^) 

for  the  primary  flame,  and 

^*MONO=  (CpPHox/^)  X*MONO  (16) 


for  the  monopropeilant  flame  (see  also  Equations  (4)  and  (6)). 

Rate  Averaging  Basis 

In  the  BDP  approach  the  local  rate  is  taken  as  the  oxidizer  rate,  averaged  over 
time,  which  requires  that  the  overall  propellant  stoichiometry  by  maintained.  Condon, 
Click,  and  King  define  the  local  rate  as  the  rate  of  a  "pseudo  propellant"  rather  than  the 
oxidizer  rate  as  in  the  original  BDP  approach.  King"''^  has  recognized  what  he  calls 
the  "end  game"  problem  in  determining  a  time  averaged  local  rate.  That  is,  one  must 
account  for  left  over  binder  after  an  oxidizer  crystal  burns  out.  In  the  BDP  approach 
this  problem  was  recognized  but  avoided  by  requiring  that  the  overall  stoichiometry  be 
maintained.  This  constraint  was  assumed  to  be  valid  because  oxidizer  effects 
dominate  the  burn  rate  of  AP-composite  propellants:  the  binder  having  a  secondary 
influence. 

The  average  propellant  burning  rate  is  a  summation  of  local  rates.  However, 
the  binder  rate  is  radically  different  from  the  oxidizer  rate  (see  Figure  18).  This 
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difference  seriously  effects  if  not  invalidates  the  typical  averaging  process.  If  a  time 
integration  is  performed  using  the  space  averaged  rate,  it  is  discovered  that  the 
oxidizer  will  invariably  burn  out  before  the  binder  does,  leaving  residual  binder. 
Accurate  modeling  requires  an  accounting  of  that  residual  binder  and  its  influence  on 
the  overall  rate.  To  date  all  models  other  than  the  SST  model  have  ignored  this. 

The  conventional  averaging  process  depends  on  two  critical  assumptions: 

1.  the  local  rate  is  the  rate  of  a  pseudo  propellant  including  oxidizer  and  binder; 

2.  the  oxidizer  dominates  the  overall  propellant  rate. 

If  either  of  these  assumptions  are  not  satisfied  when  the  conventional  approach  breaks 
down.  Based  on  the  results  of  Beckstead''3-26  and  Condon20,25_  \\  appears  that  these 
assumptions  are  reasonable  for  a  large  number  of  AP/HTPB  propellants.  However, 
they  obviously  do  not  apply  to  double  base  propellants  or  HMX/composite  propellants. 
Also  based  on  the  results  summarized  in  Reference  44,  it  appears  that  the 
assumptions  fail  for  aluminized  propellants  and  for  propellants  with  very  wide  oxidizer 
distributions.  In  other  words,  it  appears  that  wherever  the  rate  of  the  binder  has  a 
significant  influence  on  the  propellant  rate,  the  conventional  assumptions  will  fail. 

In  the  SST  model^.s  the  equation  for  the  average  pseudo  propellant  burn  rate  is 
of  the  form 


r.  -  2  distances  _  Z  distances 

Z  times  Z  (distances/rates  +  delay  times) 


(17) 


The  summations  for  distances  and  times  are  performed  over  the  burning  path,  and 
therefore,  relate  to  characteristic  oxidizer  and  binder  dimensions.  This  allows  the 
binder  to  have  a  significant  influence  on  the  rate.  Instead  of  a  single  local  rate,  both 
oxidizer  and  binder  rates  are  accounted  for.  Also,  for  each  burning  time  associated 
with  a  given  dimension,  there  is  an  associated  transient  time.  In  the  past  these  delays 
or  transient  times  have  often  been  neglected.  If  they  are  neglected,  the  time  averaging 
approach  can  be  similar  to  surface  averaging.  If  the  delay  times  are  properly  included 
the  results  are  dramatically  different. 

The  ignition  delay  time  is  the  time  it  takes  an  oxidizer  crystal  to  ignite,  assuming  that 
the  ignition  process  starts  when  a  crystal  first  becomes  exposed  to  the  burning  surface. 
As  the  burning  surface  sweeps  over  the  crystal,  heat  is  fed  to  the  crystal,  which  heats 
up,  begins  decomposing,  and  finally  ignites.  The  time  required  for  this  process  to 
occur  is  the  ignition  delay  .time.  It  has  been  included  in  previous  models  based  on  an 
empirical  data  correlation.  An  improved  correlation  that  is  now  being  used  has  been 
developed  and  reported  in  previous  reportsSS, 

The  delay  time  associated  with  the  binder  is  a  binder  burn  through  time.  It  is  related 
to  King's  "end  game",  attempting  to  account  for  unburned  binder.  After  a  crystal  burns 
out,  a  pocket  of  binder  remains  and  must  be  consumed  by  pyrolysis,  sluffing,  or  some 
other  mechanism.  An  approach  to  describing  this  phenomenon  is  discussed  later. 
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In  a  burning  propellant,  as  the  surface  approaches  a  particle  of  AP  or  AN,  the 
particle  will  be  heated  to  a  temperature  characteristic  of  the  binder  pyrolysis  rate  (i.e. 
the  medium  in  which  it  is  immersed).  Thus,  when  the  particle  is  exposed  to  the 
propellant  burning  surface,  its  surface  temperature  will  be  that  of  the  binder,  not  that  of 
a  steady  burning  particle.  Because  both  AP  and  AN  have  different  surface 
temperatures  than  HTPB  at  a  given  rate,  .when  they  are  first  exposed  to  the:  surface 
they  may  not  ignite  immediately.  Each  particle  will  be  exposed  to  the  surface 
conditions,  heat  up  enough  to  reach  ignition,  ignite,  go  through  a  short  transient  time, 
and  ultimately  reach  a  steady  state  burning  rate,  dependent  on  the  surroundings.  The 
SST  model  attempts  to  account  for  the  discrepancies  in  surface  ternperature  and  the 
associated  transient  ignition  heatup  and  delay  times.  The  SST  equation  for  the  overall 
propellant  burning  rate  is 


£NiDi'(l  +25i) 


r= 


I  Ni  d; 


1  ■  ^igni 

roxi  □; 


+ 


Tb 


(18) 


where: 

Ni  = 

Di  = 

5i  = 

figni  = 
tbi  = 
roxi  = 
rb  = 


number  of  ith  oxidizer  crystals  per  unit  length 
statistical  intersection  diameter  of  ith  oxidizer  size 
distance  from  oxidizer  to  center  of  adjacent  binder 
ignition  delay  time  of  ith  oxidizer  size 
binder  burnthrough  time  of  ith  oxidizer  size 
rate  of  ith  oxidizer 
binder  rate 


The  four  terms  in  the  denominator  represent,  respectively,  1)  the  time  to  burn 
through  an  oxidizer  particle,  2)  the  ignition  delay  time  (or  transient  heatup  time)  of  the 
corresponding  oxidizer  particle,  3)  the  time  associated  with  burning  the  binder 
associated  with  the  oxidizer  particle  and  4)  the  transient  time  associated  with  burning 
through  any  remaining  unburned  binder. 


Multimodal  and  Mixed  Oxidizer  Propellants 

A  significant  problem  in  modeling  multimodal  and  mixed  oxidizer  propellants  is 
accounting  for  the  effects  and  influences  caused  by  the  fact  that  different  sizes  and 
types  of  particles  can  and  do  have  significantly  different  burn  rates  (i.e.  interaction 
effects-the  particles  do  not  burn  independently).  What  undoubtedly  happens  is  that  f 
finer  particles  along  the  periphery  of  a  coarse  particle  burn  faster  than  the  coarse 
particle.  The  net  result  is  that  the  overall  propellant  rate  is  much  closer  to  the  rate  of 
the  pseudopropellant  corresponding  to  the  fines  rather  than  that  of  the  coarse.  The 
same  thing  can  happen  when  two  components  with  significantly  differing  rates  burn 
near  each  other.  This  can  be  visualized  in  two  ways.  Either  the  faster  burning 
components  are  igniting  the  slower  burning  particles  all  along  their  periphery, 
eliminating  the  normal  ignition  delay  time  and  coupling  with  their  burn  rate,  or  the 
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slower  burning  particles  can  be  visualized  as  being  undercut  by  the  faster  burning 
components  which  "short  circuit"  the  effect  of  the  slower  burning  particle.  The  problem 
is,  how  to  model  this  coupling  process. 

Within  the  computer  program  the  rates  of  individual  size  fractions  are  monitored 
relative  to  the  rate  of  the  finest  fraction  (assumed  to  have  the  highest  rate).  In  this 
manner  the  fraction  of  the  coarser  oxidizer  sizes  that  is  effectively  contributing  to  the 
overall  rate  is  calculated  as 


7C_  Dpi 

froxi  =  .^Amallest.  (^9) 

^  +  tigni 
•0x1 

The  fraction  of  oxidizer  term  multiplies  the  1/rox  term  in  Equation  (18)  when  HMX  is  the 
oxidizer  in  a  double  base  propellant,  and  multiplies  the  1/rox  and  tjgn  terms  for  the 
coarse  fractions  of  multimodal  propellants  as  discussed  below.  Equation  (19) 
accounts  for  the  interaction  influence  between  particles  of  drastically  differing  sizes. 
To  allow  for  the  potential  interaction  of  different  types  of  oxidizers,  the  computer 
program  has  been  set  up  to  allow  the  AP  to  have  precedence  over  the  other  oxidizers. 


Oxidizer  Ignition  Delay  Time 

In  most  models,  varying  the  ignition  delay  time  parameters  has  little  effect  on 
the  net  burning  rate  due  to  the  way  the  term  enters  into  the  Sox  calculation.  A 
fundamental  approach  was  developed  as  part  of  the  SST  model  to  calculate  the 
ignition  delay  time  of  a  crystal  in  a  propellant.  A  linear  regression  analysis  of 
Shannon's  data^^  for  AP  was  performed  and  the  ignition  delay  time  has  been 
calculated  from  the  following  equation's 

tign  =  Cjgn  (20) 


Equation  (20)  has  been  programmed  into  the  computer  program  and  is  being  used  for 
all  oxidizer  calculations. 

The  ignition  delay  time  associated  with  AN  is  not  clear  because  of  the  melting 
characteristics  of  AN.  Under  normal  conditions  AN  melts  at  443  K.  In  a  combustion 
environment  window  bomb  photography  shows  AN  melting  on  the  surface  of  burning 
propellants.  Thus,  the  AN  particles  may  start  to  lose  their  individual  size  and  shape 
characteristics.  Very  small  particles  may  melt  in  their  entirety,  particularly  at  low 
pressures  or  for  low  burning  rate  propellants.  Large  particles  may  maintain  their 
integrity  and  burn  with  the  thin  molten  layer  on  their  surface.  It  is  not  clear  that  for  the 
times  associated  with  combustion,  whether  the  molten  AN  is  able  to  move  enough  in  a 
lateral  direction  to  mix  with  the  surrounding  binder  significantly.  Window  bomb 
photography32  indicates  that  the  lateral  flow  of  molten  AN  may  occur.  The  general 
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melting  characteristics  of  AN  would  prohibit  AN  crystals  from  protruding  above  the 
surface  and  would  thus  reduce  the  effects  of  an  ignition  delay  concept.  Because  of 
these  reasons  the  coefficient  on  the  ignition  delay  time  was  reduced  by  a  factor  of  four 
below  that  used  for  AP54. 


Binder  Burnthrouah  Time 

In  wide  distribution  propellants  the  binder  is  distributed  disproportionately  so 
that  the  fines  are  very  fuel  rich  and  the  coarse  is  oxidizer  rich.  The  fines,  being 
surrounded  in  their  'sea  of  binder',  are  thermally  insulated  from  the  rest  of  the 
propellant.  The  binder  absorbs  heat  from  the  fines,  reducing  their  rate.  When  they 
burn,  they  do  so  quickly  and  leave  large  amounts  of  binder,  which  pyrolyzes  slowly, 
thus  impeding  the  overall  rate.  As  the  percentage  of  fines  increases  or  the  size 
decreases,  this  problem  is  compounded,  leaving  more  binder  unburned  and 
magnifying  the  impeding  effect.  If  the  size  of  the  fines  is  increased  or  the  percentage 
reduced,  the  propellant  will  burn  more  efficiently  and  the  rate  will  be  more  predictable. 

Oxidizers  that  are  not  rich  in  oxygen,  such  as  HMX  or  AN,  will  undergo  a  similar 
phenomenon.  If  the  oxidizer  cannot  consume  ail,  or  a  large  majority  of  the  binder 
surrounding  an  individual  crystal,  then  the  remaining  binder  must  somehow  pyrolyze 
without  the  aid  of  an  oxidizing  species. 

This  concept  shows  up  in  the  model  in  two  ways:  (1 )  the  loss  of  energy  from  the 
oxidizer  to  the  binder  (AT|oss  term  in  the  surface  temperature  equations),  and  (2)  the 
binder  burnthrough  time.  The  time  for  binder  to  cook  off  (after  a  crystal  has  burned 
past  it)  should  be  related  to  the  thickness  of  the  heated  binder.  If  the  thermal  wave 
penetrates  deeper  than  the  interparticle  distance,  (as  is  the  case  for  small  particles) 
the  preheated  binder  will  be  more  prone  to  pyrolyze.  The  burnthrough  time  should 

also  be  proportional  to  the  thermal  response  time  of  the  binder  (cxt/rb^)-  Thus,  if  the 
response  time  is  short  (i.e.  if  the  binder  can  conduct  heat  more  rapidly  than  it  burns 
away),  then  the  burnthrough  time  ought  to  be  short  also.  The  equation  for  the  binder 
burnthrough  time  in  the  SST  model  was  taken  as 

tb  =  f  (^)  (£)  (21) 

fb  To  0 


Where  0/F  is  the  oxidizer  to  fuel  ratio,  at  is  the  thermal  diffusivity  and  rb  is  the  binder 
rate.  The  binder  rate  in  the  denominator  causes  the  term  to  be  significant  for 
propellants  that  burn  slowly  or  that  have  an  inherently  slow  binder  rate  such  as  AN 
propellants. 

Calculated  Diffusion  Flame  Heights 

The  classic  work  in  this  area  is  the  Burke-Schumann  analysis^s  published  in 
1928.  Virtually  all  subsequent  work  is  either  based  on  the  Burke-Schumann  analysis 
or  at  least  references  that  work  extensively.  The  Burke-Schumann  equation  was 
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programmed  and  used  to  calculate  the  diffusion  flame  standoff  distance  for  many  of 
the  Series  III  propellantsH  Up  to  fifty  terms  had  to  be  used  in  the  summation  to  obtain 
flame  profile  calculations  for  some  of  the  propellants.  Flame  profiles  for  several 
propellants  were  reported  previously^^. 

The  most  significant  results  of  that  previous  study  are  summarized  as  follows: 

(1)  The  prevailing  concept  of  the  primary  flame  always  closing  over  the  oxidizer 
crystal  is  apparently  incorrect.  Only  the  finest  oxidizer  size  fraction  will  do  so.  The 
flame  shape  for  coarser  fractions  within  a  propellant  will  close  over  the  binder  rather 
than  the  oxidizer. 


(2)  The  flame  heights  are  virtually  independent  of  burn  rate  for  the  range  of 
conditions  covering  most  propellants.  This  observation  can  be  a  great  aid  in 
modeling.  It  infers  that  a  simpler  mathematical  solution  can  be  applied  (rather  than  a 
complete  flame  analysis),  which  reduces  the  complexity  of  the  flame  height  calculation 
significantly. 

(3)  Flame  heights  (or  diffusional  distances)  of  the  oxidizer  rich  flames  over  the 
binder  can  be  very  small  relative  to  the  oxidizer  size.  If  the  kinetics  also  lead  to  narrow 
reaction  zones,  this  could  be  interpreted  as  a  surface  reaction  relative  to  the  binder 
with  little  or  no  heat  feed  back  to  the  oxidizer. 


As  part  of  the  SST  model  a  correlation  was  developed  to  predict  flame  heights, 
x*.  instead  of  using  a  lengthy  computer  program  to  compute  the  heights  each  time  they 
are  needed  within  the  burning  rate  model.  A  general  correlation  for  the  diffusional 
flame  height  as  a  function  of  the  propellant  properties  can  be  expressed  as 


X*  _^Yst^a2^  gfi  \a3 
b  CCoxi^ 


(22) 


where  over  the  binder 
ai  =  0.834 
a2  =  0.858 
as  =  1.871 


over  the  oxidizer 
ai  =  0.472 
a2  =  0.929 
as  =  1.604 


This  expression  gives  a  single  continuous  expression  for  diffusional  distances.  An 
arbitrary  multiplier,  AFH,  is  Included  in  the  computer  program  to  allow  for  uncertainties 
between  the  diffusional  calculations  and  the  primary  flame  calculations.  The  multiplier 
is  normally  set  equal  to  3.0.  This  expression  has  been  programmed  and  is  part  of  the 
computer  model. 


Competing  Flames  •  # 

From  the  general  discussion  of  the  model  and  the  primary  flame  it  is  apparent 
that  there  is  a  range  of  pressures  where  the  monopropellant  flame  and  the  primary 
flame  are  competing  for  the  reactive  oxidizer  specie.  Mathematically  this  is  accounted 

for  by  the  term  Pf  in  the  surface  energy  balances  Equations  (9)  and  (11),  where  Pf  is 
the  fraction  of  the  oxidizing  reactants  that  react  in  the  primary  flame.  At  low  pressures 
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and/or  for  small  particle  sizes,  pp  is  unity.  In  other  words,  the  combined  primary  flame 
mixing  and  reacting  distances  are  less  than  the  monopropeliant  flame  height,  and 
therefore  all  of  the  oxidizing  species  react  directly  with  the  binder  products  rather  than 
with  themselves  in  the  monopropellant  flame.  As  pressure  and  burning  rate  increase, 
the  monopropeliant  flame  height  is  reduced  and  the  primary  flame  reactive  distance  is 
reduced,  but  the  primary  diffusive  distance  is  essentially  held  constant.  Thus,  near  the 
center  of  the  crystal  where  the  diffusive  distance  is  greatest,  the  oxidizing  species  react 

preferentially  in  the  monopropeliant  flame,  and  Pf  is  between  zero  and  one.  At  higher 
pressure,  the  monopropeliant  flame  standoff  becomes  very  small  as  does  the  prirnary 
flame  reacting  distance.  Under  these  conditions  only  a  fraction  of  the  original 
oxidizing  species  react  with  the  binder  products  in  the  primary  flame  near  the  binder- 
oxidizer  interface,  and  Pf  approaches  a  small  number.  These  effects  are  illustrated  in 
Figure  31. 

The  Pf  term  plays  a  dominant  role  in  determining  burn  rates.  The  primary  flame 
is  a  high  energy  flame  and  can  force  the  propellant  burn  rate  much  higher  than  the  low 
energy  monopropeliant  flame.  To  determine  Pf.  the  ideal  situation  would  be  to 
calculate  (or  measure)  the  actual  fraction  of  the  oxidizer  species  that  react  in  the 
primary  flam©  vsrsus  th©  total  availabi©.  This  would  roquir©  th©  accurat©  solution  of 
the  diffusion  equations  in  conjunction  with  the  monopropeliant  flame  equatioiis. 
Unfortunately  a  rigorous  approach  to  do  this  was  considered  beyond  the  scope  of  the 
original  program.  In  retrospect  it  might  have  been  practical  to  attempt  such  a  solution. 
The  characteristics  of  the  Pf  calculation  are  very  dominant  within  the  model  in 
determining  the  burning  rate  dependence  on  pressure  and  particle  size. 

For  flames  closing  over  the  binder,  Pf  is  the  fraction  of  oxidizer  reacting  in  the 
primary  fiame  relative  to  the  totai  oxidizer  available.  The  amount  of  oxidizer  reacting  in 
the  primary  flame  is  proportional  to  the  stoichiometric  ratio,  and  the  total  oxidizer 

available  is  proportional  to  the  local  0/F  ratio.  The  ratio  of  was  used  in  the 

definition  of  Pf  to  incorporate  this  proportionality. 

The  final  equation  to  calculate  pF  is 

Bp  =  XaP _ te:  (23) 

x'pD  +  XpF  0/F 


The  determination  of  x*pd  from  Eouation  22,  assumes  that  all  the  oxidizer  is 
available  for  reaction  in  the  diffusion  flame.  Once  the  monopropeliant  flame  becomes 
active,  it  robs  oxidizing  species  from  th©  primary  flam©  and  a  reduced  amoum  of 
oxidizer  is  then  available.  Actually  the  effective  amount  available  is  the  total  oxidizing 
species  multiplied  by  Pf-  Therefore,  within  the  model  an  effective  0/F  ratio  is 
calculated  by  multiplying  Pf  times  the  overall  0/F  for  each  oxidizer  fraction.  Thus,  the 
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Figure  31 .  Different  flame  configurations  that  can  occur  durino  the  burning  time  of  an 
individual  particle.  Approximate  values  of  Pf  and  pp  that  correspond  to  the  flame 
structure  are  noted. 
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calculation  of  x*pd  and  Pf  is  an  iterative  process.  The  calculation  modifies  the  effective 
geometry  that  is  used  in  the  calculations,  leading  to  smaller  effective  diffusion 
distances  as  Pf  decreases.  This  is  logical  and  should  have  been  incorporated  in  the 
original  BDP  model.  As  the  monopropellant  burns  oxidizing  species,  those  gases  are 
not  available  for  reaction  in  the  primary  flame.  Those  gases  that  are  left,  are  only  over 
the  outer  part  of  the  oxidizer  crystal,  and  thus  have  a  shorter  diffusion  path  to  reach  the 

binder  species.  Therefore,  the  actual  diffusion  distance  is  dependent  on  Pf  and  on  the 
amount  of  oxidizer  reacting  in  the  monopropellant  flame.  Incorporating  this  approach 
into  the  x*PD  determination  leads  to  an  effective  diffusion  distance  which  was  used  in 
the  model. 

Heat  Transfarfrom  The  Primary  Flame 

Once  the  amount  of  oxidizer  that  participates  in  th  primary  flame  is  determined, 
then  the  amount  of  the  heat  that  is  transferred  to  the  oxidizer  versus  the  amount 
transferred  to  the  binder  must  be  determined.  This  effect  is  accounted  for  by  the  term 
pP  in  Equations  (9)  and  (11),  where  pp  is  the  fraction  of  heat  transferred  from  the 
primary  flame  to  the  oxidizer.  When  a  particle  initially  ignites  and  the  flame  closes 
over  the  oxidizer,  as  shown  in  Figure  31 ,  the  value  of  Pp  should  be  greater  than  0.5, 
but  less  than  1.0.  As  the  particle  burns  back  an  increasing  amount  of  the  heat  will  be 
transferred  to  the  binder,  particularly  when  the  flame  closes  over  the  binder  rather  than 
the  oxidizer.  During  this  time  the  value  of  pp  should  be  relatively  small  (but  greater 
than  0.0).  This  is  obviously  a  transient  effect  that  is  not  accounted  for  within  the 
framework  of  the  assumptions  leading  to  Equations  (9)  and  (11).  This  general  effect 
has  been  addressed  previously58,59^  but  a  quantitative  descriptions  of  the  transient 
heat  transfer  is  beyond  the  scope  of  the  current  project.  Thus  an  integrated,  average 
value  of  pp  was  estimated  in  a  purely  empirical  manner.  A  functional  form,  based  on 
a  hyperbolic  tangent,  was  programmed  into  the  code  that  would  typically  be  ~0.5  and 
would  never  exceed  one  nor  fall  below  zero. 

Effect  of  Catalysts 

Modeling  the  effects  of  catalysts  is  another  difficult  area.  Previous  efforts  in  this 
area  by  Beckstead  and  Cohen, so  BurnsideST  BecksteadS,  and  Condon  and  Glick62 
have  been  limited.  Burnside  did  little  modeling  but  did  get  valuable  experimental  data 
on  catalysis.  He  concluded  that  the  burn  rate  augmentation  due  to  FeaOs  is  related  to 
the  specific  surface  of  the  catalyst,  pressure  and  AP  specific  surface.  Mechanistically, 
he  postulates  that  the  catalyzed  reaction  is  promoted  by  adsorption  of  reactive  species 
on  the  heterogeneous  catalyst.  Flanigan's  work  so  on  catalytic  mechanisms,  also 
supports  these  ideas.  In  a  series  of  definitive  small-scale  studies,  Flanigan  and  co- 
workersSO  produced  a  detailed  description  of  the  Fe203  catalytic  process.  Similar 
studies  by  MussoS4  support  these  conclusions.  In  these  studies,  it  was  clearly  shown 
that  the  Fe203  catalyzes  decomposition  of  a  vapor  phase  decomposition  product  of 
ammonium  perchlorate.  Since  the  final  flame  is  too  far  removed  to  influence  rate  and 
AP  flame  is  geometrically  inaccessible  to  the  Fe203,  the  most  logical  position  for  the 
catalyst  reaction  is  in  the  primary  diffusion  flame.  This  basic  assumption  has  been 
carried  over  to  apply  to  the  catalysis  of  AN  propellants. 
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The  model  has  been  modified  to  introduce  catalytic  rate  augmentation  by  the 
heterogeneous  mechanism  in  the  primary  flame  reaction.  The  net  rate  is  assumed  to 
be  the  sum  of  the  normal  rate  plus  a  catalytic  effect.  The  equation  for  the  catalyzed 
rate  is  assumed  to  have  the  form  of  a  general  heterogeneous  reaction. 

Teat  =  —  Acat©  -Ecat/RT pf  PL  (24) 

A-j  +  Ccat 

Ccat  =  catalyst  concentration  x  specific  surface  area 
Ai  =  specific  surface  dependent  constant 
Acat  =  prefactor  for  catalytic  reaction 
Ecat  =  activation  energy  for  catalytic  reaction 

Scat  =  reaction  order  of  catalytic  reaction 
TpF  =  temperature  of  primary  flame 

To  evaluate  the  effectiveness  of  this  format,  a  reaction  order  of  1  was  assumed 
to  simulate  adsorption  (5cat  =  -1  would  correspond  to  desorption,  Scat  =  0  would 
correspond  to  surface  reaction):  Ecat  was  taken  to  be  the  same  as  the  basic  primary 
flame  reaction;  and  Acat  was  taken  as  some  number  much  greater  than  Arf,  the 
primary  flame  frequency  factor.  The  prefactor  ratio  (Acat/ApF)  was  adjusted  to  give 
rates  in  the  right  order  of  magnitude.  Ccat  was  taken  as  the  weight  fraction  of  catalyst 
times  the  specific  surface  area  of  the  catalyst  (meter2/gm).  The  constant  Ai  was 
selected  to  give  an  asymptotic  rate  at  ~5  percent  catalyst.  The  net  reaction  rate  for  the 
catalyzed  primary  flame  becomes 

ratepF  "  Appe  -Epf/RT pf  [1  +  — ^ ]  (25) 

Ai  +  Ccat  ApF 


Using  this  equation,  rates  were  calculated  with  the  SST  model®  for  the  1  percent 
FeaOs  propellants  from  Miller's  Series  VI  propellants44.  The  calculated  rates 
compared  reasonably  well  to  the  experimental  data.  Based  on  this  general  approach 
the  same  equations  have  been  applied  to  the  AN  propellants  that  are  being  modeled 
as  part  of  this  program. 
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AP  COMPOSITE  PROPELLANT  RESULTS 

The  primary  objective  of  the  modeling  effort  was  to  model  AN  based  composite 
propellants  with  mixed  oxidizers.  The  separate  surface  temperature  and  the  time 
averaging  approaches  were  used  to  accommodate  the  description  of  several  oxidizers 
within  the  same  model.  As  has  been  outlined  above,  the  model  is  designed  to 
accommodate  AP,  HMX,  AN,  and  SN  in  any  combination.  Because  of  the  very 
extensive  data  base  of  AP  propellants,  the  approach  used  was  to  evaluate  the  model 
by  comparing  the  calculated  results  for  AP  propellants  varying  composition,  pressure, 
particle  size,  etc,  with  known  data.  For  a  given  propellant,  or  family  of  propellants,  the 
various  parts  of  the  model  can  then  be  examined  individually  to  determine  their 
influence  on  the  overall  rate  within  the  given  formulation  or  test  variables. 
Optimization  of  the  unknown  model  parameters  is  then  performed  to  try  to  match  the 
experimental  behavior.  The  various  combustion  properties  such  as  surface 
temperatures,  flame  standoff  distances,  fraction  of  binder  reacting,  etc.  are  all 
monitored.  If  unrealistic  values  are  calculated,  then  the  source  of  the  discrepancy  is 
traced  within  the  model.  That  particular  part  of  the  model  is  then  explored  to  determine 
what  must  be  done  to  produce  more  realistic  calculations. 

The  data  of  Miller,  et  ap4  for  non-metallized  AP/HTPB  propellants  was  used  as 
the  comparative  data  base,  because  of  the  large  number  of  propellants  with  the  sarne 
consistent  binder  and  a  wide  range  of  compositions,  pressure,  particle  size,  and  initial 
temperatures.  The  data  used  were  from  the  Series  III  and  Series  XI  non-metallized 
propellants44.  The  available  data  consist  of  19  propellants  from  the  Series  111  and  17 
propellants  from  the  Series  XI.  These  propellant  families  covered  a  range  of 
compositions  from  monomodal,  to  tetramodal,  from  77%  to  87.37%  AP,  and  from  0.7  to 
400  micron  AP.  Test  data  were  used  for  500, 1000,  and  2000  psi  (3.4,  6.8,  13.6  MPa). 
Eleven  of  the  Series  III  propellants  also  had  temperature  sensitivity  data.  A  total  of  151 
datum  points  were  used  in  the  overall  parameter  optimization.  Based  on  this 
comparison  procedure,  model  parameters  such  as  the  primary  flame  prefactor,  the 
coefficient  in  the  primary  flame  standoff  distance,  the  coefficient  in  the  expression  for 
Pp  or  pp,  etc.,  were  optimized  to  provide  the  best  fit  of  the  model  to  the  data  base. 
After  optimization  47%  of  the  calculated  points  were  within  10%  of  the  experimental 
data  and  74%  were  within  20%,  with  a  correlation  coefficient  of  0.92. 

With  the  optimized  parameters  the  model  was  compared  directly  with  the  Series 
III  data  at  1000  psi  (6.8  MPa)  and  the  results  are  shown  in  Figure  32.  The  agreement 
is  very  reasonable  with  only  a  few  of  the  widest  distribution  propellants  falling  outside 
of  the  10%  predictive  goal.  Datum  points  lying  outside  of  10%  bounds  have  been 
labeled  with  their  identification  number. 
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Experimental  Rate 


Figure  32.  A  comparison  of  calculated  and  experimental  burning  rates  for  Miller's 
Series  III  propellants  at  1000  psi  (6.8  MPa).  Datum  points  Ivina  outside  of  +10%  have 
been  labeled  with  their  identification  number. 

The  correct  prediction  of  the  pressure  exponent  is  a  significant  challenge  in 
modeling.  The  pressure  dependence  is  caused  by  two  main  factors:  (1)  the  reaction 
order  of  the  primary  flame,  and  (2)  the  relationship  between  the  primary  flame  and  the 
monopropellant  flame.  A  second  order  primary  flame  reaction  was  assumed,  as  a 
second  order  gas  phase  reaction  seems  most  consistent  with  high  temperature,  gas 
phase  kinetics.  The  second  factor  contributing  to  pressure  dependence  is  much  more 
difficult  to  evaluate.  The  relationship  between  the  primary  diffusion  flame  and  the 
premixed  monopropellant  flame  is  very  complex.  Within  the  model  this  interaction 
enters  into  the  equations  through  two  terms,  Pp  and  Pp.  Pp  represents  the  fraction  of 
the  oxidizer  decomposition  products  that  react  in  the  primary  diffusion  flame.  The  Pp 
term  represents  the  fraction  of  energy  fed  back  from  the  primary  flame  to  the  oxidizer. 
The  exact  nature  of  these  terms  and  how  they  vary  with  propellant  composition  and 
test  variables  is  an  area  that  is  not  well  understood,  but  still  needs  extensive  study  in 
the  future.  ' 

The  pressure  exponent  calculations  are  compared  with  the  Series  ill  data  at 
1000  psi  (6.8  MPa)  in  Figure  33.  In  the  calculations  that  were  made  the  pressure 
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exponent  was  consistently  low  for  all  of  the  propellants.  For  most  of  the  propellants 
simulated,  the  pressure  exponent  was  within  20%  of  the  measured  values,  but  several 

datum  points,  particularly  for  the  propellants  containing  400|i  AP  were  more  than  20% 
in  error.  The  very  wide  distribution  propellants  containing  400jj,  AP  have  been 
particularly  difficult  to  model  in 'the  past  also7-8.i3,i7^  due  to  the  extreme  variation  in 
the  particle  sizes. 


Experimental  Exponent 

Figure  33.  A  comparison  of  calculated  and  experimental  pressure  exDonents_fcr 
Miller's  Series  III  nrooellanfs  at  1000  osi  (6.8  MPal  Datum  points  lying  outside_of 
+20%  have  been  labeled  with  their  identification  number. 

The  model  was  also  used  to  calculate  the  temperature  sensitivity  of  the  Series 
III  propellants  where  there  were  data  available.  The  calculations  comparing  the 
calculated  values  to  the  experimental  values  are  shown  in  Figure  34.  The  agreement 
is  very  reasonable  with  only  three  calculations  lying  outside  of  the  20%  range  and  two 

of  those  calculations  were  for  propellants  containing  400u.  AP.  The  calculations  are 
also  compared  with  data  in  Figure  35  where  the  temperature  sensitivity  has  been 
plotted  versus  the  size  of  the  fine  AP  fraction  for  propellants  containing  400ii  and  200^i 

coarse  AP.  Other  than  the  calculation  for  400)1  coarse  and  20|i  fine  AP  the  calculated 
trends  are  in  the  right  direction,  but  the  model  underpredicts  the  observed  trend  in  the 
data. 
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Figure  34.  A  comparison  of  calculated  and  experimental  temperature  sensitivity  for 
Miller's  Series  III  propellants  at  1000  psi  (6.8  MPa).  Datum  points  Ivina  outside  of 
+20%  have  been  labeled  with  their  identification  number. 


Figure  35.  Calculated  temperature  sensitivity  for  Miller's  Series  III  propellants  at  1000 
Dsi  (6.8  MPat.  compared  with  experimental  data  as  a  function  of  fine  AP  fraction  size. 
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AN  COMPOSITE  PROPELLANT  RESULTS 

The  primary  objective  of  the  program  was  to  model  AN  based  composite 
propellants  with  mixed  oxidizers.  As  has  been  outlined  above,  the  model  is  designed 
to  accommodate  AP,  HMX,  AN,  and  SN  as  the  mixed  oxidizers.  The  unknown  model 
parameters  were  evaluated  by  comparing  the  calculated  results  for  AP  propellants 
with  known  data  varying  composition,  pressure,  particle  size,  etc.  Those  parameters 
that  are  anticipated  to  be  independent  of  the  chemical  nature  of  the  oxidizer  have 
been  frozen  in  the  model  with  the  values  determined  from  the  AP  propellant 
optimization.  These  include  the  the  functional  form  of  the  equations  describing  the 
terms,  pp,  Pp.  x*,  etc,  and  the  coefficients  and  limits  on  those  parameters.  On  the  other 
hand,  parameters  such  as  the  primary  flame  activation  energy  and  prefactor  will 
obviously  vary  with  oxidizer  type.  Because  there  is  not  a  consistent  set  of  data  for  AN 
propellants,  such  as  that  available  for  AP  propellants,  the  unknown  parameters  were 
determined  by  comparing  calculated  trends  with  what  appeared  to  be  reasonable. 

AN/HTPB  Propellant 

AN/HTPB  propellants  can  be  expected  to  have  a  much  lower  burning  rate  than 
AP/HTPB  propellants.  This  is  because  of  the  lower  flame  temperature  of  the  AN 
propellants.  The  flame  temperature  calculations  presented  in  Figure  12  show  that  at 
an  86%  loading  an  AP  propellant  has  a  primary  flame  temperature  of  2700  K,  while  an 
AN  propellant  has  a  flame  temperature  of  1700  K.  For  an  activation  energy  of  15 
Kcal/mole,  such  as  used  in  the  model  for  the  primary  flame  activation  energy,  a  1000  K 
change  in  flame  temperature  should  cause  a  change  in  rate  of  a  factor  of  five.  Thus,  if 
a  typical  AP  propellant  has  a  rate  of  0.5,  then  a  corresponding  AN  propellant  will  have 
a  rate  of  ~0.1.  This  is  approximately  the  ratio  of  burning  rates  that  are  observed  for  the 
two  types  of  propellants. 

Most  AN  based  propellants  have  catalysts  or  additives  of  some  type  to  augment 
the  rate  of  the  propellant.  For  standardization  of  the  model  however,  it  is  very 
desirable  to  have  a  propellant  with  no  additives  to  provide  a  basis  of  comparison.  One 
such  propellant  developed  by  Talley  Corporation's  was  identified  at  the  AFAL 
sponsored  workshop  on  AN  Combustion  held  in  November  1987.  That  propellant 
contained  77.5%  AN  with  an  HTPB  binder.  Its  burning  rate  at  1000  psi  (6.8  MPa)  was 
0.47  in/sec  (0.119  cm/sec),  and  it  was  reported  to  have  an  exponent  of  0.59.  Using 
these  characteristics  as  a  standard  a  series  of  calculations  were  made  calculating  the 
burn  rate  as  a  function  of  pressure  for  varying  AN  content.  The  AN  particle  size  was 

not  known,  and  was  assumed  to  be  ~50|i.  Parametric  calculations  indicated  that  the 
particle  size  has  little  effect  on  the  burn  rate.  The  results  are  presented  in  Figure  36. 
The  rates  at  1000  psi  (6.8  MPa)  vary  from  -0.05  to  0.1  (0.1  to  0.2  cm/sec),  which 
seems  very  reasonable.  The  data  at  the  reference  pressure  have  also  been  plotted  in 
Figure  37  as  both  burn  rate  and  pressure-  exponent  versus  AN  content.  The  burn  rate 
increases  almost  linearly  with  AN  content  which  also  seems  very  reasonable.  The 
pressure  exponent  is  higher  that  the  datum  point  of  0.59  at  78%  and  the  exponent 
reaches  a  maximum  at  80%  AN,  which  is  probably  not  consistent  with  data.  The 
reason  that  the  model  predicts  the  maximum  is  not  totally  clear. 
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An  unknown  parameter  of  great  significance  within  the  model  is  the  primaiy 
flame  activation  energy.  For  AP  and  HMX/HTPB  propellants  a  value  of  1 5  Kcal/mole  is 
normally  used  as  the  primary  flame  activation  energy,  based  on  flame  kinetics  work 
with  ammonia-perchloric  acid  flames^s  and  previous  propellant  modeling  work. 
Parametric  calculations  were  made  to  evaluate  the  effect  of  the  primary  flame 
activation  energy  for  AN  propellants.  Figure  38  shows  the  calculated  results  varying 
the  activation  energy  from  values  of  10  to  15  Kcal/mole.  In  each  case  the  primary 
flame  prefactor  was  adjusted  to  give  the  burning  rate  of  the  Talley  propellant  at  78% 
AN.  The  effect  of  the  primary  flame  activation  energy  is  to  cause  the  burn  rate  to 
increase  with  increasing  activation  energy  as  the  AN  content  (and  the  corresponding 
primary  diffusion  flame  temperature)  increases.  An  activation  energy  of  15  seemed  to 
give  rates  higher  than  would  be  expected,  while  a  value  of  10  seemed  to  give  more 
reasonable  results  and  was  used  in  all  of  the  AN  calculations. 


%  AN 

Figure  38.  Thfi  effect  of  the  primary  flame  activation  energy  on  burn  rate  for  varying 
AN  content. 

The  calculations  were  extended  to  calculate  the  temperature  sensitivity  of  the 
same  family  of  propellants.  The  temperature  sensitivity  results  are  shown  in  Figure  39. 
The  calculated  temperature  sensitivity  decreases  significantly  with  increasing  AN 
content  and  with  decreasing  activation  energy.  With  no  precise  data  to  compare  with 
the  values  predicted  with  an  activation  energy  of  10  again  appeared  to  give 
reasonable  values  of  temperature  sensitivity,  particularly  at  the  lower  concentration 
levels.  The  values  of  temperature  sensitivity  are  quite  high  compared  to  conventional 
AP  propellants,  but  based  on  limited  data^s  for  AN  propellants  where  typical  values  of 
7tK  values  of  0.2  to  0.3%/F  were  quoted,  which  are  about  double  those  normally 
expected  for  AP  propellants. 
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The  effect  of  particle  size  was  also  explored.  For  particle  sizes  ranging  from  50 
to  200  microns  there  was  no  significant  difference  in  the  calculated  burning  rates,  at 
the  low  solids  loadings.  At  the  higher  solids  loading  there  is  a  slight  effect  which  will 
be  discussed  later. 


Figure  39.  The  effect  of  the  primary  flame  activation  energy  on  temperature  sensitivity 
for  varying  AN  content. 


Effect  of  Catalysts 


Because  of  the  very  low  inherent  burning  rate  of  AN  propeHants,  catalysts  are 
usually  used  to  increase  the  burning  rate.  To  accommodate  this  effect  within  the 
model,  the  catalytic  activity  is  described  in  Equations  (24)  and  (25),  the  inherent 
assumptions  being  that  the  primary  flame  is  catalyzed  through  the  adsorption  of 
reactants  on  the  catalyst.  The  general  characteristics  of  the  functional  form  of  the 
equation  for  the  catalyzed  primary  flame  were  previously  evaluated  by  comparison  to 
AP  propellants^.  As  a  part  of  an  experimental  task  of  this  program,  a  propellant  (BRM- 
1)  was  made  and  burned,  which  contained  80%  150|i  AN  with  2%  ammonium 
dichromate  (AD)  and  18%  HTPB  binder.  This  datum  point  plus  the  Talley  propellant 
discussed  in  the  previous  section,  gave  two  reference  points  for  comparison  with  the 
calculations. 


Burn  rate  calculations  were  made  for  a  family  of  AN/AD/HTPB  propellants 
varying  the  catalyst  concentration  and  the  percentage  solids.  The  results  are 
presented  in  Figure  40  and  compared  to  the  two  datum  points.  The  calculations  show 
the  rate  increasing  with  increasing  catalyst  concentration  and  with  increasing 
percentage  solids.  The  trends  appear  to  be  consistent  and  reasonable.  The 
calculations  are  in  agreement  with  the  experimental  data. 
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Figure  40  .  Rnrn  rate  calculations  for  AN/AD/HTPB  propellants,  varying  the  catalyst 
concentration  and  the  percentage  solids.  Two  experimental  datum  points  are  included 
for  reference. 

The  calculated  effect  of  the  catalyst  is  not  great.  This  is  due  to  the  inherently  low 
energy  of  AN.  The  primary  flame  temperature  for  the  AN/HTPB  system  is  relatively  low 
(see  Figure  12),  and  no  amount  of  catalysis  is  going  to  increase  the  burning  rate  to 
values  comparable  to  an  AP  propellant  because  of  the  low  flame  temperature  (i.e.  the 
low  energy)  compared  to  the  flame  temperature  of  AP  propellants. 

The  corresponding  pressure  exponent  calculations  for  these  propellants  at 
1000  psi  (6.8  MPa)  are  plotted  in  Figure  41  versus  the  catalyst  concentration.  In  the 
calculations  that  were  made  the  pressure  exponent  consistently  decreases  with 
increasing  catalyst  concentration  and  is  generally  lower  for  the  higher  concentration  of 
solids.  The  decrease  in  exponent  is  a  direct  result  of  the  manner  in  which  the  catalytic 
effect  is  modeled  (see  Equation  (25)).  The  catalysis  is  assumed  to  occur  through 
adsorption  of  the  primary  flame  gases  onto  the  catalyst  surface.  The  adsorption 
process  is  normally  assumed  to  be  proportional  to  the  pressure.  In  comparison,  the 
global  kinetics  of  the  primary  flame  are  assumed  to  have  an  order  of  1 .5.  Thus,  the 
catalysis  has  a  lower  pressure  dependence  than  the  uncatalyzed  process.  The  net 
effect  of  the  catalysis  is  to  reduce  the  pressure  dependence  of  the  overall  reaction. 
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Figure  41.  Calculated  pressure  exponent  for  AN/AD/HTPB  propellants,  varying  the 
catalyst  concentration  and  the  percentage  solids. 

The  corresponding  temperature  sensitivity  calculations  for  these  propellants  at 
1000  psi  (6.8  MPa)  are  plotted  in  Figure  42  versus  the  catalyst  concentration.  In  the 
calculations,  the  temperature  sensitivity  is  virtually  independent  of  catalyst 
concentration  except  for  the  lowest  concentration  of  solids  where  it  increases  slightly. 
The  magnitude  of  the  temperature  sensitivity  appears  to  be  somewhat  high  especially 
when  compared  to  AP  propellants,  but  it  probably  not  unreasonable  for  AN 
propellants.  Anticipating  that  a  typical  propellant  would  probably  have  about  2% 
catalyst  in  its  formulation,  the  zero  and  2%  data  have  been  replotted  in  Figure  43.  The 
decrease  in  temperature  sensitivity  with  increasing  solids  loading  is  dramatic.  This 
trend  is  probably  consistent  also. 
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Calculations  were  also  made  to  evaluate  a  potential  effect  of  particle  size. 
Although  virtually  no  particle  size  effect  was  calculated  for  the  78%  AN  propellant 
discussed  previously,  there  was  an  effect  calculated  for  an  86%  solids  loading 
propellant  with  2%  catalyst.  The  calculated  results  are  shown  in  Figure  44.  The  fine 
size  was  varied  from  10  to  100  microns  and  the  coarse  to  fine  ratio  was  varied  from 
70/30  to  30/70.  Increasing  the  concentration  of  fine  particles  did  in  fact  increase  the 
burn  rate,  but  changing  the  particle  size  had  little  effect.  For  the  lowest  concentration 
of  fines  the  rate  actually  decreased  slightly  with  decreasing  size,  which  is  the  opposite 
effect  from  what  is  generally  observed  with  AP  propellants.  However,  the  AN 
environment  is  much  more  fuel  rich  than  AP  propellants  (i.e.  relative  to  CO  the 
stoichiometric  condition  for  AP  is  85%  AP  while  for  AN  it  is  90%).  Within  the  model  the 
decrease  in  rate  with  decreasing  particle  size  corresponds  to  an  excessive  energy 
loss  to  the  binder  (see  Equation  (21)).  The  corresponding  effect  has  been  seen  with 
very  wide  distribution  AP  propellants  and  in  self  extinguishing  AP  propel lants^®. 
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Figure  44.  Burn  rate  calculations  for  AN/AD/HTPB  propellants,  varying  the 
concentration  and  size  of  the  AN  fine  fraction  at  86%  solids. 


Temperature  sensitivity  calculations  were  performed  for  the  same  family  of 
propellants,  and  the  calculated  results  are  shown  in  Figure  45.  There  is  a  slight  effect 
calculated  for  the  86%  solids  loading  propellant.  Decreasing  the  concentration  of  fine 
particles  did  increase  the  temperature  sensitivity  slightly  for  decreasing  fine  particle 
size,  but  had  the  opposite  effect  for  larger  particles. 
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Figure  45.  Temperature  sensitivity  calculations  for  AN/AD/HTPB  propellants,  varying 
the  concentraticn  and  size  of  the  AN  fine  fraction  at  86%  solids. 

Effect  of  AP  on  AN  Propellants 

A  critical  aspect  in  the  development  of  the  model  is  the  ability  to  handle  mixed 
oxidizers,  specifically  AP  and  AN.  A  series  of  calculations  were  made  for  an  86% 
solids  loading  propellant  with  2%  catalyst,  varying  the  concentrations  and  particle  size 
of  the  AP.  The  AN  was  assumed  to  be  1 50  micron  coarse  with  50  micron  fine  in  equal 
amounts.  The  calculated  results  are  shown  in  Figure  46.  The  AP  size  was  varied  from 
5  to  200  microns,  and  from  0  to  20%  concentration.  Increasing  the  concentration  of  AP 
particles  increased  the  burn  rate,  and  decreasing  the  AP  particle  size  also  increased 
the  burn  rate.  Both  of  these  effects  are  in  agreement  with  experimental  observation. 
The  propellant  with  200  micron  AP  actual  showed  a  slight  decrease  in  the  burning  rate 
with  increasing  concentration.  This  trend  may  not  be  consistent  with  experimental 
data,  but  the  effect  of  200  micron  AP  would  not  be  expected  to  increase  the  rate  very 
much. 


The  pressure  exponent  calculations  for  the  same  family  of  propellants  are 
shown  in  Figure  47.  Increasing  the  concentration  of  fine  AP  particles  increases  the 
exponent,  but  the  exponent  for  the  larger  particles  passes  through  a  maximum  with 
increasing  AP  concentration.  This  trend  could  very  well  be  consistent  with 
experimental  data. 


64 


0.20 


Figure  46.  Burn  rate  calculations  for  86%  solids  AN/AP/HTPB  propellants,  varying  the 
concentration  and  size  of  the  AP. 


Figure  47.  Calculated  pressure  exponents  for  86%  solids  AN/AP/HTPB  propellants, 
varying  the  concentration  and  size  of  the  AP. 
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During  the  experimental  part  of  the  program  three  propellants  corresponding  to 
some  of  the  conditions  of  Figures  46  and  47  were  actually  formulated  and  data 
obtained.  The  data  are  compared  with  the  calculated  results  in  Figure  48.  The 
propellants  were  designated  BRM-275,277,  and  278  and  contained  150  micron 
coarse  and  50  micron  fine  AN  in  equal  amounts.  The  BRM-275  propellant  had  5%  5 
micron  AP;  the  BRM-277  propellant  had  10%  5  micron  AP;  and  the  BRM-278 
propellant  had  10%  200  micron  AP.  The  5  micron  propellant  calculations  are  in  very 
reasonable  agreement  with  the  data,  but  the  200  micron  calculations  are  lower  than 
the  corresponding  data. 


Figure  48.  Oomoarison  of  calculated  and  experimental  burn  rates  for  86%  solids 
AN/AP/HTPB  orooellants.  varying  the  concentration  and  size  of  the  AP. 

Temperature  sensitivity  calculations  were  performed  for  the  same  family  of 
propellants,  and  the  calculated  results  are  shown  in  Figures  49  and  50.  In  Figure  49 
the  calculations  are  plotted  versus  the  AP  concentration  with  particle  size  as  a 
parameter.  The  temperature  sensitivity  is  virtually  constant  with  only  the  200  micron 
AP  propellant  deviating  significantly  from  the  values  between  0.22  and  0.24  %/C. 
Although  there  are  no  experimental  data  with  which  to  make  a  comparison,  the 
magnitude  of  the  calculated  the  temperature  sensitivity  appears  to  be  reasonable.  The 
effect  of  , particle  size  is  illustrated  in  Figure  50  and  is  seen  to  be  small. 
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Figure  50.  AP  particle  size  effe< 
solids  AN/AP/HTPB  oroDellants. 


Effect  of  SN  on  AN  Propellants 


A  significant  aspect  of  the  model  is  the  ability  to  handle  AN  and  SN  as  mixed 
oxidizers.  A  series  of  calculations  were  made  for  an  86%  solids  loading  propellant 
with  2%  catalyst,  varying  the  concentrations  and  particle  size  of  the  SN  similar  to  those 
for  AP.  As  with  the  AP  calculations  the  AN  was  assumed  to  be  1 50  micron  coarse  with 
50  micron  fine  in  equal  amounts.  The  calculated  results  are  shown  in  Figure  51 .  The 
SN  size  was  varied  from  10  to  200  microns,  and  from  0  to  20%  concentration. 
Increasing  the  concentration  of  the  SN  particles  increased  the  burn  rate,  and 
decreasing  the  SN  particle  size  increased  the  burn  rate.  Both  of  these  effects  are  the 
same  trend  as  observed  with  the  AP  containing  propellant,  except  that  the  SN  is 
predicted  to  have  a  larger  effect  than  the  AP.  This  is  apparently  due  to  the  primary 
flame  temperature  of  SN  which  is  significantly  higher  than  that  of  AN  for  similar 
concentrations.  Although  there  are  no  experimental  data  available  with  which  to 
compare,  the  trends  would  seem  to  be  overpredicted.  Intuitively,  one  would  expect  AP 
to  have  a  greater  effect  than  SN. 


Figure  51.  Rum  rata  calculations  for  86%  solids  AN/SN/HTPB  propellants,  varyino  the 
concentration  and  size  of  the  SN. 

Propellants  Containing  AN/AP/SN 

A  primary  objective  of  the  program  was  to  model  AN  based  composite 
propellants  including  additional  oxidizers  of  AP  and  SN.  A  series  of  AN  propellants 
were  formulated  and  tested  as  part  of  the  experimental  part  of  the  programSI .  These 
propellants  contained  AP  and  SN  as  well  as  2%  AD  as  a  catalyst.  One  of  those 
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propellants,  BRM-270  has  been  compared  to  two  other  propellants  that  are 
comparable  with  respect  to  the  AN  fractions.  The  other  two  propellants  are  the  Talley 
propellant,  representing  a  propellant  without  additives,  and  the  BRM-1  propellant 
which  is  an  82%  solids  propellant  containing  AD  but  no  additional  oxidizer.  The  data 
and  the  model  calculations  are  compared  in  Figure  52.  As  would  be  expected,  the 
Talley  propellant  without  catalyst  and  with  a  relatively  low  amount  of  AN  has  the  lowest 
burn  rate.  Increasing  the  AN  content  and  adding  a  catalyst  increases  the  rate  from 
0.047  to  0.080  in/sec  (0.12  to  0.20  cm/sec).  Increasing  the  solids  content  to  86%  with 
10%  AP  and  10%  SN  increases  the  rate  to  0.135  in/sec  (0.34  cm/sec).  The  calculated 
rates  are  in  very  good  agreement  with  the  data  for  the  three  propellants.  However, 
there  is  some  deviation  at  the  higher  pressures.  The  BRM-1  propellant  without  any 
mixed  oxidizers  has  a  iower  exponent  than  predicted  while  the  mixed  oxidizer 
propellant  has  a  higher  exponent  than  predicted.  In  the  model  a  primary  flame 
reaction  order  of  1.5  was  assumed  which  will  give  pressure  exponents  on  the  order  of 
0.7.  Had  an  order  of  2.0  been  assumed  there  would  have  been  better  agreement  with 
the  mixed  oxidizer  propellant,  but  not  with  the  other  propellant.  Apparently  the  single 
global  reaction  used  in  the  model  is  not  adequate  to  model  the  wide  variation  of 
oxidizer/fuel  ratios  and  environments  of  the  different  propellants.  A  more 
comprehensive  kinetics  model  is  needed. 


Pressure  (psi) 


Figure  52.  Calculated  burning  rates  of  AN/HTPB  propellants  with  varying  amounts  of 
AN.  AP.  and  SN  compared  to  experimental  data. 

Temperature  sensitivity  calculations  were  also  made  for  the  same  propellants, 
although  there  is  not  any  temperature  sensitivity  data  for  comparison  purposes.  The 
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temperature  sensitivity  calculations  are  shown  in  Figure  53.  As  in  the  calculations 
above,  the  temperature  sensitivity  decreases  with  increasing  solids  loading  (with  the 
exception  of  the  78%  500  psi  point).  The  78%  solids  loading  propellant  has  the 
highest  temperature  sensitivity,  and  it  increases  slightly  with  increasing  pressure.  Both 
of  the  other  propellants  show  a  decreasing  temperature  sensitivity  with  increasing 
pressure,  and  the  86%  solids  propellant  has  the  lowest  temperature  sensitivity.  The 
calculated  temperature  sensitivity  values  are  all  higher  than  would  be  calculated  for 
similar  AP  propellants.  The  87.4%  AP  propellants  shown  in  Figure  35  had 
temperature  sensitivity  values  between  0.08  to  0.22  %/C.  In  Figure  53  the  86%  AN 
propeliant  has  temperature  sensitivity  values  from  0.29  to  0.20%/C.  It  appears  that  the 
temperature  sensitivity  of  AN  containing  propellants  will,  in  general,  be  higher  than  AP 
propeliants  due  to  the  inherently  lower  burn  rate  of  the  AN  propellants.  This  can  be 
seen  from  Equation  8.  The  flame  stand-off  distance,  5,  appears  in  the  denominator 
and  is  proportional  to  the  burning  rate  squared  (Equation  4),  which  leads  to  a  higher 
Op  for  lower  burning  rate  propellants. 


Figure  53.  Calculated  tamnerafure  sensitivity  for  AN/AP/SN  propellants,  varying  the 
oxidizer  concentrations  and  the  percentage  solids. 

The  six  mixed  oxidizer  propellants  from  the. experimental  paj^  of  the  program 
were  BRM-270,  271 , 272,  275,  277  and  278.  These  propellant  had  a  solids  content  of 
86%  with  2%  AD,  5  to  10%  AP  and  0  or  10%  SN.  Calculations  were  performed  to 
simulate  each  of  these  propellants,  and  a  comparison  of  calculated  and  experimental 
burning  rates  is  presented  in  Figure  54.  The  data  were  taken  at  four  pressures.  The 
agreement  between  the  model  and  the  data  is  reasonable,  although  the  model 
predicts  slightly  lower  rates  at  the  higher  pressures. 
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Experimental  Burn  Rate 


Figure  54,  A  comparison  of  calculated  and  experimental  burning  rates  for  the  six 
propellants  BRM-270.  271 , 272.  275.  277  and  278  at  varying  pressure. 

Formulations  corresponding  to  three  of  the  mixed  oxidizer  propellants  (BRM- 
270,  271  and  272)  were  used  to  study  the  effect  of  AN  particle  size  and  distribution. 
The  propellants  had  a  solids  content  of  86%  with  2%  AD,  10%  AP  and  10%  SN.  The 
results  are  shown  in  Figures  55  and  56  and  should  be  compared  to  the  results  shown 
in  Figures  44  and  45  for  a  propellant  without  the  mixed  oxidizers.  The  fine  size  was 
varied  from  10  to  100  microns  and  a  coarse  to  fine  ratio  of  70/30,  50/50  and  30/70  was 
used.  Increasing  the  concentration  of  fine  particles  caused  a  very  small  increase  in 
the  burn  rate,  and  changing  the  particle  size  had  a  similar  small  effect.  This  is 
essentially  the  same  as  noted  in  the  discussion  of  Figure  44  for  AN  propellants  without 
added  oxidizers.  It  seems  very  apparent  that  in  AN  propellants  the  primary  flame 
temperature  (energy  release)  is  so  low  that  the  kinetic  aspects  of  the  primary  flame 
dominate  over  the  diffusion  aspects,  virtually  eliminating  the  possibility  of  using  the  AN 
particle  size  to  tailor  the  propellant  burn  rate.  This  is  a  major  conclusion  from  the 
modeling  effort. 

For  the  smallest  particle  size  of  fines  the  rate  actually  decreased  slightly  with 
decreasing  size,  more  so  than  without  the  mixed  oxidizers.  As  noted  in  the  discussion 
of  Figure  44,  this  effect  is  felt  to  be  a  real  effect,  due  to  the  fuel  richness  of  the 
environment  surrounding  the  fine  particles.  This  effect  is  introduced  into  the  model 
through  Equation  (14). 
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Figure  55.  Burn  rate  caleulations  for  AN/AP/SN/HTPB  propellants,  varying  the 
noncentratinn  and  size  of  the  AN  fine  fraction  at  86%  solids.  Data  from  propellants 
BRM-270.  271 .  272  are  included  for  reference. 
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Figure  56.  Temperature  sensitivity  calculations  for  AN/AP/SN/HTPB  propellants^ 
varying  the  nnnnantration  and  size  of  the  AN  fine  fraction  at  86%- solids. 
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Temperature  sensitivity  calculations  performed  for  the  same  family  of 
propellants,  show  a  slight  effect  on  the  calculated  temperature  sensitivity,  but  almost 
identical  to  the  temperature  sensitivity  for  the  propellant  without  the  mixed  oxidizers. 
Decreasing  the  concentration  of  fine  particles  did  increase  the  temperature  sensitivity 
slightly  for  decreasing  fine  particle  size,  but  had  the  opposite  effect  for  larger  particles. 
The  value  of  temperature  sensitivity  remains  significantly  higher  than  corresponding 
AP  propellants. 


SUMMARY  AND  CONCLUSIONS 


Pfirformance  Calculations 

A  thermochemical  investigation  was  conducted  evaluating  the  performance  of 
AN  based  propellants.  AP/HTPB  propellants  have  a  maximum  Isp  fit  22  ^  aluminuni, 
while  AN/HTPB  propellants  have  an  equivalent  maximum  at  28%.  Unfortunately,  it  is 
difficult  to  get  good  combustion  efficiency  above  about  20%  aluminum,  which  reduces 
the  potential  practical  usage  of  AN  propellants. 

Performance  calculations  were  made  evaluating  lithium  nitrate  (LN),  sodium 
nitrate  (SN),  and  potassium  nitrate  (KN)  as  additives  to  scavenge  HCI  in  the  exhaust. 
LN  appears  to  be  the  best  HCI  scavenger,  and  has  a  higher  Isp  than  SN  or  KN. 
However,  it  is  more  hygroscopic  and  less  compatible  than  SN  or  KN,  and  SN  improves 
propellant  agglomeration  characteristics.  SN  also  has  a  higher  density  than  either  AN 
or  AP  which  makes  it  attractive  on  an  Isp-density  basis.  Therefore,  SN  ultimately 
became  the  additive  of  choice  for  the  study. 

Calculations  for  AN  propellants,  with  sufficient  SN  to  scavenge  HCI  to  less  than 
2%  of  the  exhaust  gases,  shows  that  increasing  AP  content  improves  the  Isp 
dramatically  due  to  its  inherently  higher  performance.  The  low  energy  of  AN  plus  its 
low  density  result  in  very  poor  performance  numbers  for  AN  containing  propellants. 

Combustion  Mechanisms 

Due  to  the  nitrate  in  AN,  it  was  initially  speculated  that  the  AN  flame  structure 
might  be  similar  to  double  base  propellants,  where  a  dual  stage  flame  exists. 
Therefore,  thermochemical  calculations  were  made  simulating  different  potential 
flames.  The  flame  temperature  calculated  assuming  NO  as  a  product  (analogous  to 
the  double  base  propellant  dark  zone)  is  less  than  600  K,  which  is  the  nominally 
measured  surface  temperature  of  burning  AN.  Therefore,  it  is  highly  unlikely  that  there 
is  an  inner  flame  with  NO  as  a  product  of  AN  combustion.  It  was  finally  concluded  that 
there  is  probably  a  single  AN  flame  which  is  in  equiiibrium  with  the  final  products  of 
N2,  H2O,  and  O2  find  at  a  monopropellant  flame  temperature  of  1247  K. 

A  critically  important  part  of  the  model  is  the  primary  diffusion  flame.  Initially  four 
separate  flame  structures  were  postulated  with  nitrogen  and  carbon  at  various  st^es 
of  oxidation.  Previous  work  on  AP  has  led  to  the  conclusion  that  the  formation  of  CO2 
is  relatively  slow  and  that  CO  is  the  product  of  importance  in  the  primary  flame.  It  was 
finally  concluded  that  the  AN  primary  flame  occurs  with  N2,  CO  and  H2O  as  the  final 
products.  Based  on  this  assumption,  AP/HTPB  propeliants  have  a  peak  temperature 
of  275.0  K  at  a  loading  cd  85%;  SN/HTPB  propellants  have  a  peak  temperature  slightly 
greater  than  2300  K  at  a  loading  of  84%;  while  the  peak  temperature  for  AN/HTPB 
propellants  is  barely  2000  K  and  occurs  at  a  loading  of  90%. 

A  composite  propellant  flame  structure  for  AP,  AN,  and  SN  has  been  postulated. 
AP  and  AN  each  burn  with  a  monopropellant  flame  whose  products  are  approximately 
30  and  14%  O2,  respectively.  This  leads  to  a  final  diffusion  flame  where  the  oxygen- 
rich  products  from  the  monopropellant  flame  react  with  the  binder  pyrolysis  products. 
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SN  on  the  other  hand,  does  not  burn  with  a  monopropellant  flame  but  only  burns  with 
a  single  diffusion  flame. 

Because  of  lower  flame  temperatures,  AN  propellants  have  a  much  lower  burn 
rate  than  AP  propellants.  For  example,  at  an  86%  loading  an  AP  propellant  has  a 
primary  flame  temperature  of  2700  K,  while  at  the  same  loading  an  AN  propellant  has 
a  flame  temperature  of  1700  K.  A  1000  K  decrease  in  flame  temperature  leads  to  an 
AN  propellant  burn  rate  five  times  lower  than  an  AP  propellant  rate.  This  is 
approximately  the  ratio  of  experimental  burn  rates  that  is  observed. 

Flame  temperatures  calculated  for  a  silicone  binder  were  significantly  higher 
than  for  a  corresponding  HTPB  binder  propellant  for  both  AP  and  AN.  Thus,  one 
would  expect  the  silicone  propellant  to  burn  better,  both  faster  and  more  efficiently, 
than  the  corresponding  HTPB  propellant  which  is  in  agreement  with  experimental 
observation. 

Experimental  data  for  AN  indicate  surface  temperatures  from  500  to  600  K,  with 
an  apparent  surface  activation  energy  of  ~10  Kcal/mole.  Data  for  SN  indicate  surface 
temperatures  of  ~1000  K,  and  an  activation  energy  of  ~4.5  Kcal/mole.  Apparently 
there  is  a  wide  diversity  of  surface  temperatures  and  individual  rates  occurring 
simultaneously  on  a  propellant  surface.  Because  of  its  higher  surface  temperature  for 
an  overall  propellant  rate,  an  ingredient  such  as  SN  will  be  slow  to  ignite  while  the 
thermal  profile  adjusts  to  a  configuration  compatible  with  the  SN  surface  kinetics. 
Once  ignited,  it  should  burn  very  rapidly,  possibly  leaving  an  empty  pocket  of  binder  or 
disrupting  adjacent  surface  regions  with  an  explosive  type  of  flow  due  to  its  high 
relative  rate.  AN  on  the  other  extreme  could  possible  puddle  and  flow  along  a  surface 
due  to  its  low  melting  point,  surface  temperature  and  combustion  rate,  and  long 
residence  time. 

AN  Monopropellant  Modeling 

The  literature  was  reviewed  for  pertinent  articles  on  AN  combustion  and 
modeling,  but  few  articles  were  found.  Because  of  the  extensive  previous  work  done 
with  AP,  it  was  decided  to  use  AP  as  a  basis  of  comparison  for  the  monopropellant 
model.  The  key  parameters  that  are  required  to  make  the  calculation  are  tabulated  in 
Table  4  (page  31)  with  a  comparison  of  the  values  used  for  both  AP  and  AN.  The  AP 
burn  rate,  temperature  sensitivity,  surface  temperature,  surface  heat  release,  and 
flame  standoff  distance  were  calculated  and  compared  to  experimental  data,  where 
available,  or  other  modeling  results. 

Using  the  known  physical  properties  for  AN,  the  following  combustion 
parameters  were  selected:  a  reference  burning  rate  of  0.19  cm/sec  at  68  atm,  a 
reference  surface  temperature  of  530  K,  reaction  order  of  2,  Es  of  14  Kcal/mole,  Ef  of 
25  Kcal/mole  and  a  Qc  value  of  0.26.  The  calculated  burn  rates  agree  very  well  with 

experimental  data,  and  a  Op  of  -0.3  %/K  was  calculated  (which  is  a  very  dependent  on 
the  selected  value  of  the  flame  activation  energy).  Although  AP  and  HMX  have  Op 
values  on  the  order  of  0.16  and  0.2%/K  respectively,  it  seems  logical  that  AN  would 
have  a  relatively  high  Op  due  to  its  lower  surface  temperature  and  low  burning  rate 
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(similar  to  double  base  propellants  which  have  Op  values  of  0.2  to  0.5  %/K).  The 
measured  and  calculated  surface  temperature  varies  between  500  and  600  K  for  the 
AN  monopropellant.  This  puts  AN  in  the  same  range  of  surface  temperatures  as 
double  base  propellants,  but  significantly  below  AP.  Variations  in  Es  indicate  that  a 
value  between  ~10  and  15  Kcal/mole  would  fit  the  surface  temperature  data  equally 
well.  The  molten  AN  layer  was  calculated  to  be  15  to  50  microns  over  the  pressure 
range  of  interest,  which  compares  qualitatively  with  window  bomb  movies  of  burning 
AN  propellants.  The  calculated  flame  stand-off  distances  of  25  to  250  microns  are 
much  larger  than  those  values  calculated  for  AP  or  HMX,  but  are  comparable  to  those 
observed  for  double  base  propellants. 

AN  Composite  Propellant  Model 

The  basis  of  the  AN  model  is  the  Beckstead  SST  model.  The  SST  model  is 
based  on  a  separate  energy  balance  for  each  ingredient  type  and  size,  and  the  overall 
burning  rate  is  calculated  on  a  time  averaged  basis  rather  than  the  original  BDP 
approach  of  using  a  space  averaged  approach.  Other  aspects  of  the  SST  model  that 
influence  AN  propellants  are  summarized  briefly. 

To  model  composite  propellant  combustion  one  must  assume  a  distribution  of 
fuel  among  the  different  oxidizer  fractions.  In  the  SST  model  it  has  been  assumed  that 
the  binder  is  distributed  according  to  the  oxidizer  specific  surface  area.  However, 
binder  thicknesses  for  intermediate  and  coarse  oxidizer  fractions  are  ~1%  of  the 
particle  diameter  which  is  smaller  than  the  roughness  and  non-sphericity  of  actual 
crystals.  Thus,  a  significant  amount  of  binder  can  be  buried  within  the  surface 
irregularities  of  a  crystal  in  excess  of  what  would  be  calculated  for  a  purely  spherical 
particle.  To  compensate  for  this,  Cohen  assumes  that  the  finest  oxidizer  fraction  burns 
at  the  stoichiometric  0/F  condition,  and  the  fuel  left  over  from  one  particle  size  is 
available  to  be  burned  in  the  diffusion  flame  of  the  next  larger  size  fraction.  This 
approach  appears  to  be  quite  realistic  on  a  mechanistic  basis,  and  has  been  adopted. 

Multimodal  and  mixed  oxidizer  propellants  experience  interaction  effects  and 
influences  caused  by  different  sizes  and  types  of  particles  having  significantly  different 
burn  rates.  Either  faster  burning  components  are  igniting  slower  burning  particles  all 
along  their  periphery,  or  slower  burning  particles  can  be  undercut  by  faster  burning 
components.  The  model  attempts  to  account  for  this  effect  also. 

A  critical  aspect  of  a  composite  propellant  model  is  the  description  of  the 
monopropellant  flame  and  primary  diffusion  flame  interaction.  Within  the  BDP  type  of 
approaches,  the  Pf  term  accounts  for  this  interaction,  and  is  very  dominant  within  the 
model  in  determining  the  effect  of  pressure  and  particle  size  on  burn  rate.  This  effect  is 
still  accounted  for  in  a  relatively  crude  manner,  and  a  more  rigorous  calculational 
procedure  needs  to  be  developed. 

Once  the  amount  of  oxidizer  that  participates  in  the  primary  flame  is  determined, 
then  the  amount  of  the  heat  that  is  transferred  to  the  oxidizer  versus  the  amount 
transferred  to  the  binder  must  be  determined.  This  effect  is  accounted  for  by  the  term 
pP  in  the  model.  Evaluation  of  this  effect  is  complicated  by  the  fact  that  the  primary 
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flame  flips  back  and  forth  from  being  over  the  oxidizer  to  being  over  the  binder,  during 
the  lifetime  of  a  particle.  This  general  effect  has  been  addressed  previously,  but  a 
quantitative  description  of  the  transient  heat  transfer  was  beyond  the  scope  of  the 

current  project,  and  an  integrated,  average  value  of  Pp  is  estimated  empirically. 

AP  Propellant  Results 

'  Miller's  non-metallized  AP/HTPB  propellants  were  used  as  a  comparative  data 
base,  because  of  the  large  number  of  propellants  with  the  same  binder  and  a  wide 
range  of  compositions,  pressure,  particle  size,  and  initial  temperatures.  Calculations 
were  compared  to  the  experimental  data,  and  model  parameters  were  optimized  to 
provide  the  best  fit  of  the  model  to  the  data.  Model  calculations  for  Miller’s  Series  III 
data  at  1000  psi  (6.8  MPa)  are  very  reasonable  with  only  a  few  of  the  widest 
distribution  propellants  falling  outside  of  a  10%  predictive  goal.  The  calculated 
pressure  exponents  for  most  of  data  were  within  20%  of  the  measured  values,  but  the 
calculations  were  consistently  low  with  several  datum  points  for  propellants  containing 

400|i  AP  more  than  20%  in  error.  Calculated  Op  values  were  very  reasonable  with 
only  three  calculations  lying  outside  of  the  20%  range  of  the  data  and  two  of  those 
calculations  were  for  propellants  containing  400ii  AP. 

AN  Propellant  Results 

Burn  rate  calculations  were  made  for  AN/HTPB  propellants  (with  no  additives) 
varying  AN  content  from  76%  to  88%.  The  calculated  rates  vary  almost  linearly  from 
0.05  to  0.1  in/sec  (0.1  to  0.2  cm/sec)  which  seems  very  reasonable.  Pressure 
exponents  of  0.7  to  1 .0  were  calculated,  which  appear  to  be  a  little  high. 

Increasing  the  primary  flame  activation  energy  causes  burn  rate  to  increase  as 
AN  content  (and  the  corresponding  primary  flame  temperature)  increases.  An 
activation  energy  of  15  seemed  to  give  rates  higher  than  would  be  expected,  while  a 
value  of  10  seemed  to  give  more  reasonable  results.  For  AP  and  HMX/HTPB 
propellants  a  value  of  15  Kcal/mole  is  normally  used  as  the  primary  flame  activation 
energy. 

Burn  rate  calculations  for  propellants  containing  ammonium  dichromate  as  a 
catalyst  show  the  rate  increases  with  increasing  catalyst  concentration  and  with 
increasing  percentage  solids,  which  is  in  general  agreement  with  experimental  data. 
The  calculated  effect  of  the  catalyst  is  small  due  to  the  inherently  low  energy  of  AN. 
The  low  flame  temperature  (i.e.  the  low  energy)  of  AN  propellants  prohibits  a  catalysis 
comparable  to  AP  propellants.  The  catalysis  process  has  a  lower  pressure 
dependence  than  the  uncatalyzed  process,  and  the  net  effect  is  to  increase  the  burn 
rate  slightly,  but  reduce  the  pressure  exponent. 

Virtually  no  particle  size  effect  was  calculated  for  a  78%  AN  propellant,  but  there 
was  a  slight  effect  calculated  for  an  86%  solids  loading.  For  the  lowest  concentration 
of  fines  the  rate  actually  decreased  slightly  with  decreasing  size.  The  AN  environment 
is  very  fuel  rich  and  the  decrease  in  rate  with  decreasing  particle  size  corresponds  to 


an  Gxcessiv©  ©nergy  loss  to  th©  bind©r.  Th©  corr©sponding  ©ff©ct  has  b©©n  seen  with 
very  wide  distribution  AP  propellants  and  in  self  extinguishing  AP  propellants. 

A  critical  aspect  in  the  development  of  the  model  is  the  ability  to  handle  mixed 
oxidizers,  specifically  AP.  SN,  and  AN.  Increasing  AP  concentration  increased  burn 
rate,  and  decreasing  AP  particle  size  increased  burn  rate.  Both  of  these  effects  are  in 
agreement  with  experimental  observation.  Increasing  fine  AP  concentration  increased 
the  exponent  slightly,  but  in  general,  the  exponent  did  not  change  significantly  with 
increasing  AP  concentration.  The  calculations  were  compared  with  propellants 
containing  5%  5  micron  AP,  10%  5  micron  AP;  and  10%  200  micron  AP.  The  5  micron 
propellant  calculations  are  in  very  reasonable  agreement  with  the  data,  but  the  200 
micron  calculations  are  lower  than  the  corresponding  data. 

Increasing  the  concentration  of  SN  increased  the  burn  rate,  and  decreasing  the 
SN  particle  size  increased  the  burn  rate.  Both  of  these  effects  are  the  same  trend  as 
obsen/ed  with  AP,  except  that  the  SN  is  predicted  to  have  a  larger  effect  than  the  AP. 
Although  there  are  no  experimental  data  available  with  which  to  compare,  the  trends 
would  seem  to  be  overpredicted.  Intuitively,  one  would  expect  AP  to  have  a  greater 
effect  than  SN. 

Calculations  were  performed  to  simulate  six  propellants  containing  86%  soiids 
with  2%  AD,  5  to  10%  AP  and  0  or  10%  SN.  The  agreement  between  the  model  and 
the  data  is  reasonable,  although  the  model  predicts  slightly  lower  rates  at  the  higher 
pressures.  Increasing  the  concentration  of  fine  AN  particles  caused  a  very  small 
increase  in  the  burn  rate,  and  changing  the  AN  particle  size  had  a  similar  small  effect. 
In  AN  propellants  the  primary  flame  temperature  is  so  low  that  the  kinetic  aspects  of 
the  primary  flame  dominate  over  the  diffusion  aspects,  virtually  eliminating  the 
possibility  of  using  AN  particle  size  to  tailor  propellant  burn  rate.  This  is  a  major 
conclusion  from  the  modeling  effort. 

Calculated  temperature  sensitivity  decreases  significantly  with  increasing  AN 
content  (i.e.  solids  loading)  and  with  decreasing  activation  energy.  With  no  precise 
data  to  compare  with  the  values  predicted  with  an  activation  energy  of  10  Kcal/mole 

appeared  to  give  reasonable  Op  values,  0.2  to  0.3  %/K,  particularly  at  the  lower 
concentration  levels.  The  Op  values  are  quite  high  compared  to  conventional  AP 
propellants,  but  limited  data  for  AN  propellants  indicate  tck  values  of  0.2  to  0.3%/F, 
which  are  about  double  those  normally  expected  for  AP  propellants  and  in  very 
reasonable  agreement  with  the  calculated  ap  values.  The  addition  of  AP  and  SN  to 
the  propellants  had  little  effect  on  the  calculated  Op  values.  For  the  mixed  oxidizer 
propellants,  decreasing  the  concentration  of  fine  AN  particles  did  increase  the 
temperature  sensitivity  slightly  for  decreasing  fine  AN  particle  size,  but  had  the 
opposite  effect  for  larger  AN  particles.  The  calculated  Op  value  remains  significantly 
higher  than  corresponding  AP  propellants. 
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Modeling  Needs 


The  following  is  an  abbreviated  list  of  areas  where  additional  work  is  needed  to 
further  advance  the  understanding  of  composite  propellant  combustion. 

•  Phenomenological  needs  in  composite  propellant  modeling: 

•  a  better  understanding  of  binder  pyrolysis 

•  a  better  description  of  how  the  fuel  is  distributed  relative  to  oxidizer 

particles  in  the  combustion  process 

•  an  accurate  description  of  the  interaction  of  the  diffusion  flame  and  the 

monopropellant  flame  (i.e.  (3f) 

•  a  better  understanding  of  the  dynamic  diffusion  flame  shape  during  the 

life  time  of  a  burning  particle 

•  a  description  of  the  diffusion  flame  heat  transfer  both  to  the  oxidizer  as 

well  as  the  binder  (i.e.  pp)  during  the  life  time  of  a  burning  particle 

•  a  description  of  the  interaction  between  flames  from  adjacent  oxidizer 

particles 

•  the  effect  of  aluminum  agglomeration  and  combustion  on  all  of  the  above 

•  determination  of  the  number  of  kinetic  steps  necessary  to  describe  the 

essential  combustion  features  of  typical  flames  including  the 
condensed  phase  reactions 

•  Experimental  needs  to  support  modeling  of  composite  propellants: 

•  accurate  surface  temperature  measurements  of  any  ingredient  intended 

for  use  in  composite  propellants  (i.e.  AP,  AN,  SN,  HMX,  HTPB,  etc) 

•  surface  decomposition  products  of  composite  propellant  ingredients  (i.e. 

AP,  AN,  SN,  HMX,  HTPB,  etc)  closer  to  actual  burn  rate  conditions 

•  surface  heat  release  and  heat  flux  values 

•  Op  data  of  oxidizers  intended  for  use  in  composite  propellants  (i.e.  AP, 

AN,  SN,  HMX,  etc) 

•  systematic  combustion  data  (i.e.  r,  n,  and  Op)  varying  the  oxidizer  content 

and  particle  size  to  allow  identification  of  global  kinetic  parameters 

•  systematic  combustion  data  (i.e.  r,  n,  and  ap)  varying  the  oxidizer  content 

and  particle  size  to  allow  identification  of  global  kinetic  parameters 

•  species  composition  data  through  flames  to  identify  which  gas  phase 

species  are  important  in  the  flames 

•  identification  of  key  species  needed  to  describe  the  principal  kinetic 

steps  in  both  the  condensed  and  gas  phases 
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NOMENCLATURE 


A 


a 

AP 

AN 

b 

BDP 

Cign 

CMDB 

Cp 

D' 

Do 

E 


frb 

Hex 

HMX 

HTPB 

k 

KN 

LN 

m 


mi 

MW 

n 

P 

Q 

Qc 

Qfuel 

Qf 
Ql 
Qlm 
ror  rp 
R 

RDX 
So 
SN 
SST 
.  t 
tb 
^gn 

Tf 

To 

V 

X* 

XLDB 


a 


Arrhenius  frequency  factor 

numerical  coefficients  in  flame  standoff  expression 

Ammonium  perchlorate 

Ammonium  nitrate 

Characteristic  surface  dimension 

Beckstead-Derr-Price  model 

Coefficient  in  ignition  delay  time  calculation 

Composite-modified  double-base 

Average  mean  heat  capacity  for  the  solid  and  gases 

Statistical  intersection  diameter 

Particle  diameter 

Activation  energy 

Fraction  of  binder  reacted 

Fraction  of  oxidizer  reacted 

Heat  of  explosion 

Cyclotetramethylene  tetranitramine 
Hydroxyterminated  polybutadiene 
Kinetic  rate  constant 
Potassium  nitrate 
Lithium  nitrate 

Mass  flux  associated  with  propellant  components 

Total  mass  flux  of  propellant 

Molecular  weight 

Pressure  exponent  of  burn  rate 

Pressure 

Heat  release  associated  with  combustion  steps 
Proportionality  constant  in  the  equation  for  Ql 
Heat  of  pyrolysis  of  the  fuel  binder 
Heat  release  in  monopropellant  flame 
Heat  of  gasification  of  the  oxidizer 
Heat  of  fusion  of  the  metal 
Linear  burning  rate 
Gas  constant 

Cyclotrimethylene  trinitramine 
Total  surface  area 
Sodium  nitrate 

Separate  surface  temperature 
Characteristic  time 
Binder  burnthrough  time 
Ignition  delay  time 
Temperature 

Adiabatic  flame  temperature  of  the  propellant 

Initial  temperature  of  the  propellant 

Gas  velocity 

Flarhe  standoff  distance 

Crosslinked  double-base 

Weight  fraction  oxidizer 
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at  Thermal  diffusivity 

Pap  Fraction  of  oxidizer  reactants  available  for  final  flame 

Pf  Fraction  of  oxidizer  reactants  involved  in  the  primary 

diffusion  flame 

Pp  Fraction  of  heat  transferred  from  primary  flame  to  oxidizer 

Yst  Stoichiometric  oxidizer  to  fuel  ratio  in  primary  diffusion  flame 

ATioss  Temperature  loss  differential  from  oxidizer  to  binder 

5i  Distance  from  oxidizer  to  center  of  adjacent  binder 

5  Reaction  order 

<!)ST  Stoichiometric  oxidizer  to  fuel  ratio 

Nondimensionai  flame  standoff  distance 
Op  Temperature  sensitivity  of  burn  rate 

p  Density 

X  Thermal  conductivity  of  gases 


Subscripts 

Al 

Aluminum 

b 

Fuel  binder 

Comb 

Combustion 

DB 

Double  base 

f 

Monopropellant  flame  conditions 

fuel 

Fuel  binder 

FF 

Final  flame  conditions 

g 

Gas 

i 

Particle  size  counter 

igh 

Ignition 

j 

Oxidizer  type  counter 

MONO 

Oxidizer  monopropellant  flame  conditions 

0 

Initial  conditions 

ox 

Oxidizer 

P 

Solid  propellant 

PD 

Diffusion  part  of  primary  flame 

PF 

Primary  flame 

s 

Propellant  surface  conditions 
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APPENDIX  A 


Test  Case 

#1 

#2 

#3 

#4 


Computer  Test  Cases 

Description 

AN  Monopropeilant  test  case  calculating  temperature  sensitivity. 
(Note:  A  binder  must  be  included  for  program  to  run.  Use  small 
weight  %  ~  0.01 ) 

AN  monopropeilant  data  compared  statistically  with  statistics  input, 
AN/AP/SN/HTPB  mixed  oxidizer  calculation  with  standard  output. 
AN/HTPB  test  case  with  simplified  output. 
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0,  0,  -0,  1 

’AN  MONOPROPELLANT  —  TEST  CASE’ 

,■  ,’AN’,  ,’HTPB’ 

0.0, 0.0, 0.0,  0.0, 0.0, 0.0,  1.0, 0.0, 0.0,  0.0, 0.0, 0.0 
1.0, 1.0, 1.0,  1,0, 1.0, 1.0,  999. ,1.0, 1.0,  1.0, 1.0, 1.0 
300.0,  0.0,  ,  ,  ,  0.0,  ,  ,  ,  0.0,  0.0,  0.0 
SNAMl  PSTOP=280.,  SEND 

SNAMl  NJOB  =  3  SEND 

000000000111111111122222222223333333333444444444455555555556666666666777 

123456789012345678901234567890123456789012345678901234567890123456789012 


13-MAR-89  AN  MONOPROPELLANT  —  TEST  CASE 

OXIDIZER  DEPENDENT  VARIABLES: 

TYPE - WT"« . DIA'(WT'CRT)FIS1 — 

AN  1  99.9  999.00 

BDSQ  ARATIO  ERATIO  D0XNUM 
1.0002  9999.00  9999.00  1.22 

OXIDIZER  CONSTANTS: 

"TYPE  iJ'C'OEF  PHOX  AOX  EOX  HINAP  EAP  TMONO  AAP  KAPl 

AN  0.290  1.72  4.51E+05  12000.  2.0  25000.  1191.  127.  0.328E-02 

CIGN  TMELT  POWD  TLOSS  FRD  GEOM  HIEAP  TMONOH  HIAAP 

0.00  443.0  2,00  0.00  0.0  0.000  25000.  1191.  182, 


DELTA  SFE  DEL 
0.00  1.000  0.0000 


MISC  PARAMETERS: 

EEDBP  XL'A'MB  XFH  TZERO  BPI  RHOP  BETMIN  BPMIN  BFDIV 

0.30  0.0003  3.00  219.  1.50  1.72  0.00  0.10  999999.00 

PRESSURE  DEPENDENT  RESULTS: 


OX 

“PTATWT7FSI - 

POX  TSOX  X* 

■PD  BETAP  BETAF  X*PF  X*AP 

QL  ZAP  ZPF 

DNOMl  DN0M2 

DN0M3 

DN0M4 

CNVRG 

AN 

34.0/ 

500. 

1 

0.029  492.  0. 

000  0.001  0. 

002  0. 

29  335 

.7 

79.2  4,28  0.15 

13.485  0.000 

0 . 000 

0 . 000 

2 

AN 

68.0/ 

1000. 

1 

0.058  524.  0. 

000  0.001  0. 

002  0. 

17  166 

.4 

- 

88.3  4.23  0.15 

6.799  0.000 

0.000 

0.000 

2 

AN 

136.1/ 

2000. 

1 

0.114  559.  0. 

000  0,001  0. 

002  0. 

10  81 

.6 

- 

98.6  4.07  0.15 

3.466  0.000 

0 . 000 

0 . 000 

2 

BINDER 

RB 

TSB 

FRB 

DELTB 

ZP  EXZP 

QFUEL 

HTPB 

0.0296 

850.0 

0.99 

0.0  0 

.02  10. 

433.0 

7 

HTPB 

0.0288 

850.0 

0.99 

0.0  0 

.01  10. 

433.0 

6 

HTPB 

0.0290 

850.0 

0.99 

0.0  0 

.00  10. 

433.0 

6 

BURN  RATE  SUMMARY 

TTA 

RBAR 

N 

IT  RB  ROX 

500. 

0,0292 

0.00 

5  0.0296  0.0292 

1000. 

0.0579 

0.98 

6  0.0288  0,0579 

2000. 

0.1136 

0.97 

6  0.0290  0.1136 

0.370  CPU  SECONDS  USED 
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13-MAR-89 


AN  MONOPROPELLANT 


TEST  CASE 


OXIDIZER  DEPENDENT  VARIABLES: 
TYPE  Wt  i,  DIA  (MICRONS) 

AN  1  99.9  999.00 


BDSQ  ARATIO  ERATIO  D0XNUM 
1.0002  9999.00  9999.00  1.22 


DELTA  SFE 
0 . 00  1 . 000 


DEL 

0.0000 


OXIDIZER  CONSTANTS: 

“TYPE  QCOEP  RHOX  AOX  EOX  HINAP  EAP  TMONO  AAP  KAPl 

AN  0.290  1.72  4.51E+05  12000.  2.0  25000.  1247.  127.  0.528E-02 

CIGN  TMELT  POWD  TLOSS  FRD  GEOM  HIEAP  TMONOH  HIAAP 

0.00  443.0  2.00  0.00  0,0  0.000  25000.  1247.  182. 

MISC  PARAMETERS: 

rEUEP  xLamb  XFH  TZERO  BPI  RHOP  BETMIN  BPMIN  BFDIV 

0.30  0.0003  3.00  298.  1.50  1.72  0.00  0.10  999999  00 

PRESSURE  DEPENDENT  RESULTS: 

P(AiM)/PSI  PUX  TSOX  X*PD  BETAP  BETAF  X*PF  X*AP 

QL  ZAP  ZPF  DNOMl  DN0M2  DN0M3  DNOM4  CNVRG 


34.0/  500. 

68.0/  1000. 
136.1/  2000. 


0.037  496.  0.000  J 

-57.4  4.35  0,15  10.503 

0.074  529.  0.000  J 

-66,9  4,31  0.15  5.290 

0.148  567.  0.000  « 

-77.9  4.26  0.15  2.665 


BINDER  RB  TSB  FRB 

HTPB  0.0307  850.0  0.99 

HTPB  0.0309  850,0  0.99 

HTPB  0.0312  850.0  0.99 

BURN  RATE  SUMMARY 

“PTS  PBXR  N  IT  RB  ROX 

500.  0.0375  0.00  6  0.0307  0.0375 

1000.  0.0744  0.99  6  0.0309  0.0744 

2000,  0.1477  0.99  6  0.0312  0.1477 


001 

0. 

002 

0 

.27 

267 

.8 

0. 

000 

0 

.000 

0 

.  000 

001 

0, 

002 

0, 

.16 

133, 

.0 

0. 

000 

0. 

.  000 

0 

.000 

001 

0. 

001 

0. 

09 

66. 

.0 

0. 

000 

0. 

,000 

0. 

,  000 

DELTB 

ZP 

EXZP 

QFUEL 

0 

.0 

0 

.02 

10. 

433 

,0 

6 

0 

.0 

0 

.00 

10. 

433 

.0 

6 

0 

.0 

0 

.00 

10. 

433 

.0 

7 

0.0307  0.0375 

0.0309  0.0744 


0.180  CPU  SECONDS  USED 
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13-MAR-89  AN  MONOPROPELLANT  —  TEST  CASE 

OXIDIZER  DEPENDENT  VARIABLES: 

TYFE  Wr . }4  DiA  (MICRONS) — 

AN  1  99.9  999.00 


BDSQ  ARATIO  ERATIO  D0XNUM  DELTA  SFE  DEL 

1.0002  9999.00  9999.00  1.22  0.00  1.000  0.0000 

OXIDIZER  CONSTANTS; 

"T7PH  (JCOEF  RHOX  AOX  EOX  HINAP  EAP  TMONO  AAP  KAPl 

AN  0.290  1.72  4.51E^05  12000.  2.0  25000.  1282.  127  0  695E-02 

hieap  'tmonoS  h!aap 

0.00  443.0  2.00  0.00  0.0  0.000  25000.  1282.  182. 

MISC  PARAMETERS: 

CiUBK  XLAMB  AFH  TZERO  BPI  RHOP  BETMIN  BPMIN  BFDIV 

0.30  .  0.0003  3.00  347.  1.50  1.72  0.00  0  10  999999  00 

PRESSURE  DEPENDENT  RESULTS:  999999.00 

^(ArM)/PSI  RTJX  TSOX  X*PD  BETAP  BETAF  X*PF  X*AP 

AM  OA  a/  cr,n,  DNOMl  DN0M2  DN0M3  DN0M4  CNVRG 


PCAiM)/PSI 

34.0/  500. 


68.0/  1000. 


136.1/  2000 


0.043  496.  0.000  0.001  0.002  0.26  235.4 

-43.3  4.39  0.15  9.093  0.000  0.000  0.000  2 

0.086  530.  0.000  0.001  0.002  0.15  116.8 

-53.2  4.35  0.15  4.581  0.000  0.000  0.000  2 

0.171  569.  0.000  0.001  0.001  0.09  57.9 

-64.5  4.30  0.15  2.309  0.000  0.000  0.000  "  2 


BINDER  RB  TSB  FRB 

HTPB  0.0317  850.0  0.99 

HTPB  0.0324  850.0  0.99 

HTPB  0.0327  850.0  0.99 

BURN  RATE  SUMMARY 

PTA  KBAR  N  it  RB  ROX 

500.  0.0433  0.00  7  0.0317  0.0433 

1000.  0.0859  0.99  6  0.0324  0.0859 

2000.  0.1705  0.99  6  0.0327  0.1705 


DELTB 

0.0 

ZP 

0.01 

EXZP 

10. 

QFUEL 

433.0 

5 

0,0 

0 . 00 

10. 

433,0 

6 

0,0 

0 . 00 

10. 

433.0 

7 

0.180  CPU  SECONDS  USED 
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.y 


TEMPERATURE 

SENSITIVITY 

DATA  BETWEEN 

-65.  AND  165. 

DEG  F 

PSI 

500, 

1000. 

2000. 

RATE  (IN/SEC) 
0.038 
0.074 
0.148 

EXPONENT 

0 . 00 
0.99 
0.99 

SIG  P  (1/C) 
0.00307 
0.00309 
0.00317 

PI  K  (1/F) 
1.00000 
0.14850 
0.15789 

AVE  TS(K) 

496.6 

530.7 
569.6 

COMBUSTION  INSTABILITY  CHARACTERISTICS  USING  THE  DENISON  AND  BAUM  MODEL 

TDIFF  =  0.000300 


P=1000.  A; 

=  3.9  B= 

1.40 

P=2000. 

A=  4.1  B: 

=1.17 

P=  0. 

A=  0.0  8 

rr0 , 00 

FREQ 

OMEGA 

REAL  PC 

IMAG 

OMEGA 

REAL  PC 

IMAG 

OMEGA 

REAL  PC 

IMAG 

100. 

34.04 

0.67 

-0.44 

8.64 

1.26 

-0.49 

0 . 00 

0 . 00 

0 . 00 

200. 

68.08 

0.48 

-0.36 

17.29 

0.91 

-0.55 

0 . 00 

0 . 00 

0  •  00 

300. 

102.12 

0.39 

-0.31 

25.93 

0.73 

-0.51 

0.00 

0.00 

0.00 

400 , 

136.16 

0.33 

-0.28 

34.57 

0.62 

-0.47 

0.00 

0 . 00 

0 . 00 

500 . 

170.20 

0.30 

-0.26 

43.21 

0.55 

-0.44 

0.00 

0 . 00 

0 . 00 

600 . 

204.24 

0.27 

-0.24 

51.86 

0.50 

-0.41 

0 . 00 

0 . 00 

0 , 00 

700 . 

238.28 

0.25 

-0.22 

60.50 

0.46 

-0.39 

0 . 00 

0.00 

0 . 00 

800 . 

272.32 

0.24 

-0.21 

69,14 

0.43 

-0.37 

0.00 

0 . 00 

0 . 00 

900 , 

306.37 

0.22 

-0.20 

77.79 

0.40 

-0.35 

0.00 

0.00 

0 . 00 

1000. 

340.41 

0.21 

-0.19 

86.43 

0.38 

-0.33 

0.00 

0 . 00 

0 . 00 

2000 . 

680.81 

0.15 

-0.14 

172.86 

0.26 

-0.25 

0.00 

0 . 00 

0 . 00 

3000 . 

1021.22 

0.12 

-0.11 

259.29 

0.21 

-0.20 

0 . 00 

0 . 00 

0 . 00 

4000 , 

1361.62 

0.10 

-0.10 

345.72 

0.18 

-0.18 

0.00 

0 . 00 

0  •  00 

5000 . 

1702.03 

0.09 

-0.09 

432.14 

0.16 

-0.16 

0 . 00 

0.00 

0 . 00 

6000 . 

2042,44 

0.09 

-0.08 

518.57 

0,15 

-0.15 

0.00 

0.00 

0 . 00 

7000 . 

2382.84 

0.08 

-0.08 

605.00 

0.14 

-0.14 

0.00 

0 . 00 

0 . 00 

8000 , 

2723.25 

0.07 

-0.07 

691.43 

0.13 

-0.13 

0 . 00 

0 , 00 

0  •  00 

9000 , 

3063.66 

0.07 

-0.07 

777.86 

0.12 

-0.12 

0.00 

0.00  . 

0.00 

IN* 
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13-MAR-89 


AN  MONOPROPELLANT 


TEST  CASE 


OXIDIZER  DEPENDENT  VARIABLES: 

TYPE - WT  «  DIA'CMICROT^ 

AN  1  99.9  999.00 

BDSQ  ARATIO  ERATIO  D0XNUM 
1.0002  9999.00  9999.00  1.22 


DELTA  SEE  DEL 
0,00  1.000  0.0000 


OXIDIZER  CONSTANTS: 

“TYPE - 5EUEF - PROX  AOX 

AN  0.290  1.72  4.51E+05 

CIGN  TMELT  POWD 
0,00  443.0  2.00 


EOX  HINAP  EAP  "TMONO  AA'P"  KAPl 
12000.  2.0  25000.  1191.  127.  0.32SE-02 

TLOSS  FRD  GEOM  HIEAP  TMONOH  HIAAP 

0  00  0.0  0.000  25000.  1191,  182, 


CS'OBP — YEAmP  AFH  tzero  bpi  rhop  betmin 

0.30  0.0003  3.00  219.  1.50  1.72  0.00 

PRESSURE  DEPENDENT  RESULTS: 

“DY 

P(ATM)/FSi  ROX  TSOX 

QL  ZAP  ZPF  DNOMl 

AN 

34,0/ 

500. 

1 

0.029  492. 

-79.2  4,28  0. 

0.000 

15  13.485 

AN 

68.0/ 

1000. 

0 . 000 

15  6.799 

1 

0.058  524. 

-88.3  4.23  0. 

AN 

136.1/ 

2000. 

0.000 

15  3.466 

1 

0.114  559. 

-98.6  4.07  0. 

AN 

204.1/ 

3000. 

0 . 000 

16  2.262 

1 

0.174  584. 

-105.8  3.96  0. 

AN 

272.1/ 

4000. 

. 

1 

0,239  604. 

-111.6  3.87  0. 

0 . 000 

15  1.648 

AN 

340,1/ 

5000. 

1 

0.303  620. 

-116.3  3.81  0, 

0.000 

15  1.300 

BINDER 

RB 

TSB  FRB 

HTPB 

0.0286 

850,0  0.99 

HTPB 

0.0288 

850.0  0.99 

HTPB 

0.0290 

850.0  0.99 

HTPB 

0.0310 

850.0  0.99 

HTPB 

0.0312 

850.0  0.99 

HTPB 

0.0366 

850.0  0.99 

BURN  RATE  SUMMARY 

ROX 

PTA 

- RBAR” 

N  IT  RB 

500. 

0,0292 

0.00  6  0.0286  0 

.0292 

1000. 

0.0579 

0.98  6  0.0288  0 

.0579 

2000. 

0,1136 

1.00  6  0.0290  0 

.1136 

3000. 

0.1738 

1.07  5  0.0310  0 

.1741 

4000 . 

0.2387 

1.08  5  0.0312  0 

,2389 

5000. 

0.3012 

1.04  4  0.0366  0 

.3029 

BPMIN  BFDIV 

0.10  999999.00 

BETAP  BETAF  X*PF  X*AP 

DN0M2  DN0M3  DN0M4  CNVRG 

0.001  0.002  0.29  335.7 

0.000  0.000  0.000  2 

0.001  0.002  0.17  166.4 

0 . 000  0 . 000  0 . 000  2 

0.001  0.002  0.10  81.6 

0.000  0.000  0.000  2 

0.001  0.002  0.07  51.9 

0.000  0.000  0.000  2 

0.001  0.001  0.06  37.0 

0.000  0,000  0,000  2 


0.001  0.001  0.05  28.7 

0.000  0.000  0.000  2 


DELTB 

ZP 

EXZP 

QFUEL 

0,0 

0,02 

10, 

433.0 

6 

0,0 

0.01 

10. 

433.0 

6 

0.0 

0.00 

10. 

433.0 

6 

0.0 

0.00 

10. 

433.0 

8 

0.0 

0.00 

10. 

433.0 

8 

0.0 

0 .00 

10. 

433.0 

9 

0.320  CPU  SECONDS  USED 


92 


13-MAR-89 


AN  MONOPROPELLANT 


TEST  CASE 


OXIDIZER  DEPENDENT  VARIABLES: 


99.9 


999.00 


BDSQ  ARATIO  ERATIO  D0XNUM 
1.0002  9999.00  9999.00  1.22 


DELTA  SFE  DEL 

0.00  1.000  0,0000 


OXIDIZER  CONSTANTS; 

■  TYPE  CfnJEF  l?HOX  AOX  EOX  HINAP  EAP  -TMONO-  AAP  -  KAPl 

AN  0.290  1.72  4.51E+05  12000.  2.0  25000.  1247.  127.  0.528E-02 

CIGN  TMELT  POWD  TLOSS  FRD  GEOM  HIEAP  TMONOH  HIAAP 

0.00  443,0  2.00  0.00  0.0  0.000  25000.  1247.  182. 

MISC  PARAMETERS: 

CsUbp  xlamb  ?:fh  TZERO  BPI  RHOP  BETMIN  BPMIN  BFDIV 

0.30  0.0003  3.00  298.  1.50  1.72  0.00  0.10  999999.00 

PRESSURE  DEPENDENT  RESULTS: 

PCArM)/PSI  f?UX  TSOX  X*PD  BETAP  BETAF  X*PF  X*AP 

QL  ZAP  ZPF  DNOMl  DN0M2  DN0M3  DN0M4  CNVRG 


34.0/  500. 

68.0/  1000. 
136.1/  2000. 
204.1/  3000, 
272.1/  4000. 
340.1/  5000. 


0.037  496.  0.000  0.001  0.002  0.27  267.8 

-57.4  4.35  0.15  10.503  0,000  0.000  0.000  2 

0,074  529.  0.000  0.001  0.002  0.16  133.0 

-66.9  4.31  0.15  5.290  0.000  0.000  0.000  2 


0.148  567.  { 

■77.9  4.26  0.15 


0.000  0.001 


0.001 


0.09 


66,0 


2.665 


0.000 


0.000 


0.000 


0.228  593.  0.000  0.001  0.001  0.07  42.3 

-85.6  4.22  0.15  1.726  0.000  0.000  0.000 

0.313  615.  0.000  0.001  0.001  0.06  30.2 

-91.8  4.14  0.15  1.257  0.000  0.000  0.000 

0.397  631.  0.000  0,001  0.001  0.05  23.4 

-96.7  4.06  0.15  0.992  0.000  0.000  0.000 


BINDER 


RB 

TSB 

FRB 

DELTB 

ZP 

EXZP 

QFUEL 

0.0302 

850.0 

0.99 

0.0 

0.02 

10. 

433.0 

5 

0.0309 

850.0 

0.99 

0.0 

0.00 

10. 

433,0 

6 

0.0312 

850.0 

0.99 

0.0 

0.00 

10. 

433.0 

7 

0.0334 

850 . 0 

0.99 

0.0 

0.00 

10. 

433.0 

8 

0.0337 

850.0 

0.99 

0.0 

0.00 

10. 

433.0 

8 

0.0398 

850.0 

0.99 

0.0 

0.00 

10. 

433.0 

10 

BURN  RATE  SUMMARY 


500. 

1000. 

2000. 

3000. 

4000. 

5000. 


N 

IT 

RB 

ROX 

0.00 

7 

0.0302 

0.0375 

0.99 

6 

0.0309 

0.0744 

1.02 

6 

0.0312 

0.1477 

1.08 

5 

0.0334 

0.2281 

1.08 

5 

0.0337 

0.3133 

1.04 

4 

0.0398 

0.3968 

0.340  CPU  SECONDS  USED 
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13-«1AR-89 


AN  MONOPROPELLANT 


TEST  CASE 


OXIDIZER 

'  DEPENDENT  VARIABLES: 

TYFE 

AN  1 

WT  % 
99.9 

DIA (MICRONS) 
999.00 

BDsq 

1.0002 

ARATIO 

9999.00 

ERATIO 

9999.00 

D0XNUM 

1.22 

DELTA  SFE  DEL 

0.00  1.000  0.0000 


OXIDIZER  CONSTANTS: 

TYFE —  Q'COEF - RHOX  AOX  EOX  HINAP  EAP  TMONO  AAP  KAPl 

AN  0.290  1.72  4.51E+05  12000.  2.0  25000.  1282.  127.  0.695E-02 

CIGN  TMELT  POWD  TLOSS  FRD  GEOM  HIEAP  TMONOH  HIAAP 

0.00  443.0  2.00  0.00  0.0  0.000  25000.  1282.  182. 


MISC  PARAMETERS: 

rSUB'P  XllAMB  AFH  TZERO  BPI  RHOP  BETMIN  BPMIN  BFDIV 

0.30  0.0003  3.00  347.  1.50  1.72  0.00  0.10  999999.00 

PRESSURE  DEPENDENT  RESULTS: 

“U3<  P  (ATW")7BBI  BUX  TSOX  X*PD  BETAP  BETAF  X*PF  X*AP 

QL  ZAP  ZPF  DNOMl  DN0M2  DN0M3  DN0M4  CNVRG 

AN  34.0/  500. 

1  0.043  496.  0.000  0.001  0.002  0.26  235.4 

-43.3  4.39  0.15  9.093  0.000  0.000  0.000  2 


AN 


AN 


AN 


68.0/  1000. 

1 

136.1/  2000. 

1 

204.1/  3000. 

1 


0,086  530.  0.000  0.001  0.002 

-53.2  4.35  0.15  4.581  0.000  0.000 


0.15  116.8 

0 . 000  2 


0.171  569.  0.000  0.001  0.001 

-64.5  4.30  0.15  2.309  0.000  0.000 


0.09  57.9 

0 . 000  2 


0.263  597.  0.000  0.001  0.001 

-72.5  4.26  0.15  1.495  0.000  0.000 


0.07  37.1 

0 . 000  2 


AN 


AN 


272.1/  4000. 

1 

340.1/  5000. 

1 


0.364  619.  0.000  0.001  0.001 

-79.0  4.23  0.15  1.082  0.000  0.000 


0.05  26.6 

0.000  2 


0.463  637.  0.000  0.001  0.001 

-84.2  4.21  0.15  0.850  0.000  0.000 


0.05  20.8 

0 . 000  2 


BINDER 

RB 

TSB 

FRB 

DELTB 

ZP 

EXZP 

QFUEL 

HTPB 

0.0315 

850.0 

0.99 

0.0 

0.01 

10. 

433.0 

4 

HTPB 

0.0324 

850.0 

0.99 

0.0 

0.00 

10. 

433.0 

6 

HTPB 

0.0327 

850.0 

0.99 

0.0 

0 , 00 

10. 

433.0 

7 

HTPB 

0.0351 

850.0 

0.99 

0.0 

0 . 00 

10. 

433.0 

8 

HTPB 

0.0353 

850.0 

0.99 

0.0 

0 . 00 

10. 

433,0 

8 

HTPB 

0.0420 

850.0 

0.99 

0.0 

0,00 

10. 

433.0 

10 

BURN  RATE  SUMMARY 

ROX 

RBAR 

N 

IT 

RB 

500. 

0.0433 

0 . 00 

8 

0.0315 

0.0433 

1000. 

0.0859 

0.99 

6 

0.0324 

0.0859 

2000. 

0.1705 

1.02 

6 

0.0327 

0.1705 

3000.^ 

0.2629 

1.09 

5 

0.0351 

0.2633 

4000  . 

0.3632 

1.10 

5 

0.0353 

0.3637 

5000. 

0.4607 

1.07 

4 

0,0420 

0.4633 

0. 

.300  CPU  SECONDS 

USED 
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TEMPERATURE  SENSITIVITY  DATA  BETWEEN  -65.  AND  165.  DEG  F 


PSI 

RATE 

(IN/SEC)  EXPONENT 

SIG  P  (1/C) 

PI  K  (1/F)  AVE  TS(K) 

500 

0.038 

0.00 

0.00307 

1 . 00000 

496.6 

1000 

0.074 

0.99 

0 . 00309 

0.14818 

530.7 

2000 

0.148 

1.02 

0.00317 

-0.09697 

569.6 

3000 

0.228 

1.08 

0.00323 

-0.02156 

597.3 

4000 

0.313 

1.08 

0.00328 

-0.02388 

619.6 

5000 

0.395 

1.04 

0.00332 

-0.04587 

637.4 

COMBUSTION  : 

INSTABILITY  CHARACTERISTICS  USING  THE 

DENISON  AND  BAUM  MODEL 

TDIFF 

=  0. 

000300 

P: 

=1000.  A: 

=  3.9  B= 

1.40  P=2000.  A 

=  4.1  B: 

=1.17 

P=3000. 

A: 

:  4.2 

B=1.05 

FREQ 

OMEGA 

REAL  PC 

IMAG  OMEGA 

REAL  PC 

IMAG 

OMEGA 

REAL  PC  IMAG 

100. 

34.04 

0.67 

-0.44  8.64 

1.30 

-0.50 

3.63 

1.82 

-0.18 

200. 

68.08 

0.48 

-0.36  17.29 

0.94 

-0.56 

7.27 

1.54 

-0.58 

300. 

102.12 

0.39 

-0.31  25.93 

0.75 

-0.53 

10.90 

1.26 

-0.68 

400. 

136.16 

0.33 

-0.28  34.57 

0.64 

-0.49 

14.54 

1.07 

-0.68 

500. 

170.20 

0.30 

-0.26  43.21 

0.57 

-0.45 

18.17 

0.94 

-0.66 

600. 

204.24 

0.27 

-0.24  51.86 

0.51 

-0.42 

21.81 

0.84 

-0.63 

700. 

238.28 

0.25 

-0.22  60.50 

0.47 

-0.40 

25.44 

0.77 

-0.60 

800. 

272.32 

0.24 

-0.21  69.14 

0.44 

-0.38 

29.07 

0.71 

-0.58 

900. 

306.37 

0.22 

-0.20  77.79 

0.41 

-0.36 

32.71 

0.66 

-0.55 

1000. 

340.41 

0.21 

-0.19  86.43 

0.39 

-0.34 

36.34 

0.62 

-0.53 

2000. 

680.81 

0.15 

-0.14  172.86 

0.27 

-0.25 

72.68 

0.42 

-0.39 

3000. 

1021.22 

0.12 

-0.11  259.29 

0.22 

-0.21 

109.03 

0.34 

-0.33 

4000 . 

1361.62 

0.10 

-0.10  345.72 

0.19 

-0.18 

145.37 

0.29 

-0.28 

5000. 

1702.03 

0.09 

-0.09  432.14 

0.17 

-0.16 

181.71 

0.26 

-0.25 

6000. 

2042.44 

0.09 

-0.08  518.57 

0.15 

-0.15 

218.05 

0.23 

.  -0.23 

7000. 

2382.84 

0.08 

-0.08  -  605.00 

0.14 

-0.14 

254.39 

0.22 

-0.22 

8000 . 

2723.25 

0.07 

-0.07  691.43 

0.13 

-0.13 

290.74 

0.20 

-0.20 

9000. 

3063.66 

0.07 

-0.07  777.86 

0.13 

-0.12 

327.08 

0.19 

-0.19 

INtt 
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0,  0,  1,  0 

'AN  MONOPROPELLANT  STATISTICAL  DATA  —  TEST  CASE’ 

,  ,’AN’,  ,’HTPB’ 

,  ,  ,  ,  ,  ,  1.0, 0.0, 0.0,  ,  ,0.0 

1.0, 1.0, 1.0,  1.0, 1.0, 1.0,  999. , 1 .0, 1 .0,  1.0, 1.0, 1.0 
298.0,  ,  ,  ,  ,  ,  ,  ,  ,  0.0,  0.0,  0.0 

,  298.,  985.,  .067 

,  ,  ,  ,  ,  ,  1.0  ,0.0, 0.0,  ,  ,0.0 
1.0,1.,!.,  l.,l.,l.,  999.  ,1.0, 1.0,  1.0,1. ,1. 

,  298.  ,  1050.  ,  .079 

0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  1.0  ,0.0, 0.0, 0.0, 0.0, 0.0 

1.0, 1.0, 1.0, 1.0, 1.0 


1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  999.0 

,  298.  ,  988.  ,  .083 

0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  1.0 
1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  999.0 

,  298.  ,  1500.  ,  .102 
0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 ,  1.0 
1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  999.0 

,  298.  ,  1500.  ,  .114 

0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 ,  1.0 
1 .0,1 .0, 1 .0, 1 .0, 1 .0, 1 .0,  999.0 

,  298.  ,  1750.  ,  .126 

0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  1.0 

1 .0, 1 .0, 1 .0, 1 .0, 1 .0, 1 .0,  999.0 

,  298.  ,  2000.  ,  .154 

0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  1.0 
1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  999.0 

,  298,  ,  2000.  ,  .157 

0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  1.0 
1.0, 1.0, 1.0, 1,0, 1.0, 1.0,  999.0 

,  298.  ,  2500.  ,  ,177 

0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  1.0  , 
1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  999.0 

,  298.  ,  2500.  ',  .209 

0 . 0,0 .0, 0 . 0, 0 . 0,0 . 0, 0 .0,  1.0  , 
1.0,1.0„1. 0,1. 0,1. 0,1-0,  999.0 

,  298.  ,  3000.  ,  .197 

0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 ,  1.0  , 
1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  999.0 

,  298.  ,  3000.  ,  .209 

0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  1.0  , 
1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  999.0 

,  298.  ,  3000.  ,  .213 

0.0, 0,0, 0.0, 0.0, 0.0, 0.0,  1.0  , 
1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  999.0 

,  298.  ,  3500.  ,  .256 

0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  1.0  , 
1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  999.0 

,  298.  ,  3500.  ,  .268 

0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  1.0  , 
1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  999.0 

,  298.  ,  4000.  ,  .315 

0.0,0.0,000,0.0,0.0,0.0,  1.0  , 
1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  999.0 

,  298.  ,  4000.  ,  .350 

0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  1.0  , 
1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  999.0 

,  298.  ,  4800.  ,  .374 

0.0, 0.0, 0.0, 0.0, 0,0, 0.0,  1.0  , 
1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  999.0 

,  298.  ,  4800.  ,  .437 

0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 ,  1.0  , 
1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  999.0 

’END’ ,0.0, 0.0, 0.0 
SNAMl  ISTAT=2,  ' QCOEF (3) =0 . 28 , 
SNAMl 
SNAMl 


,0.0, 0.0, 0.0, 0.0, 0.0 
,1.0, 1.0, 1.0, 1.0, 1.0 

,0.0, 0.0, 0.0, 0.0, 0.0 
,1.0, 1.0, 1.0, 1.0, 1.0 

,0.0, 0.0, 0.0, 0.0, 0.0 
,1.0, 1.0, 1.0, 1.0, 1.0 

,0.0, 0.0, 0.0, 0.0, 0.0 
,1.0, 1.0, 1.0, 1.0, 1.0 

,0.0,0.0,00  0,0.0,0.0 
,1.0, 1.0, 1.0, 1.0, 1.0 

,0.0,0.0,00  0,0.0,0.0 
,1.0, 1.0, 1.0, 1.0, 1.0 

,0.0,0.0,0.0,0.0,00  0 
,1.0, 1.0, 1.0, 1.0, 1.0 

,0.0,0.0,0.0,00  0,0.0 
,1.0, 1.0, 1.0, 1.0, 1.0 

,0.0,0.0,0.0,0.0,00  0 
,1.0, 1.0, 1.0, 1.0, 1.0 

,0.0, 0.0, 0.0, 0.0, 0.0 
,1.0, 1.0, 1,0, 1.0, 1.0 

0.0,0.0,0.0,0.0,00  0 
,1.0, 1.0, 1.0, 1.0, 1.0 

0.0,0.0,00  0,0.0,0.0 
,1.0, 1.0, 1.0, 1.0, 1.0 

0.0, 0.0, 0.0, 0.0, 0.0 
,1.0, 1.0, 1.0, 1.0, 1.0 

0.0, 0.0, 0.0, 0.0, 0.0 
,1.0,1 .0,1. 0,1. 0,1.0 

0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 
,1.0,1 .0,1 .0,1 .0,1.0 

0.0, 0.0, 0.0, 0.0, 0.0 
,1.0, 1.0, 1.0, 1.0, 1.0 

0.0,0.0,0.0,0.0,00  0 
,1.0, 1.0, 1.0, 1.0, 1.0 

SEND 

NJOB  =  3 
NJOB  =  3 


SEND 

SEND 
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13-MAR-89 


AN  MONOPROPELLANT  STATISTICAL  DATA 


TEST  CASE 


OXIDIZER  DEPENDENT  VARIABLES: 

TYPE - m"yr  — 

AN  1  99.9  999 . 00 

OXIDIZER  CONSTANTS: 

“TYPE  QCOEF  PROX  AOX  EOX  HINAP  EAP  TMONO  AAP  KAPl 

AN  0.290  1.72  4.S1E+05  12000.  2.0  25000.  1247.  127.  0.528E-02 

CIGN  TMELT  POWD  TLOSS  FRD  GEOM  HIEAP  TMONOH  HIAAP 

0.00  443.0  2.00  0.00  0.0  0.000  25000.  1247.  182. 


MISC  PARAMETERS : 


tsUBP 
0.30 
PS  I 

^LAMB 

0.0003 

TZERO 

Tfh 

3.00 

R8AR 

TZERO 

298. 

CNVRG 

BPI 

1.50 

DATA 

RHOP  BETMIN  BPMIN  BFDIV 

1.72  0.00  0.10  999999.00 

%ERROR  CNVRG  AN  DIA  WT% 

985. 

298. 

.0.073 

6 

0.067 

9.4 

2 

999.100.00 

1050. 

298. 

0.078 

4 

0.079 

-1.3 

2 

999.100.00 

988  . 

298. 

0.074 

6 

0.083 

-11.4 

2 

999.100.00 

1500. 

298. 

0.111 

5 

0.102 

8.8 

2 

999.100.00 

1500. 

298. 

0.111 

6 

0.114 

-2.5 

2 

999:100.00 

1750. 

298. 

0.129 

4 

0.126 

2.3 

2 

999.100.00 

2000. 

298. 

0.147 

4 

0.154 

-4.4 

2 

999.100.00 

2000. 

298. 

0.148 

5 

0.157 

-5.8 

2 

999.100.00 

2500. 

298. 

0.183 

4 

0.177 

3.5 

2 

999.100.00 

2500. 

298. 

0.184 

4 

0.209 

-11.8 

2 

999.100.00 

3000. 

298. 

0.227 

4 

0.197 

15.2 

2 

999.100.00 

3000. 

298. 

0.227 

4 

0.209 

8.7 

2 

999.100.00 

3000. 

298. 

0.227 

4 

0.213 

6.7 

2 

999.100.00 

3500. 

298. 

0.270 

4 

0.256 

5.4 

2 

999.100.00 

3500. 

298. 

0.271 

5 

0.268 

1 .0 

2 

999.100.00 

4000 . 

298. 

0.312 

4 

0.315 

-0.9 

2 

999 . 100 . 00 

4000 . 

298. 

0.313 

5 

0.350 

-10.7 

2 

999.100.00 

4800. 

298. 

0.378 

4 

0.374 

1-2 

2 

999.100.00 

4800. 

298. 

0.379 

5 

0.437 

-13.2 

2 

999.100.00 

SQUARE  OF  THE  CORRELATION  COEFFICIENT  FOR  THE  DATA  IS  0.95949 
STANDARD  ERROR  OF  ESTIMATE  IS  0.0045 
NUMBER  OF  DATA  POINTS  IS  19 

73%  OR  14  POINTS  ARE  WITHIN  10%  100%  OR  19  POINTS  ARE  WITHIN  20% 

MAXIMUM  DEVIATIONS  ARE  +  15.2%,  AND  -13.2% 


1.050  CPU  SECONDS  USED 
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13-MAR-89 


AN  MONOPROPELLANT  STATISTICAL  DATA 


TEST  CASE 


OXIDIZER  DEPENDENT  VARIABLES: 

TYPE - WT  «  DIACWICROrJSl 

AN  1  99  9  999  00 

BDSQ  ARATIO  ERATIO  D0XNUM  DELTA  SFE  DEL 
1.0002  9999.00  9999.00  1.22  0.00  1.000  0.0000 


OXIDIZER  CONSTANTS: 

. TYPE — - RHOX  AOX 

AN  0.280  1.72  4.51E+05 

CIGN  TMELT  POWD 
0.00  443.0  2.00 


EOX  HINAP  EAP  TMONO  AAP  KAPl 

12000.  2.0  25000.  1247.  150.  0.622E-02 

TLOSS  FRD  geom  hieap  tmonoh  hiaap 

0.00  0.0  0.000  25000.  1247.  218. 


MISC  PARAMETERS: 


rSUBT 

0.30 

PSI 

xlaMB 

0 . 0003 
TZERO 

"Xfh 

3.00 

REAR 

TZERO 

298. 

CNVRG 

BPI 

1.50 

DATA 

RHOP  BETMIN  BPMIN  BFDIV 

1.72  0.00  0.10  999999.00 

%ERROR  CNVRG  AN  DIA  WT% 

985. 

298. 

0.075 

6 

0.067 

12.4 

2 

999.100.00 

1050. 

298. 

0.080 

4 

0.079 

1.2 

2 

999.100.00 

988. 

298. 

0,076 

6 

0.083 

-9.0 

2 

999.100.00 

1500. 

298. 

0.113 

5 

0.102 

10.5 

2 

999.100.00 

1500. 

298. 

0.113 

6 

0.114 

-1.0 

2 

999.100.00 

1750. 

298. 

0.130 

4 

0.126 

3.5 

2 

999.100.00 

2000. 

298. 

0.148 

4 

0.154 

-3.6 

2 

999.100.00 

2000. 

298. 

0.149 

5 

0.157 

-5.1 

2 

999.100.00 

2500. 

298. 

0.184 

4 

0.177 

3.7 

2 

999.100.00 

2500. 

298. 

0.185 

4 

0.209 

-11.6 

2 

999.100.00 

3000, 

298. 

0.227 

4 

0.197 

15.0 

2 

999.100.00 

3000. 

298. 

0.227 

4 

0.209 

8.6 

2 

999.100.00 

3000. 

298. 

0.227 

4 

0.213 

6.5 

2 

999.100,00 

3500. 

298. 

0.269 

4 

0.256 

5.1 

2 

999.100.00 

3500. 

298. 

0.270 

5 

0.268 

0.7 

2 

999.100,00 

4000  . 

298. 

0.311 

4 

0.315 

-1.2 

2 

999.100 '00 

4000 . 

298. 

0.312 

5 

0.350 

-10.9 

2 

999.100.00 

4800. 

298. 

0.377 

4 

0.374 

0.9 

2 

999.100.00 

4800, 

298. 

0.378 

5 

0.437 

-13.5 

2 

999.100i^0 

SQUARE  OF  THE  CORRELATION  COEFFICIENT  FOR  THE  DATA  IS  0.95804 
STANDARD  ERROR  OF  ESTIMATE  IS  0.0046 
NUMBER  OF  DATA  POINTS  IS  19 

685!  OR  13  POINTS  ARE  WITHIN  10%  100%  OR  19  POINTS  ARE  WITHIN  20% 

MAXIMUM  DEVIATIONS  ARE  +  15.0%,  AND  -13.5% 
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0,  0,0,1 

’AN/AP/SN  MIXED  OXIDIZER  TEST  CASE' 

'AP',  , 'AN' , 'SN' , 'HTPB' 

0.0001  ,0.0, 0.0,  0.0, .0  , .0  ,  0.39,0.39,0.0,  0.0001,0.0,0.0 

5. ,1.0, 1.0,  20. ,1.0, 1.0,  50.,  50. ,1.0,  200. ,1.0, 1.0 

298.0,  0.0,  ,0.0  ,0.0,  0.0,  ,  0.12,  ,  0.0,  0.0,  0.0 

SNAMl  ALF(3 , 1)=0.4,  ALF  (3 ,2)  =0 .4 ,  DZERO  (3 , 1")  =150  .  ,  CATCON  (3)  =0 . 02  ,  SEND 
SNAMl  ALF(3,1)=0.32,  ALF (3 , 2) =0 . 32 /  ALF (1 , 1)  =0 . 1 , ALF (4 , 1) =0 . 1 ,  SEND 

SNAMl  NJ08  =  3  SEND 

SNAMl  NJOB  =  3  SEND 

000000000111111111122222222223333333333444444444455555555556666666666777 

123456789012345678901234567890123456789012345678901234567890123456789012 
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AN/AP/SN  MIXED  OXIDIZER  TEST  CASE 


OXIDIZER  DEPENDENT  VARIABLES: 


T7FE 

WT"}r 

DIA(WiCkuNs) 

AP  1 

0,0 

5.00 

BDSQ 

ARATIO 

ERATIO 

D0XNUM 

DELTA 

SFE 

DEL 

8.2904 

0.30 

5.70 

0.939E-04 

0.90 

0.002 

0.2205 

AN  1 

40.0 

150.00 

BDSQ 

ARATIO 

ERATIO 

D0XNUM 

DELTA 

SFE 

DEL 

1.2143 

8.92 

2.83 

0.426 

6.94 

0.785 

0.05P  '  ^ 

2 

40.0 

50.00 

1.6430 

2.97 

10.40 

0.426 

4.01 

0.214 

0.0  '98, 

SN  1 

0.0 

200 . 00 

BDSQ 

ARATIO 

ERATIO 

D0XNUM 

DELTA 

SFE 

D'  EL 

1.2113 

11.89 

5.20 

0.810E-04 

9.19 

0 . 000 

0.'  0563 

OXIDIZER  CONSTANTS: 

TYPE  '  ■ 

QCOET 

fTHOX 

AOX 

EOX  HINAP  EAP  TMONO 

AAP  KAPl 

AP 

0 . 200 

1.95  3 . 

22E-k08 

30000 . 

1.7  15000.  1413 

.  365, 

.  1 

,  .75 

CIGN  TMELT 

POWD 

TLOSS 

FRD 

GEOM 

HIEAP 

TMONOH 

HI.-  AAP 

4  , 

.32  865.0 

1.70 

1.00 

900.0 

0.200 

15000. 

1413  . 

12  i  . 

TYPE 

QCOEF 

RHOX 

AOX 

EOX  HINAP  EAP  TMONO 

AAP  KAPl 

AN 

0.290 

1.72  4. 

51E+05 

12000. 

2.0  25000.  1247 

.  127. 

0. 

528  •E-02 

CIGN  TMELT 

POWD 

TLOSS 

FRD 

GEOM 

HIEAP 

TMONOH 

HT  -AAP 

0. 

,50  443.0 

1.70 

1.00 

100.0 

0.300 

25000. 

1247. 

ir82. 

TYPE 

QCOEF 

RHOX 

AOX 

EOX  HINAP  EAP  TMONO 

AAP  KAP  1 

SN 

0.260 

2.26  1. 

15E+02 

4500. 

2.0  25000.  999 

.  800. 

0. 

27  lE-02 

CIGN  TMELT 

POWD 

TLOSS 

FRD 

GEOM 

HIEAP 

TMONOH 

H  lAAP 

0. 

10  580.0 

1.70 

1.00 

100.0 

0.300 

25000. 

999. 

9  7.3 

HTPB  BINDER  CONSTANTS: 

QFUcL 

. PhOF" 

IF 

EF 

CSUBPB 

ALFF 

433.0 

0.90 

6 .74E+02 

10200 . 

0.43 

0.180 

CATALYST  PROPERTIES 
-ITR - - EOT 

ACAT  CATK 

AAA 

CATSSA 

XNCl 

AN 

0.020  10000 

.  183.6 

0 , 1000 

0 . 1200 

1.00 

MISC 

rEDBP 

PARAMETERS : 
XLAK^B  ?fFH 

1)  0.22 

2)  1.00 

TZERO  BPI  RHOP 

BETMIN 

BPMIN 

BFDIV 

0.30 

0.0003  3.00 

298.  1.50  1.49 

0.30 

0.10 

0.10 

PRIMARY  FLAME  CONSTANTS 
TPP  (K) 


QPF  KPF  TPF(K)  XNFT 

261.  0.784  1169.  2.00 

273.  0.507E-01  1266.  1.50 

491.  0.232  2050.  1.50 

626.  0.327  2385.  1.50 

PRESSURE  DEPENDENT  RESULTS: 


APF 

500. 

2.70 


2 

2 


70 

70 


EPF 
15000. 
10000 . 
10000 . 
10000 . 


THT 


AP 


AN 


SN 


AP 


P(ATM)/PSr 

34.0/  500 


P!UX  TSOX  X*PD  BETAP 
QL  ZAP  ZPF  DNOMl  DN0M2 


68.0/  1000 


0.130  897.  1.446  0.542 

-119.8  5.99  0.93  3.032  0.213 

0.036  493.  70.075  0.658 

-56.7  5.78  2.72  9.862  0.237 

0.137  562.  7.108  0.499 

-76.6  5.78  5.56  2.872  0.124 

0.057  831.  65.694  0.542 

-138.6  5.99  2.92  6.907  0.065 


0.286 


943  . 


1.446  0.542 


X«PD 

STOIC 

FROXM 

1.45 

5.70 

1.00 

70.07 

10.40 

1.00 

7.11 

10.40 

1.00 

65.69 

5.20 

1.00 

BETAF 

X*PF 

X^i^AP 

DNOM3  DNOM4  CNVRG 

0 

.999 

11.70 

2050.3 

3. 

773  267. 

284 

2 

0 

.950 

103.45 

497.0 

0. 

971  0. 

991 

2 

0 

.950 

71.41 

1907.7 

1  . 

681  2. 

970 

2 

0 

.999 

54.29 

1043.5 

0. 

964  0. 

742 

2 

0 

.999. 

4.93 

1128.3 

102 


AN 


SN 


AP 

AN 


SN 


-129.0  5.99  0.91  1.377  0.126  2.593  125.983  '  2 

1  0.059  517.  70.075  0.658  0.950  62.65  204.2 

-63.4  5.05  2.91  5.773  0.134  0.667  0.467  2 

2  0.213  589.  7.108  0.499  0.950  44.71  741.4 

-84.4  5.78  5.42  1.848  0.073  1.155  1.400  2 

1  0.082  1026,  65.694  0.542  0.999  26.42  374.8 

-189.3  5.99  2.62  4,808  0.039  0.662  0.350  2 

136.1/  2000. 

1  0.600  991.  1.446  0.542  0.999  2.05  592.0 

-138.7  5.99  0.88  0.656  0.079  1.871  65.854  2 

1  0.096  542.  70.075  0,658  0.950  36.25  83.7 

-70.8  3.43  3.16  3.577  0.085  0.481  0.244  2 

2  0.337  620.  7.108  0.499  0.950  27.75  293.0 

-93.3  5.77  5.38  .  1.169  0,046  0.834  0.732  2 

1  0.107  1259.  65.694  0.542  0.999  11.84  122.7 

-249.8  5.99  2.36  3.672  0.024  0.478  0.183  2 


BINDER 

RB 

TSB 

FRB 

DELTB 

ZP 

EXZP 

QFUEL 

HTPB 

0.0460 

922. 

9 

0.05 

0.0 

0.02 

8. 

433.0 

4 

0.81 

51.4 

1.07 

391 . 

433.0 

0.27 

17.9 

0.16 

1158. 

433.0 

0.99 

0.0 

0.28 

0. 

433.0 

HTPB 

0.0669 

1001 . 

3 

0.05 

0.0 

0.01 

8. 

433.0 

4 

0,81 

50.9 

0.92 

464  . 

433.0 

0.27 

17.3 

0.14 

1164  . 

433.0 

0.99 

0.0 

0.28 

0. 

433.0 

HTPB 

0.0928 

,1085. 

6 

0.05 

0.0 

0.01 

9. 

433.0 

3 

0.81 

47.3 

0.76 

541  . 

433.0 

0.27 

15.4 

0.11 

1182. 

433.0 

0.99 

0.0 

0.31 

0. 

433.0 

BURN  RATE  SUMMARY 

^T7^ 

RBaR 

N 

IT 

RB 

ROX 

500. 

0.0461 

0.00 

8 

0.0460 

0.1299 

0.0357  0.1371  0 

.0570 

1000. 

0.0789 

0.73 

7 

0.0669 

0.2859 

0.0587  0.2131  0 

.0819 

2000 . 

0.1267 

0.68 

7 

0.0928 

0.5999 

0.0963  0.3367  0 

.1072 

0.790  CPU  SECONDS  USED 
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OXIDIZER  DEPENDENT  VARIABLES: 
Tr%  biAm<:R0N3) — 


AP 

1 

10.0 

5.00 

F. 

AN 

1 

32.0 

150.00 

BDSQ 

3.0896 

ARATIO 

1.04 

2 

32.0 

50.00 

BDSQ 

1.0614 

ARATIO 

31.11 

SN 

1 

10.0 

200.00 

1.1843 

10.37 

BDSQ 

1.0606 

ARATIO 

41.48 

OXIDIZER  CONSTANTS: 

"i  XPE  QLoep  RT70X 


an/ap/sn  mixed  oxidizer  test  case 


1-04  5.70  0.100 

*»RATIO  ERATIO  60XNUM 
31.11  3.58  0.363 

10.37  10.40  0.363 

\RATIO  ERATIO  D0XNUM 
41.48  5.20  0.864E-01 

EOX  HINAP  EAP  TMONO 


delta 

0.48 

SFE 

0.441 

DEL 

0.1180 

DELTA 

3.72 

SFE 

0.408 

DEL 

0.0304 

2.15 

0.141 

0.0526 

DELTA 

4.92 

SFE 

0.011 

DEL 

0.030; 

AP  0.2ee  1.9S  3.22E.08  35^0.  ij  15000  U13  365*"  f*^J 

^32  sste  90^0  fZ  «aEa1"-tmonoh''h1**p 

TYPE  QCOEF  RHOX  AOX  ^  EOX^^HINAP  PAP 

AN  0.290  1.72  4.61E.05  12000.  T0  2^000  1,^*" 

IaI^j  h^"!a?= 

0’^0°l0"  2"”^6^  1.^s“e\02  ^^0^0.  "2^^^";i0*^0  "‘-'“/‘PC; 

0.?0  H||AP  TMONOH  HIAAp" 

HTPB  binder  CONSTANTS : 

^  CSUBPB  ALFF 

433.0  0.90  6.74E+02  10200.  0.43  0.140 

CATALYST  PROPERTIES : 

9^  CONC  ECAT  ACAT  CATK  AAA  TATCCA  vkin 
AN  0  020  10PiPiC^  1 R  n  CATSSA  XNCl 

<o.a^io  10000.  183.6  0.1000  0.1200  1.00 

1)  0 . 24 

2)  1.00 

MISC  PARAMETERS: 

^5°^ 


AAA  CATSSA  XNCl 

0.1000  0.1200  1.00 


0.30 


BFDIV 

0.10 


x*pd  stoic 

1.56  5.70 

97.81  10.40 

7.11  10.40 

64.35  5.20 


FROXM 

1.00 

1.00 

1.00 

1.00 


npc  ,^r.r-  P^^'^ARY  flame  CONSTANTS: 

26P.  0.9^"  ’TIIFTIS  7^  J:L  ^TOIC  FROXM 

288.  0.S68E-01  1303  1  50  I  70  Jiff!'  5.70  1.00 

A99.  0.232  205!:  I.TI  \i%  1-00 

616.  0.327  2385.  1  50  97C5  7.11  10.40  1.00 

PRESSURE  DEPENDENT  RESULTS-'  ’  10000.  64.35  5.20  1.00 

PT^im)/psi  ^x  tsox  X.PD  BETAP  BETAF  X.PF  X*ap 
AP  34.0/  500.  DN0M3  DN0M4  CNVRG 

AN  ^  -118^5^^^5  72^^P(  or  ^  ‘  0.950  7.81  1825.4 

iitJ.b  5.72  0.95  3.405  0.138  0  695  <5  R44  o 

0.035  493.  97.814  0.631  0.950’  88.54’  492  9 

2  Vl28  Iss  7  0.035  2 

SN  -75^/^%  78^®i  B7^'^?®a7a®'‘'f®  *5.60  1779.4 

1  ®  3.079  0.080  0.310  0  104  o 

^  64.354  0.515  0.950  50.24  1015  3  ^ 

AP  68.0/  1000. ^•“*5  5-5A2  0.178  ^02^  -^  2 

1  0.251  935.  1.555  0.515  pi  qr«  ^  op,  o.,„  _ 


68.0/  1000. 


1.555  0.515  0.950 


3.28  992.5 


104 


AN 


SN 


AP 


1 

2 

1 


136.1/  2000. 


-127.5  5.71  0.96  1.566  0.082  0.445  3.780  2 

0.058  516.  97.815  0.631  0.950  54.18  200.5 

-63.2  4.84  2.70  3.185  0.047  0.115  0.014  2 

0.198  584.  7.108  0.499  0.950  41.59  689.6 

-83.0  5.78  4.72  1.987  0.047  0.198  0.042  2 

0.079  1003.  64.354  0.515  0.950  24.27  362.0 

-183.3  5.71  2.58  4.977  0.025  0.114  0.011  2 


1 

0.514  981 

1.555  0.515 

0,950 

1. 

32  507.4 

AN 

-136.6 

5.72 

0.95  0.766  0.051 

0. 

322 

1.977 

1 

0.085  536 

.  97.815  0.631 

0.950 

28. 

53  73.8 

-68.9 

2.67 

2.67  2.025  0.028 

0. 

083 

0,007 

2 

0.315  615 

7.108  0.499 

0,950 

25. 

97  274.2 

SN 

-91.9 

5.78 

4.74  1.249  0.030 

0. 

144 

0.022 

1 

0.104  1230 

.  64.354  0.515 

0.950 

10. 

97  119.5 

-242.4 

5.71 

2.26  3.770  0.016 

0. 

082 

0.005 

BINDER 

R8 

TSB  FRB 

DELTB 

ZP 

EXZP 

QFUEL 

HTPB 

0.1337 

1194.3  0.17 

25.4 

0.22 

1604  . 

433.0  3 

0.99 

22.4 

8.00 

0. 

433.0 

0.95 

10.2 

1.34 

237. 

433.0 

0.99 

0.5 

2.46 

4 . 

433.0 

HTPB 

0.2086 

1362.6  0.17 

26.5 

0.12 

1778. 

433.0  4 

0.99 

22.9 

8,00 

0. 

433.0 

0.95 

10.4 

1.36 

230. 

433.0 

0.99 

0.6 

2.76 

3. 

433.0 

HTPB 

0.2883 

1520.0  0.17 

27.0 

0.07 

1874. 

433.0  2 

0.99 

22.3 

8.00 

0. 

433.0 

0.95 

10.1 

1.11 

293  . 

433.0 

0.99 

0.6 

3.06 

2. 

433.0 

BURN  RATE  SUMMARY 

RBAR 

N  IT 

RB 

ROX 

500. 

0.0756 

0.00  6 

0.1337 

0.1156  0.0354  0.1279  0 

.0555 

1000. 

0.1281 

0.71  6 

0.2086 

0.2515  0.0576  0.1982  0 

.0791 

2000. 

0.2023 

0.66  6 

0.2883 

0.5142  0.0849  0.3151  0 

.1044 

0.690  CPU  SECONDS  USED 
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0,  0,  0, 

’AN/HTPB 


0 


. ’AN’, 


0.0  ,0.0, 0.0, 
5  .  ,  1 .0, 1 .0, 
298.0,  0.0,  ,0 
SNAMl 
SNAMl 
SNAMl 
SNAMl 


TEST  CASE’ 

, ’HTPB’ 

0 . 0  ,  .  0  ,  .  0  , 

20. ,1.0, 1.0, 

0  ,0.0,  0.0, 
CATCON(3)=0.01, 
ATCON(3)=0.02, 


0.39,0.39,0.0,  0.0, 0.0, 0.0 

150. ,50. ,1 .0,  200. ,1.0, 1.0 
0.12,  ,  0.0,  0.0,  0.0 
ALF(3,1)=0.385,  ALF (3 , 2) =0 . 385 , 
ALF(3,1)=0.38,  ALF(3,2)=0.38 

NJOB  =  3 
NJOB  =  3 


SEND 

SEND 


000000000111111111122222222223333333333444444444455555555556666666666777 

123456789012345678901234567890123456789012345678901234567890123456789012 
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AN/HTPB 


rEST  CASE 


U  JmW  A. 

T7P't“ 

WT  7.  “ 

OIA (MICRONS) 

AN  1 

39.0 

150.00 

2 

39.0 

50.00 

TZERO=' 

298. 

BURN  RATE  SUMMARY 

IT 

“FTA 

REAR 

N 

500. 

0.0263 

0.00 

9  0 

1000. 

0.0456 

0.82 

.9  0 

2000. 

0.0814 

0.84 

8  0 

RB 


ROX 

0.0270 

0.0497 

0.0888 


13-MAR-89 


AN/HTPB 


0.1108 

0.1862 

0.3100 

0.540  CPU  SECONDS  USED 
TEST  CASE 


TYFI - 

-wnr~\ 

DIA (MICRONS) 

AN  1 

38.5 

150.00 

2 

38.5 

50.00 

TZER0=298. 

BURN  RATE  SUMMARY 

PIA 

RBAR  ■ 

N 

IT 

RB 

500. 

0.0285 

0 . 00 

9 

0.0299 

1000. 

0 , 0505 

0.80 

8 

0.0428 

2000. 

0.0859 

0.77 

8 

0.0595 

13-MAR 

-89 

AN/HTPB 

OXIDIZER 

DEPENDENT  VARIABLES: 

ROX 

0.0322 

0.0559 

0.0970 


AN 


1 

2 


38.0 

38.0 


150.00 

50.00 


0.1289 
0.2081 
0.3357 

0.370  CPU  SECONDS  USED 

TEST  CASE 


TZER0=298  . 

BURN  RATE  SUMMARY 

PTA” 


500.  0.0307 
1000.  0.0539 
2000.  0.0890 


N 

0.00 

0.77 

0.72 


IT 

9 

8 

8 


RB 

0.0307 

0.0441 

0.0598 


ROX 

0.0364 

0.0614 

0.1045 


0.1442 

0.2272 

0.3588 

0.430  CPU  SECONDS  USED 
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APPENDIX  B 


Composite  Propellant  Modeling  Bibliography 
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Composite  Propellant  Modeling  Bibliography 

1988 

Deur,  J.  M.  and  Price^.W.,  "A  Surface  Coupled  Flamelet  Approach  to  Dynamic  Response  in 
Heterogeneous  Propellant  Combustion",  AIAA-88-2938,  (1988), 

NOTE:  Apply  the  KLLEF  model  to  unstable  combustion.  A  very  unrealistic  physical  model. 
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Beckstead,  M.  W. ,  Derr,  R.  L.  and  Price,  C.  F.,  "The  Combustion  of  Solid  Monopropellants 
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Describing  Combustion  of  Solid  Composite  Propellants  Containing  Ammonium  Nitrate  , 
Combustion  and  Flame,  Vol.  3,  (1959),  pp.  301-318. 


112 


APPENDIX  C 


Aluminum  Combustion  References 


113 


Aluminum  Combustion  References 


1987 
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Ivanov,  G.V  Surkov,  V.G.,  Karmadonov.  LN. ,  and  Viktorenko,  A.m.,  ’The  Role  of  Condensed-Phase 
ShoTSi  vorarNo'  “  CP-bustio.,  Explosioa  S 

low in“ .hen lorm 

voi.  irTj  Qo  D °{ Metaiized  Propellants”,  Fundamentals  of  Solid-Propellant  Combustion, 
Volume  90,  Progress  in  Astronautics  and  Aeronautics,  Chapter  14,  (1984)  po  479-514 
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CombSrAS^'^' '22*?fl  wi.'Sp  fSTi  sf  •  Solid-Propellant 

.  NOTE:  Data  on  both  size  and  fraction  of  A1  agglomerated.  Agglomerate  size  decreases  with  pressure  and 

wi  h  increasing  fine  oxidizer  size,  and  with  increasing  fine  fraction.  Agglomerate  size  is  maximum  for  -15-2^ 
micron  Al.  Price  ignores  the  Aerojet  work  is  ludximum  ror  io  2J 


(1983^°pp  7^o’725  ^  Model  for  Aluminum  Agglomeration  in  Composite  Propellants",  AIAA  J.,  21 , 5, 

,4-  0^  al  agglomerate  size  and  concentration  based  on  pocket  size,  melting  of  AL  (invorselv 

of  1400Kr  temperature  to  an  fgnition  temperaLe 

Frolov  Yu.  V..  and  Nikolshii,B.E.,  "Combustion  Characteristics  of  Metallized  Compositions”  Combusti^  ' 
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is  probably  domln^™^  experiments;  Conclusion  is  there  are  many  contributing  factors  but  oxide  melting 


-n  &^akov,V.A..  and  Razdobreev,A.A.,  "Dispersity  of  Condensed  Products  of  Combustion  of 

an  Alum.num  Drop  ,  Combustion,  Explosion  S  ShockWaves,  18,  2,  (1982),  pp  16-19. 

NOTE;  Measured  size  distribution  by  dropping  burning  particles  onto  a  glass  plate  in  air  at  atmospheric 
pressure.  Tney  ooservered  1)  most  oxide  particles  are  1-2  microns.  2)  As  Do  increases  the  maximum  o^xide  size 
also  increases,  3)  Most  of  the  oxide  is  in  the  form  of  hollow  spheres.  4)  flame  location  is  at  a  radius  of  ■-3Do. 

’/?coo^  Aluminum  Level  on  Solid  Rocket  Motor  Performance”,  19th  JANNAF  Combustion 

Meeting,  ii,  ( iS82),  pp  17-26. 

motor  Compares  SPP  calculations  for  HTPB  and  NEPE  propellants  containing  AP,  HMX,  and  Al  with 
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Renie,  J.P.,  Lilley,  J.S.,  Frederick,  R.A. ,  and  Osborn,  J.R.,  "Aluminum  Particle  Combustion  in  Composite 
Solid  Propellants",  AIAA-82-1 110, 18th  Joint  Propulsion  Conference,  (1982), 

NOTE:  Describes  current  agglomeration  models.  Describes  PEM  Al  combustion  model:  no  agglomeration 
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Price,  E.W.,  Park,  C.J.,  Sigman,  R.K.,  and  Sambamurthi,  J.K.,  'The  Nature  and  Combustion  of 
Agglomerates",  18th  JANNAF  Combustion  Meeting,  III,  (1981),  pp  121-145. 
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